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ABSTRACT 


The  transition  matrix,  relating  the  scattered  and  incident  acoustic  waves 
for  a  thin  elastic  spherical  shell  in  a  free-fleld  environment,  has  been 
evaluated  using  the  ATILA  finite-element  code.  A  three-dimensional 
finite-element  model  of  a  0.5-m  outer  radius,  1-cm  thick  spherical  steel  shell 
surrounded  by  water  was  developed.  The  ATILA  code  was  used  to  calculate 
the  scattered  pressure  over  the  surface  of  the  shell  for  incident  waves 
represented  as  products  of  radial  Hankel  functions  and  spherical  harmonics. 
The  chosen  driving  frequency  was  474  Hz,  corresponding  to  a  value  of  ka=l, 
where  k  is  the  wavenumber  of  sound  in  water  and  a  is  the  radius  of  the  shell. 
The  ATILA  results  were  compared  with  the  results  of  analytical  thin  shell 
theory,  and  were  found  to  agree  for  a  model  which  divided  the  spherical 
shell  surface  into  72  approximately  equal  area  triangular  regions.  Also 
computed  for  each  component  was  the  modal  acoustical  impedance  of  the 
shell.  These  results  agreed  within  two  percent  for  the  zeroth  order 
component  and  thirteen  percent  for  the  first  order  components. 
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I.  INTRODUCTION 


The  utilization  of  sonar  systems  in  operations  at  sea  depends  on  the 
ability  to  forecast  their  performance  during  the  design  and  production  phases 
of  their  development.  For  low  frequency  active  sonar  systems,  arrays  of 
large  transducers  are  required  to  produce  the  necessary  power  output. 
Because  of  their  size  and  power  output,  interactions  among  the  elements  of 
an  array  may  need  to  be  considered  in  modelling  a  given  array. 

Professors  S.  R.  Baker  of  the  Physics  Department  and  C.  L.  Scandrett 
of  the  Mathematics  Department  of  the  Naval  Postgraduate  School  are 
engaged  in  a  joint  research  effort  with  the  objective  of  developing  the 
models  necessary  to  predict  the  performance  of  arbitrary  dense  active  sonar 
arrays. 

The  approach  used  in  this  research  is  based  on  the  T-Matrix  method. 
[Ref.  3]. 

The  T-matrix  method  uses  the  superposition  of  spherical  harmonics  to 
represent  the  total  radiated  pressure  from  a  transducer.  The  radiated 
pressure  from  one  such  transducer  for  each  incoming  harmonic  will  be 
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calculated  in  an  effort  to  produce  the  transition  matrix  or  T-matrix.  The 
elements  of  the  T-matrix  may  then  used  to  calculate  scattered  pressures 
resulting  from  the  interaction  between  transducers. 

The  T-matrix  method  is  useful  in  computing  the  scattered  pressure 
resulting  from  an  arbitrary  incident  pressure  on  a  given  scatterer.  This  thesis 
is  concerned  with  the  application  of  the  T-matrix  method  on  an  array  of 
thin-shelled  elastic  spheres. 

Once  the  theoretical  results  from  scattering  of  spherical  harmonics 
incident  on  a  sphere  have  been  obtained,  the  calculations  from  the  T-matrix 
method  can  be  compared  for  each  harmonic  with  the  computation  from  the 
Finite  Element  Code  ATILA  for  the  purpose  of  validation. 

Furthermore,  the  acoustic  impedance  calculated  from  the  thin  shell 
theory  (Chapter  II)  can  be  compared  with  the  acoustic  impedance  computed 
by  ATILA's  results. 

The  remainder  of  this  thesis  is  divided  into  five  chapters.  Chapter  II 
describes  the  theory  involved  in  the  T-Matrix  method,  the  Pritchard 
Approximation,  the  finite  element  analysis  of  a  structure  excited  by  an 
impinging  wave,  and  the  theory  of  forced  vibration  of  a  spherical  shell. 

Chapter  III  describes  the  spherical  model  used  in  the  ATILA  code.  Chapter 

2 


IV  presents  and  discusses  the  results.  Chapter  V  presents  the  conclusions. 
Appendix  A  contains  a  copy  of  the  FORTRAN  program  used  to  generate  the 
spherical  harmonic  impinging  wave.  Appendix  B  presents  a  copy  of  the 
code  used  in  ATILA.  Appendix  C  presents  a  FORTRAN  code  [Ref.  17]  to 
calculate  elements  of  the  Transition  Matrix. 
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II.  THEORY 


A.  "T-MATRIX"  METHOD 

The  "T-Matrix"  method  is  a  procedure  for  computing  the  acoustic  field 
due  to  multiple  radiating  and/or  scattering  bodies.  An  outline  is  given  below 
of  how  this  method  may  be  applied  to  compute  the  acoustic  pressure  due  to  a 
pair  of  piezoelectric  spherical  shell  transducers. 

If  a  voltage  is  applied  to  a  thin  piezoelectric  spherical  shell  in  a 
free-field  environment,  a  displacement  of  the  surface  of  the  shell  can  be 
obtained  by  the  canonical  equations  [Ref.  1],  This  movement  generates  an 
acoustic  pressure  field  in  the  surrounding  medium.  This  is  indicated  in 
Figure  1  ,  where  a  sphere  excited  by  a  voltage  V  produces  a  radial  velocity  u 
over  the  surface  of  the  sphere.  The  perturbation  of  the  sphere  on  the 
environment  yields  a  radiated  pressure  pR(r,0,(j)). 

The  acoustic  pressure  at  a  field  point  z ,  located  at  spherical  coordinates 
r,  0,  and  equals  pR(r,0,<J>),  since  the  sphere  is  in  a  free-field  environment. 
Numerical  values  of  pR(r,0,0)  can  be  found  by  application  of  a  finite-element 
analysis  code  such  as  ATILA  [Ref.  2]. 


4 


Figure  1 .  Single  Radiation  from  a  Plezolectric  Sphere. 


Next,  consider  the  case  of  an  array  of  two  thin  piezoelectric  spherical 
shell  transducers  that  are  close  enough  to  be  considered  "in  the  near  field"  of 
one  another  (see  Figure  2). 
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Outgoing  waves  at  sphere  1  are  composed  of  a  radiated  wave,  pf, 
resulting  from  the  applied  voltage  V,  and  a  scattered  wave,  pj,  due  to  the 

incident  wave  from  sphere  2,  p1, .  The  scattered  pressure  field  can  be 
mathematically  represented  by  outgoing  spherical  Hankel  functions  and 

surface  harmonics  in  the  sphere  1  coordinate  system  as  [Ref.  3] 

p?(rI,el,(j),)  =  I^S^P?nX)(krOPir(cos0,)ei^'  (1) 

where  hi^kri)  =  jn(krO  -  iyn(kri )  is  the  Hankel  function  of  the  second 
kind  (note  that  an  e1<ot  harmonic  time  dependence  has  been  assumed;  for  e*j<1* 
harmonic  time  dependence  hit )(krO  =jn(kri)  +  iyn(kri)  must  be  used  instead 
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of  h<2),  P?nm  is  the  amplitude  of  the  n,  mth  scattered  wave  component,  k  is 
the  radial  wave  number,  P^cosGi)  is  the  associated  Legendre  function,  and 


r,  ,  0,  and  0,  are  the  spherical  coordinates  of  the  field  point  z  with  origin  at 
center  of  sphere  1 .  The  radiated  pressure  field  is  similarly  represented: 

pf(r  1 ,  e  Zms-n  Pfnmhn2)(kri  )Pftcos  6 ,  (2) 

The  incident  pressure  pi  is  assumed  to  be  of  the  form: 

pi  (r  1 , 0 1 ,  <|>  | )  =  Zio  £m=-n  Pinm  jn(kr  1  )PJT(COS  0 1  )e in*'  (3) 

This  form  is  valid  provided  r,  is  less  than  the  distance  between  the 
origins  of  spheres  1  and  2. 

By  definition,  the  coeficients  of  the  scattered  pressure  Pfnm  are  related 
to  the  coefficients  of  the  incident  pressure  Pinm  by 

(Pf}  =  [T]{Pi  }  (4) 

Where  {}  denotes  a  column  vector  of  coefficients  and  [T]  is  the 
so-called  "transition  matrix"  or  T-Matrix ,  a  property  of  the  transducer. 

Similar  expressions  to  the  above  apply  to  the  corresponding  fields  on 
sphere  2,  i.e., 

pf(r2,02,<t>2)  =  X  ^PL,hLl,(kr2)PH'(cose2)e1"^  (5) 

{P2}  =  [T]{P2}  .  (6) 
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Note  that  it  is  assumed  here  for  the  sake  of  simplicity  that  the  two 
transducers  are  identical,  so  that  their  T-matrices  are  also;  in  the  general 
case  this  need  not  be  so. 

To  proceed  fruther  it  is  necessary  to  be  able  to  transform  from  one 
coordinate  system  into  another.  This  is  done  using  an  "addition  theorem", 
by  which  a  spherical  wave  relative  to  one  origin  is  expressed  as  a  series  of 
spherical  waves  relative  to  another  [Ref.  3  and  references  therein].  Thus,  for 
example,  the  incident  pressure  on  sphere  1  can  be  expressed  as  the  sum  of 
the  radiated  and  scattered  pressures  from  sphere  2,  transformed  into  the 
sphere  1  coordinate  system  using  the  addition  theorem.  In  matrix  form,  this 
can  be  written 

{P|}  =  [G2i]{P2+Pf}.  (7) 

where  [G21]  represents  the  matrix  of  coefficients  that  transforms  outgoing 
spherical  harmonics  at  sphere  2  into  standing  wave  spherical  harmonics  at 
sphere  1.  Similarly,  the  incident  pressure  on  sphere  2  is  related  to  the 
radiated  and  scattered  pressure  on  sphere  1  by 

{Pi}  =  [G,2]{Pf  +  P?}.  (8) 
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Using  Equations  (4)  and  (6)  in  Equations  (7)  and  (8),  we  have  a  system 
of  two  equations  in  two  unknowns: 


{P?}  =  [T][Gj,]{Pjtt  +  Pf}, 


(9) 


{Pf}  =  [T][G,2]{Pf  +  P?}. 

Equations  (9)  and  (10)  may  be  written 

I  -TG2,  Ifpfl  r  0  TG2,  lfpf  ( 
_-TG,2  I  J  [  P|  J  |_TGi2  0  J|p£ 


(10) 


(11) 


where  [I]  is  the  identity  matrix. 

If  {Pf },  {?2 },  and  [T]  can  be  found,  such  as  by  using  the 
finite-element  code  ATILA,  then,  since  [G,2]  and  [G21]  are  known  from  the 

positions  and  orientations  of  the  transducers  through  application  of  the 
addition  theorem,  the  scattered  pressures  {P?}  and  {P2}  can  be  determined. 
Then  the  resulting  pressure  at  the  field  point  z,  p(z),  is  given  by 

p(z)  =  pf(zi)  +  p?(zi)  +  p2(z2)  +  pf(z2)  ,  (12) 


with  the  forms  for  the  radiated  and  scattered  pressure  fields  given 
previously,  where  zx  represents  the  coordinates  of  the  field  point  z  relative  to 
the  origin  of  the  ith  sphere. 


The  above  development  may  be  generalized  to  an  array  composed  of  an 
arbitrary  number  N  of  identical  elements.  For  the  ith  element, 

{P?}  =  [T]{Pj}.  (13) 

Here  the  scattered  pressure  from  the  ith  transducer  is  represented.  The 
components  of  the  radiated  pressure,  may  be  found  using  such  a: 
ATILA  finite-element  code.  Then  the  acoustic  pressure  at  a  field  point  z  is 

the  sum  of  the  radiated  and  scattered  pressures  of  each  transducer. 

p(z)  =  I  £i[pf(z«)+P?(Zi)]  .  (14) 


B.  THE  ’’PRITCHARD  APPROXIMATION” 

The  analytical  calculation  of  mutual-radiation  impedance,  Zij5  is  made 
through  the  solution  of  the  wave  equation  for  the  sound  pressure  produced 
by  one  transducer  and  by  integrating  that  pressure  over  the  radiating  surface 
of  another  transducer  [Ref.  4]. 

In  Pritchard's  original  paper  [Ref.  5],  the  mutual-radiation  impedance 
was  approximated  for  the  case  of  flat  circular  pistons  on  an  infinite  rigid 
baffle.  For  small  pistons  and  large  separation,  the  Pritchard  equation  can  be 
reduced  to 
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Z,2  =  i(ka)^  +  i^)pcA. 


(15) 


Here  an  e"“  harmonic  time  dependance  is  assumed.  This 
approximation  neglects  secondary  scattered  pressures. 

It  can  be  noted  that  Z,2  has  a  real  and  an  imaginary  part, 

Z12  =  R12 +jXl2  .  (16) 

Here,  Z,2  is  the  mutual-radiation  impedance  between  transducers  1  and 
2,  d  is  the  distance  between  the  centers  of  two  transducers  (see  Figure  3),  a 
is  the  piston  radius,  k  is  the  acoustic  wave  number,  pc  is  characteristic 
impedance  of  the  environment,  and  A  is  the  area  of  the  pistons. 

Piston  Transducer  1  Piston  Transducer  2 


Figure  3.  The  Circular  Pistons  derived  by  Pritchard. 


The  self  radiation  resistance  R^,  for  small  ka,  can  be  approximated  by 
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R«  =  j(ka)2pcA.  (17) 

The  self  radiation  resistance  is  the  real  part  of  the  transducer 


self-radiation  impedance,  written 

Re(ZSdf)  =  j(ka)2pcA.  (18) 

Therefore,  for  small  ka,  the  mutual-radiation  impedance  can  be  expressed  as 

Zi2  =  +i““)  =  Re(2«lf)^.  (19) 

The  so-called  "Pritchard  Approximation"  refers  to  the  application  of 
Equation  (19)  to  find  the  mutual  radiation  impedance  of  an  array  of  identical 
transducers,  not  restricted  to  baffled  pistons,  i.e..  Equation  (19)  has  been 
applied  to  volumetric  arrays.  To  calculate  the  mutual-radiation  impedance, 
the  wave  number,  the  distance  between  transducers,  and  the  self-radiation 
resistance  of  the  transducer  must  be  known  in  order  to  apply  Pritchard's 
approximation. 
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C.  THE  ATILA  CODE  AND  THE  HARMONIC  ANALYSIS 
1.  The  Code 

ATILA  is  a  finite  element  code  developed  at  Institut  Superieur 
d'Electronique  du  Nord  (ISEN)  in  France,  and  is  intended  for  finite  element 
analysis  of  underwater  transducers.  It  uses  the  variational  formulation  of  the 
governing  mathematical  equations  of  fluid-structure  interactions  [Refs. 
6,7,8,9,10]. 

ATILA  can  solve  a  variety  of  problems,  including:  Static  Analysis, 
Modal  Analysis,  and  Harmonic  Analysis. 

Static  Analysis  yields  the  displacement  field,  the  pressure  field,  and 
the  electrical  potential  of  an  elastic  or  piezoelectric  structure.  Modal 
Analysis  computes  the  eigenfiequencies  and  normal  modes  of  an  elastic, 
piezoelectric  or  magnetostrictive  structure.  Harmonic  Analysis  corresponds 
to  a  forced  vibration  problem  where  the  excitation  of  the  structure  comes 
from  an  incident  wave  or  from  the  voltage  applied  across  the  electrical 
terminals  at  a  prescribed  frequency.  The  analysis  can  be  for  a  radiation  or  a 
scattering  problem,  for  the  given  frequency  in  which  the  displacement  field, 
electrical  potential,  electrical  impedance  of  the  structure,  pressure  field. 
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reflection  and  transmission  coefficients,  and  transmitting  voltage  response 
are  found. 

For  scattering  problems,  an  incident  wave  is  defined  by  the  user.  A 
default  function  is  provided  with  ATILA  which  creates  a  plane  wave 
traveling  in  the  negative  x  axis  direction.  One  can  excite  the  structure  with 
an  arbitrary  incident  wave  by  adding  a  proper  function  at  the  end  of  the  main 
FORTRAN  program.  Appendix  A  contains  the  FORTRAN  program  used  to 
generate  the  incident  spherical  harmonic  waves  on  the  structure  of  the 
transducer. 

ATILA  calculates  either  the  resulting  "total  pressure"  or  only  the 
"scattered  pressure."  The  total  pressure,  for  the  code,  is  the  sum  of  the 
pressure  generated  by  the  electrical  potential  and  the  pressure  of  the 
scattered  wave,  plus  the  pressure  of  the  impinging  wave.  In  the  case  of 
scattered  pressure,  the  pressure  of  the  incident  wave  is  not  included . 

The  ATILA  library  has  46  different  types  of  elements,  used  to  model 
elastic,  composite,  piezoelectric,  magnetostrictive,  and  magnetic  materials 
as  well  as  fluids,  solid-fluid  interfaces,  and  radiation  dampers.  Most  of  the 
elements  use  the  same  polynomial  (quadratic)  interpolation  for  both 
geometry  and  field  variation  (isoparametric  elements). 
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2.  Harmonic  Analysis  of  an  Elastic  or  Piezoelectric  Structure 

Excitation  by  an  Impinging  Wave 

ATILA  transforms  equations  of  the  Harmonic  Analysis  of  a 
Radiating  Piezoelectric  Transducer  into  a  matrix  form.  The  motion  equation 
is  used  for  piezoelectric  and  elastic  structures,  Poisson's  equation  is  used  for 
piezoelectric  material,  and  Helmholtz's  equation  is  used  for  fluids  [Ref. 
2,11]. 

On  the  interfaces  between  the  structure  and  the  fluid,  kinematics  and 
dynamic  continuity  conditions  hold  under  the  assumption  of  no  fluid 
cavitation  at  the  fluid-structure  interface.  A  radiation  condition  is  applied  to 
the  external  fluid  boundary.  Furthermore,  to  represent  the  condition  of  an 
impinging  wave,  the  pressure  and  the  flux  fields  are  separated  into  incident 
d  scattered  parts. 

In  ATILA,  the  governing  equations  are  written  in  matrix  form  as 


[Kuu]-co2[M]  [Ku<>] 
[Kuo]T  [Koo] 

-[L]T  [0]T 


-[L] 

[0] 

[H]  [M,] 

A  B 


’  u  ' 

$ 

= 

_  pes  . 

i[G]Pes  + 


F  —  [L]P  i 
— Q 

Vi 


par 


(20) 


where: 
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u 

0 

Pi 

p« 

F 

Q 


K. 

[KJ 

PW 

[M] 

[M,] 

[H] 

[L] 

[G] 

[0] 

(0 

P 

c 

A 

B 

[]T 


vector  of  the  nodal  values  of  the  components  of  the 
displacement  field 

vector  of  the  nodal  values  of  the  electrical 
potential 

vector  of  the  nodal  values  of  the  incident  pressure 
field 

vector  of  the  nodal  values  of  the  scattered  pressure 
field 

vector  of  the  nodal  values  of  the  applied  forces 

vector  of  the  nodal  values  of  the  electrical  charges 

vector  of  the  nodal  values  of  the  integrated  normal 

derivative  of  the  incident  pressure  on  the  surface 

boundary  S  (the  externally  applied  pressure  field  is 

proportional  to  the  externally  applied  flux) 

stiffness  matrix 

piezoelectric  matrix 

dielectric  matrix 

consistent  mass  matrix 

consistent  fluid  (pseudo-)  mass  matrix 

fluid  (pseudo-)  stiffness  matrix 

coupling  matrix  at  the  fluid  structure  interface 

(connectivity  matrix) 

complex  linear  operator  that  is  frequency 

dependent 

zero  matrix 

angular  frequency 

fluid  density 

fluid  sound  speed 

constant  of  the  material 

constant  of  the  material 

the  superscript  T  means  the  matrices  transpose. 


The  incident  flux  field  can  be  expressed  with  the  nc^al  values  of  the 
incident  pressure  normal  derivative, 
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Here,  [D]  is  a  matrix  that  the  code  dispenses  for  the  damping 
elements. 

The  constants  A  and  B  are 

A  =  p2c2co2  (22) 

and 

B  =  p2c2  .  (23) 

Note  that  ATILA  assumes  a  harmonic  time  dependance  of  e)0*.  The 
incident  pressure  is  provided  by  the  user  through  a  FORTRAN  program  (see 
Appendix  A),  and  for  each  input  frequency  the  code  outputs  the  complex 
displacement,  the  complex  pressure,  rotational  and  electric  potential  fields  at 
each  node,  and  the  complex  electrical  impedance  and  admittance.  For  a 
more  detailed  discusion  of  ATILA's  operation,  the  reader  is  advised  to 
consult  the  manual  [Ref.  2]. 


D.  FORCED  VIBRATION  OF  A  SPHERICAL  SHELL. 

1.  Theory 

The  mathematical  equations  governing  the  true  vibrations  and 
deformation  of  thin  elastic  shells  were  first  derived  by  Love  in  1888  [Ref. 
12].  It  assumes  the  following  postulates: 

♦  The  shell  is  thin,  compared  to  the  smallest  radius  of  curve  of  shell 

♦  Deflections  of  the  shell  are  small,  relative  to  the  shell's  thickness 

♦  There  is  no  transverse  normal  stress  acting  on  planes  parallel  to  the 
middle  surface  of  shell 

♦  There  are  no  changes  to  the  normals  of  the  reference  surface  after  and 
no  changes  in  length  during  deformation 

The  equations  are 

uee  +  uecot0  +  (l  +\))we-ucot29-\)u-^(3-'U)^^V4,  +  ^(l  +  'o)-^ve<»  + 

=  (25) 

|(1  -  v>)[vM  +  v9cot  6]  +  1  -  u)(2  -  csc20)v  +  ( 1  +v)-^w9  +  j(3-u)^u^ 

-^v#  +  i(l+u)-Lu<*  =  j7V„,  (26) 
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and 


-( 1  +  \))[ue  +  u  cot  0  +  ~qV^  +  2w]  +  *^a2aa  =  4wtt , 

cp 

where  the  subscripts  0,  <j>,  and  t  indicate  partial  differentiation  and 

v  Poisson's  ratio, 

h  shell  thickness. 

Cp  :  shell's  "plate"  velocity. 

E  Young's  modulus. 

aa  :  outward  normal  stress  applied  to  shell. 

a  spherical  shell  radius. 

0,4>,p  :  spherical  coordinates  (see  Figure  1). 
u,v,w  :  components  of  displacement  (see  Figure  4). 


Figure  4.  Components  of  Displacement. 


(27) 
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Cp 


To  calculate  the  shell's  plate  velocity,  Cp, 

_  m/g 

V  Ps(t— u*) 


(27) 


is  used  [Ref.l]. 

Here,  ps  is  the  density  of  shell's  material. 

By  a  series  of  manipulations  detailed  in  Love's  work,  u  and  v  are 
eliminated  from  the  three  equations,  resulting  in 


Wee  +  cot0We  + 


(Q2+l-v)(Q2-2(l+t))] 

Q2-(l-u2) 


J2j-1w  = 
sin20  J 


-a' 


Eh[Q2-(l 


+o,„cote+(£J2  + 1  -'0-^)0. 


(28) 


where  the  frequency  dependent  term,  Q,  is 

n  =  f,  (29) 

Cp 

and  all  (})  dependence  is  of  the  form  e,m*.  A  harmonic  time  dependence  of  the 
form  e'“‘has  been  assumed. 

Here, 

0)  =  2jtf,  (30) 


and  f  is  the  excitation  frequency. 
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If  the  new  independent  variable  Tj  is  introduced: 


T|  -  cos©. 


The  derivations  with  respect  to  0  can  be  replaced  by 


-J 


Writing  Equation  (28)  in  terms  of  (32)  and  (33): 


( 1  - q 2)W„T1  - qWr,  +  ~==r(~J  1  -T\2  Wt,) 

Ji-ir 


(fi2+i-u) 


i)  [Q2-2  (I+t>)1  _  m2  1  .  _ 
-  (1-u2)  1  — ri2  _ 


-a2(l-\)2)  ' 

Eh[«2-(l-u2)]y 


-(l-q2)aann  "^an  +j==(-Jl -q2  aaT1)  +  {(Q2  +  1  -d)--^i  aa , 
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and  it  can  be  seen  that  w  and  o#  are  being  operated  on  by  Legendre's 
differential  equation. 

Representing  the  applied  stress  ca,  and  the  normal  component,  w,  in 
terms  of  the  Legendre  polynomial  yields 

W  =  Xn,m  w„mPn,(T|)e'm^  (35) 

and 

CTa  =  Xn,mFnniPS'(Ti)e^,  (36) 

where  P™(n)  are  the  associated  Legendre  eigenfunctions  weighted  by  the 
excitation  stress  amplitude  F^,  and  the  normal  displacement  amplitude  wm. 

Legendre's  equation  can  be  written 

(1  -  2ti|jPS*(t|)  -  -^PiTai)  =  -x„p?(n).  (37) 

where 


A.„  =  n(n  +  1). 


(38) 


and  by  substitution  of  (36)  and  (37)  into  (40)  using  (39),  one  is  left  with  the 
algebraic  equation  related  Fm  to  w^: 


(Q2+I-v)  [Q2-2(1+v)J 


Q2-<1-v2) 


Wnm  — 


-df 


hCpPs[£22-  (1— v2)] 


[Q2  +  l-v-X„]F 


nm 


(39) 
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In  analogy  with  mechanical  impedance  [Ref.  13],  the  modal 


mechanical  impedance,  Zn,  can  be  defined  as: 

y  F  nm 

Z'n  “ 


(40) 


where  the  time  dependence  is  given  by  e 

It  is  then  possible  to  write  Equation  (40)  with  modal  mechanical 
impedance,  as 


hCpp, 

[  [Q2-2  (l+\»)  (02+l-\>-X„)  -X„(1+u)2 

afl 

1  Q2+I-t>-XI1 

(41) 


Now,  for  the  scattering  problem,  the  pressure  directed  radially 
inward  will  be  the  sum  of  the  incidental  spherical  wave  and  the  outgoing 
spherical  wave. 

Fnm  =  -{j  „(ka)I„m  +  hi2)(ka)Rnm } ,  (42) 

where  the  term  jn(ka)Inm  refers  the  incident  spherical  wave,  p‘,  and  the  term 
hn2)(ka)Rnm  is  the  scattered  spherical  wave,  ps.  The  constants  1^  and  R^  are 
amplitudes  of  incident  and  scattered  waves  respectively. 

pi  =  InmInmjn(kr)P!f(T1)ei^,  (43) 


ps  =  ZnmRnmh(n2)(kr)P!T(Ti)e^, 


(44) 
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w  =  lnm  WnmP^(Tl)eln^. 


(45) 


Equations  (43),  (44),  and  (45)  express  the  pressures  and  radial 
displacements  in  terms  of  spherical  harmonics. 

Applying  Equation  (42)  in  Equation  (40)  yields 

Inm  jn(ka)  +  Rnmhn  (ka)  =  — iOlZnWnm .  (46) 

Furthermore,  applying  Euler's  Equation  [Ref.  13]  at  the  surface  of 

sphere, 

!J~— pifrJ—  W 

where,  pf  is  the  density  of  the  fluid. 

The  pressure  at  surface,  p,  is  the  sum  of  the  incident  pressure,  p‘,  and 
the  scattered  pressure,  p\  as  in 

P  =  Pi+Ps-  (48) 

Now,  applying  Equations  (43),  (44),  (45)  and  (46)  on  (47), 

Inm  y  n(ka)  Rnmhn  (ka)  —  COPf  Cf  Wnm  •  (49) 

CJ\/ 

where  j„(ka)  and  hn  (ka)  are  the  first  derivatives  of  the  spherical  Bessel 

and  Hankel  functions. 

From  Equation  (49), 
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(50) 


Wnm  —  (jjpfQf  {inm  jn(ks)  +  R-nmhn  ^(ka)  j . 

Now,  Equation  (50)  can  be  utilized  in  Equation  (46),  yielding 
Pf  Cf  [Inm  jn(ka)  +  Rnmhn  (ka)]  =  —  iZn[Inm  jn(ka)  Rnmhn  ^(ka)].  (51) 

Finally,  we  have  an  equation  for  calculation  of  the  constant  term  R^, 


R 


nm  — 


iZ,Qn(ka)+PfCfja(ka) 
iZnh^CkaJ  +  PfCfhfflca) 


lnm 


(52) 


The  values  of  R^,  for  each  n  and  m  will  be  the  diagonal  elements  of 
the  T-matrix  of  the  thin  spherical  shell;  each  1^,  equals  one. 

2.  Thinness  Criteria 

The  theory  of  thin  elastic  shells  is  based  upon  the  postulate  that 
shells  are  thin.  No  exact  definition  of  thinness  is  available,  however. 

From  Junger  &  Feit  [Ref.  14]  we  have  a  definition  of  "thick" 
plates,  expressed  as 

h>B-  (53) 

where  A.s  is  the  shear  wavelength.  An  ad  hoc  treatment  to  account  for  thick 
shells  that  takes  into  account  bending  stresses  introduces  the  parameter  P, 
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which  is  a  function  of  the  thickness,  h,  and  the  middle  surface  of  the  shell. 


an,,  as  in 


(54) 


A  rule  of  thumb  suggested  by  Kraus  [Ref.  15],  is  that  the  thickness 
should  be  less  than  one  tenth  of  the  radius  of  curvature  of  the  reference 
surface,  i.e.,  p2<  1/1200. 

On  the  other  hand,  when  Kraus  examines  the  dynamic  analysis  of 
shells  [Ref.  15],  he  affirms  that  the  theory  based  upon  Love's  postulate  gives 
a  reliable  result  in  the  range 

0<O<fl„  (55) 


where  Qs  is  the  dimensionless  frequency  of  the  first  thickness  shear  mode: 


Os  = 


mSa<n 

CT  ’ 


(56) 


where 


Cl),  = 


JlCj 

~  9 


(57) 


and 


i 


(58) 
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Applying  Equation  (27)  and  (29)  in  Equation  (55)  yields 


n  ^  /p»(‘-u2)  ,|V1  ^  ±  /mEELa  o 

U<y  £  (J&&m  <  10V  E  ^s^m- 


Now,  after  simplifying  Equation  (59)  and 


which  can  be  rewritten  as 


(59) 

substituting  Equation  (57), 

(60) 


h<  0.05A,*,  (61) 

and  is  exactly  the  same  definition  offered  by  Junger  &  Feit  [Ref.  14]  listed 
in  Equation  (53). 

At  this  point,  the  thinness  criterion  for  applying  the  theory  of  thin 
elastic  shell  can  be  considered  to  be  the  expression  in  Equation  (61).  This 
criterion  is  obeyed  in  all  subsequent  calculations. 


IU.  SPHERICAL  THREE-DIMENSIONAL  MODEL 


A.  INTRODUCTION 

In  order  to  apply  the  T-matrix  method  to  modeling  a  dense  sonar  array, 
the  radiated  and  scattered  pressure  from  an  individual  transducer  must  be 
calculated.  Using  finite  element  analysis,  the  pressure  scattered  from  a 
transducer  due  to  an  incident  spherical  harmonic  wave  can  be  computed 
very  accurately.  To  employ  the  Finite  Element  Code  ATILA,  a 
three-dimensional  model  of  an  elastic  spherical  shell  is  used. 


B.  THE  MODEL 

A  0.5  m  outer  radius,  one  centimeter-thick  spherical  steel  shell  in  water 
is  modeled.  An  ATILA  model  is  realized  using  3D  elastic  elements,  3D 
fluid  elements,  2D  solid-fluid  interface  elements,  and  2D  radiation  surface 
elements.  The  thickness  was  chosen  in  an  attempt  to  comply  with  aspect 
ratio  constraints  [Ref.  2]  for  elements  used  by  ATILA,  and  to  be  within  the 
theoretical  conditions  of  thinness  (Chapter  II). 
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In  order  to  calculate  the  limit  of  thickness  to  which  we  can  apply  the 
theory  of  thin  elastic  shells  (Chapter  II),  the  following  parameters  of  the 
model  (for  steel)  are  applied  in  Equation  (58) 

E  =  0.125xl012  (Y dung's  modulus) 
ps  =  7500  Kg/m3 ,  and 
\)  =  0.33  (Poisson's  ratio). 

The  value  of  CT  is  then  3282.83  m/sec2,  which  for  a  chosen  operating 
frequency  of 474  Hz  (corresponding  to  ka=l)  results  in  an  upper  limit  for  the 
shell  thickness  of  0.346  m,  and  therefore  a  thickness  of  one  centimeter  is 
verifiably  inside  the  limit. 

The  complete  ATILA  code  for  the  elastic  shell  is  presented  in 
Appendix  B.  Note  that,  to  write  an  ATILA  code  for  a  piezoelectric  sphere, 
the  same  model  can  be  used  as  for  the  elastic  sphere.  In  this  instance,  the 
ATILA  code  can  Handle  the  piezoelectric  sphere  by  changing  the  elements 
from  elastic  to  piezoelectric  type. 

The  model  is  composed  of  2746  nodes  and  288  elements.  The  mesh  is 
shown  in  Figures  5  and  6. 
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Figure  6.  The  Model,  a*  Seen  from  Above. 


The  material  properties  of  the  shell  are  listed  on  the  first  page  of 
Appendix  B.  The  format  is  described  in  the  ATTLA  user's  manual  [Ref.  2]. 

The  elements  used  in  the  model  followed  the  pattern  from  the  manual 
[Ref.  2],  and  are  presented  in  Table  1. 


TABLE  1.  ATILA  ELEMENTS 

Region 

Element 

Geometry 

Elastic  shell 

SHEL06C 

6-node  triangular 

Interface  solid  -  fluid 

TRIA12I 

12-node  triangular 

Fluid 

PRIS15F 

15-node  triangular  prism 

Radiation  Surface 

TRIA06R 

6-node  triangular 

32 


IV.  RESULTS 


A.  THE  T-MATRIX  ELEMENTS  AND  THE  ATILA  CODE. 

The  definition  of  the  T-matrix  adopted  in  this  work  relates  outgoing 
spherical  waves  (i.e.  Hankel  functions)  to  incident  standing  waves  (i.e., 
spherical  Bessel  functions).  However,  using  ATILA,  calculation  of  the 
coefficients  of  the  scattered  pressure  for  each  harmonic,  {Ps  },  requires  that 
the  incident  wave  [Appendix  A]  be  an  incoming  traveling  wave  rather  than  a 
standing  wave. 

To  use  the  results  of  the  ATILA  code  in  calculations,  Equation  (46)  and 
its  respective  derivative.  Equation  (49)  is  adapted  to 

Inm  hn  (ka)+Rnin  hn  (ka)  =  —  icoZn  Wnm,  (62) 

and  its  first  derivative, 

Inm  hn  (ka)+ Rnm  hn2)/(ka)  =  (OpfCf  W  nm  •  (63) 

The  coefficients  Rnm  are  _the  amplitudes  of  the  scattered  wave 
components  and  the  coefficients  Wnm  are  the  amplitudes  of  the  components 
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of  the  radial  displacement  coefficients  for  an  incoming  wave  with 
coefficients  Inm  • 


Equation  (52)  (Chapter  II)  is  used  to  calculate  the  coefficients  from 
an  incident  standing  wave,  from  where  the  T-Matrix  elements  are  obtained. 
For  incident  incoming  traveling  waves,  the  elements  can  be  calculated 
as  follows. 

From  Equation  (62)  and  (63), 


Rnm= 


uUV 


iZnhn  (ka}+  pfCf  hn  (ka) 

I  i/.„hi2)/(ka)  +  pfCfhi2)(ka)  J 


run  > 


We  can  relate  Rnm  and  Rnm  .  We,  first  .define 


Xnm  — 


iZn(i  y„(ka))  +  pfCf  (i  yn(ka)) 

L  iZ„h!,2)/(ka)  +  PfCf  hn2)  (ka)  J 


flnm> 


(64) 


(65) 


where  yn(ka)  is  the  spherical  Bessel  function  of  the  second  kind. 

Applying  Equation  (65)  on  Equations  (52)  and  (64),  respectively  : 

Rnm-X„m=-1  (66) 


and 


Rnm  — Xnrn  —  R 


‘•nm 


Solving  Equations  (66)  and  (67)  in  terms  of  Rnm: 


Rnm  —  J  {Rnm 


(67) 


(68) 
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Equation  (68)  is  used  to  calculate  the  coefficients,  and  hence  the 
T-Matrix  elements  for  incident  standing  waves,  from  the  Rnm  coefficients 

calculated  when  ATILA  code  uses  incoming  incident  waves  to  compute  the 
scattered  pressure. 

The  goal  of  this  research  is  to  compare  the  T-Matrix  elements 

computed  by  scattered  field  pressure  from  the  ATILA  code  with  those 

determined  by  theoretical  means  (Appendix  C). 

The  T-Matrix  for  an  elastic  spherical  shell  is  a  diagonal  matrix  where 
the  non-zero  elements  equal  the  R^  coefficients,  for  Inm  =Inm=  1 .  In  this 

research,  calculations  of  Rnm  were  performed  up  to  the  harmonic  of  second 
order,  and  the  T-Matrix  has  the  following  format: 


Roo 

0 

0 

0 

0 

0 

0 

0 

0  ' 

0 

Ri-1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Rio 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Rn 

0 

0 

0 

0 

0 

0 

0 

0 

0 

R2-2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

R2-1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

R20 

0 

0 

0 

0 

0 

0 

0 

0 

0 

R21 

0 

0 

0 

0 

0 

0 

0 

0 

0 

R22  , 

In  the  T-matrix,  for  instance,  the  element  R^,  corresponds  to  the 
spherical  harmonic  with  n  equal  to  two  and  m  equal  to  negative  one.  The 
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T-Matrix  elements  were  theoretically  computed  using  thin  shell  theory 

» 

(Appendix  C)  and  the  results  are  presented  in  Table  II. 


TABLE  II.  THE  ANALYTIC  DIAGONAL  ENTRIES  OF  THE  T-MATRIX 

ELEMENT 

REAL 

IMAG 

MAGNITUDE 

PHASE 

PART 

PART 

(DEGREES) 

T„  =R<x> 

-1.3465E-02 

1.1525E-01 

1.1604E-01 

9.6663E+01 

“R1-1 

-6.0255E-03 

7.7390E-02 

7.7624E-02 

9.4452E+01 

Tj3  =R10 

-6.0255E-03 

7.7390E-02 

7.7624E-02 

9.4452E+01 

^4*  =R11 

-6.0255E-03 

7.7390E-02 

7.7624E-02 

9.4452E+01 

^*55  =R2-2 

-6.1773E-04 

-2.4847E-02 

2.4854E-02 

-9.1424E+01 

T„  ~R2-1 

-6.1773E-04 

-2.4847E-02 

2.4854E-02 

-9.1424E+01 

^77  =R20 

-6.1773E-04 

-2.4847E-02 

2.4854E-02 

-9.1424E+01 

=R21 

-6.1773E-04 

-2.4847E-02 

2.4854E-02 

-9.1424E+01 

"^99  =R22 

-6.1773E-04 

-2.4847E-02 

2.4854E-02 

-9.1424E+01 

Incident  waves  of  the  form 


Pinc(r,  0,  <}>)  =  hiVracos  0)ein*eivrt 

[Appendix  A]  for  particular  values  of  n  and  m  were  applied  to  the 
finite-element  model  spherical  shell  using  the  ATILA  code  and  the  resulting 
scattered  pressures  were  calculated  .  For  each  incident  wave  of  a  particular 
spherical  harmonic  component  (one  n  and  m),  the  resulting  reflection 
coefficients  for  all  the  spherical  harmonic  components  (all  n  and  m)  were 
calculated  from  nodal  values  of  the  scattered  pressure  on  the  shell  surface  by 
a  least-squares  fitting  procedure  employing  singular  value  decomposition 
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[Appendix  C].  The  results  for  the  diagonal  elements  of  the  T-matrix  are 
given  in  Table  III.  The  last  two  columns  of  Table  in  give  the  error  in  the 
ATILA  results  for  the  diagonal  components  compared  to  the  theoretical 
results  of  Table  II. 


TABLE  Hi.  THE  ATILA  DIAGONAL  T-MATRIX  ELEMENTS 

MAG. 

PHASE 

ELEMENT 

REAL 

IMAG 

MAGNITUDE 

PHASE 

ERROR 

ERROR 

PART 

PART 

(DEGREES) 

(PERCENT) 

(DEGREES) 

T„ 

-1.2624E-02 

1.1429E-01 

1.1498E-01 

9.6303E+01 

-0.91 

-0.36 

T* 

-6.8902E-03 

8.2376E-02 

8.2664E-02 

9.4781  E+01 

6.49 

0.33 

T„ 

-5.7422E-03 

8.3113E-02 

8.331  IE-02 

9.3952E+01 

7.33 

-0.5 

-6.8868E-03 

8.2363E-02 

8.2649E-02 

9.4766E+01 

6.47 

0.31 

T„ 

-2.8487E-03 

-1.5203E-02 

1.S468E-02 

-1.0061E+02 

-37.8 

-9.19 

T„ 

-2.7899E-03 

-1.5159E-02 

1.5413E-02 

-1.0043E+02 

-38 

-9.01 

T„ 

-2.4889E-03 

-1.5287E-02 

1.5489E-02 

-9.9247E+01 

-37.7 

-7.82 

T„ 

-2.7847E-03 

-1.S164E-02 

1.5418E-02 

-1.0041E+02 

-38 

-8.99 

T* 

-2.8227E-03 

-1.5182E-02 

1.5443E-02 

-1.0053E+02 

-37.9 

-9.11 

The  T-matrix  for  a  spherically  symmetric  scattered  should  be  diagonal, 
that  is,  there  should  be  only  one  spherical  harmonic  component  of  the 
scattered  wave  for  a  single-component  incident  wave,  and  it  should  be  the 
same  component  as  the  incident  wave.  Nonzero  off-diagonal  components 
computed  for  the  T-matrix  of  the  spherical  elastic  shell  can  be  termed 
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'‘leakage"  (in  some  abstract  sense).  The  existence  of  leakage  indicates  a 
weakness  in  the  model. 


Nonzero  off-diagonal  components  of  the  T-matrix  which  were  greater 
than  10"14  of  the  diagonal  element  were  found  for  a  few  of  the  ATTLA 
results.  These  are  given  in  Table  IV  below. 


TABLE  IV. 

SIGNIFICANT  OFF-DIAGONAL  T-MATRIX  ELEMENTS 

element 

real 

imag 

mag 

phase 

mag  T, 
rel  to  T, 

phase  T, 
rel  to  T„ 

t71 

-7.6679E-04 

5.3764E-04 

9.3650E-04 

1 .4496E+02 

8.14E-03 

48.7 

2.2172E-05 

7.2191E-06 

2.331 8E-05 

1.8035E+01 

2.82E-04 

-76.7 

-8.9863E-06 

-3.8094E-07 

8.9944E-06 

-1.7757E+02 

1.09E-04 

-272.3 

T* 

1.6194E-05 

1.4535E-05 

2.1761  E-05 

4.191 1E+01 

1.41  E-03 

142.5 

T„ 

-1.1214E-05 

4.1657E-08 

1.1214E-05 

1.7979E+02 

7.28E-04 

280.2 

T„ 

5.5328E-03 

-4.3539E-03 

7.0405E-03 

-3.8200E+01 

4.55E-01 

61 

3.5575E-05 

3.6843E-05 

5.1215E-05 

4.6003E+01 

3.32E-03 

146.4 

T. 

9.8906E-03 

7.9364E-03 

1.2681E-02 

3.8744E+01 

8.21  E-01 

139.3 

It  is  seen  that  there  is  significant  leakage  from  the  2,0  to  the  0,0 
component  T17,  and  from  the  2,2  to  the  2,-2  component  T59,  for  the 
finite-element  model  employed.  This  could  indicate  that  either  the  mesh  is 
too  coarse,  particularly  that  of  the  shell,  or  that  the  radius  of  the  radiation 
damping  elements  which  terminate  the  fluid  mesh  is  too  small.  These 
possibilities  remain  to  be  investigated. 
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B.  THE  ACOUSTIC  IMPEDANCE 


To  further  validate  the  three-dimensional  model  used  in  the  ATILA 
code,  the  acoustic  impedance  computed  analytically  and  with  ATILA  output 
were  also  compared.  The  theoretical  acoustic  impedance,  for  each 
harmonic,  is  given  by  Equation  (42). 

On  the  other  hand,  output  from  ATILA  can  be  used  to  determine  the 
impedance  from  Equation  (41),  where  will  be  the  total  field  pressure. 
Total  pressures  are  output  from  ATILA  when  the  pressure  command  in 
ATILA  code  (Appendix  B)  is  "PRESSURE  TOTAL"  [Reference  2].  The 
terms  are  also  computed  from  the  ATILA  output.  Averaged  values  of 
Fnm  and  Wnrn  were  used  to  compute  an  approximation  to  the  acoustic 
impedance  of  each  harmonic.  The  results  are  presented  in  Table  VI. 


TABLE  V.  RESULTS  OF  ACOUSTIC  IMPEDANCE 

Theoretical  Values 

From  ATILA 

Error  and  Number  of  Points 

4=0.00-1  0.8396E+07 
Z,=0.00+i  0.6945E+06 

Z^.SBtME+OS-i  0.82305E+07 
z^^eessE+CM+i  0.60223E+06 

2%  and  146 

13%  and  138 
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V.  CONCLUSIONS 


The  transition  matrix,  or  “T-matrix”,  relating  the  scattered  and  incident 
acoustic  waves  for  a  thin  elastic  spherical  shell  in  a  free-field  environment 
has  been  evaluated  using  the  ATILA  finite-element  code.  A 
three-dimensional  finite-element  model  of  a  0.5-m  outer  radius,  1-cm  thick 
spherical  steel  shell  surrounded  by  water  was  developed  [Appendix  A].  In 
this  model  the  spherical  shell  surface  is  divided  into  seventy-two 
approximately  equal  area  triangular  regions.  The  AULA  code  was  used  to 
calculate  the  scattered  pressure  over  the  surface  of  the  shell  for  incident 
waves  represented  as  products  of  radial  Hankel  functions  and  spherical 
harmonics.  This  was  done  for  all  components  through  the  second  order  in 
the  Hankel  function.  The  chosen  driving  frequency  was  474  Hz, 
corresponding  to  a  value  of  ka=l,  where  k  is  the  wavenumber  of  sound  in 
water  and  a  is  the  radius  of  the  shell.  The  ATILA  results  were  compared 
with  the  results  of  analytical  thin  shell  theory.  For  the  diagonal  elements  of 
the  T-matrix  ,  the  two  results  were  found  to  agree  within  one,  seven,  and 
thirty-eight  percent,  respectively,  for  the  zeroth,  first,  and  second  order 
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components.  The  ATILA  results  also  produce  a  few  nonzero  off-diagonal 
T-matrix  elements  which  do  not  appear  in  the  theoretical  results,  and  which 
should  not  be  present  according  to  symmetry  considerations.  The  reason  for 
the  appearance  of  nonzero  off-diagonal  elements  in  the  ATILA  results  is  not 
known  at  this  time  and  warrants  further  investigation.  Also  computed  for 
each  harmonic  component  was  the  modal  acoustical  impedance  of  the  shell. 
These  results  agreed  with  the  results  of  thin  shell  theory  within  two  percent 
for  the  zeroth  order  component  and  thirteen  percent  for  the  first  order 
components. 
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APPENDIX  A 


FUNCTION  INCPRE(X,Y,Z,K) 


PROGRAM  BY  ARTHUR  LOBO  DA  COSTA  RUIZ  12/23/93 


FUNCTION: 

COMPUTES  THE  INCIDENT  SPHERICAL  HARMONIC  PRESSURE  AT  THE  POINT 
(X,Y,Z),FOR  THE  WAVENUMBER  K 

VARIABLES  INPUT: 

X,Y,Z:  CARTESIAN  COORDINATES  OF  THE  POINT 
K:  WAVENUMBER 

VARIABLES  OUTPUT: 

RADIUS  :R  IN  METERS 

PHI  IN  DEGREE  (HORIZONTAL  ANGLE) 

THETA  IN  RADIANS  (AZIMUTAL  ANGLE) 

INCPRE  :  PRESSURE  AT  POINT  IN  PASCAL 


DOUBLE  PRECISION  K,X,Y,Z,R, PHI, THETA, KR 
REAL*8  PMN(-2:2,0:2),LOUT 
INTEGER  N.M.NMAX 
COMPLEX*16  INCPRE, H(0:2),HOUT 

* 

*  N  AND  M  ARE  ORDERS  FOR  HANKEL  AND  LEGENDRE  FUNCTIONS 


N=0 

M=0 

*  HOUT  AND  LOUT  ARE  HANKEL  AND  LEGENDRE  OUTPUTS 

*  REF  TON  AND  M 

NMAX=2 

*  TRANSFORM  CARTESIAN  COORDINATES  (X,Y,Z) 

*  INTO  SPHERICAL  COORDINATES 

*  R(RADIUS),  PHI(AZIMUTAL  ANGLE)  AND  THETA(POLAR  ANGLE) 

R=DSQRT  (X*X+Y*  Y+Z*Z) 

PHI=DATAN2D(Y,X) 

IF((X.EQ.O).AND.(Y.EQ.O))  PHI=O.ODO 
THETA=DACOS(Z/R) 

KR=K*R 

* 

*  NMAX  IS  THE  MAXIMUM  NUMBER  OF  HARMONICS 
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o  o 


CALL  HANKEL(KR ,  NMAX,  H) 


SUBROUTINE  HANKEL  RETURNS  SPHERICAL  HANKEL 
FUNCTIONS  JN  AND  I  YN 


CALL  LEGNDR(TH ETA, NMAX, PMN) 
LOUT=PMN(M,N) 


*  SUBROUTINE  LEGNDR  RETURNS  ASSOCIATE  LEGENDRE  FUNTION 

* 


INCPRE=H(N)*LOUT*DCMPLX(DCOSD(M*PHI),DSIND(M"PHI)) 

IF  ((R.LE.0.501  ).AND.(R.GT.0.499))  THEN 

PRINT  *,'X,Y,Z  =",X,Y,Z 

PRINT  *,"PHI,THETA  PHI, THETA 

PRINT  YK,R,KR  =",K,R,KR 


«***************a 


PRINT  ",'REAL  HANKEL  *",HOUT 
PRINT  ‘."LEGENDRE  =’,LOUT 
PRINT  MNCPRE 
PRINT  ", 

ELSE 

CONTINUE 
ENDIF 
RETURN 
END 


SUBROUTINE  HANKEL(X.NMAX.H) 

IMPLICIT  REAL*8  (A-H.O-Z) 

COMPLEX*16  H(0:NMAX) 

C 

C  GIVEN  THE  VARIABLE  X.  AND  THE  MAXIMUM  ORDER  NMAX, 

C  THIS  ROUTINE  GENERATES  THE  SPHERICAL  HANKEL  FUNCTION  HN 
C  FOR  ALL  N  FROM  0  TO  NMAX  (INCLUSIVE) 

C  INPUT: 

C  X  =  DOUBLE  PREC.  VARIABLE  (RADIUS) 

C  NMAX  =  INTEGER  MAXIMUM  ORDER  OF  BESSEL  FUNCTIONS  DESIRED 
C  OUTPUT: 

C  H(N)  =  ARRAY  OF  SPHERICAL  HANKEL  FUNCTIONS  HN(X),  WHERE 
C  HN  =  JN  +  I  YN  C  C  THIS  ROUTINE  IS  BASED  ON  THE  RECURSION  FORMULAE 
C  FROM  ABRAMOWITZ  &  STEGUN:  10.1.10  &  10.1.15,  PP.438-9 
C  THE  F  S  ARE  THE  COEFFICIENTS  OF  ORDER  N  &  -(N+1 ), 

C  THE  FO1  S  ARE  OLD  F  S,  FOR  RECURSION 
C 

IF  (  X  .LE.  0.0D0 )  THEN 
H(0)  =  DCMPLX(1 .000,-1 .0035) 

DO  2  N  =  1 ,  NMAX 
H(N)  =  CMPLX(0.0D0,-1.0D35) 

2  CONTINUE 
RETURN 
END  IF 
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o  o 


SX  «  OSIN(X) 

CX  «  DCOS(X) 

XINV  *  1 ,000/X 
MIN  *-1.000 
FN  *  XINV 
FMN  *  0.000 
FNO  =  FMN 
FMNO  *  FN 
DO  4  N  *  0,  NMAX 

H(N)  =  CMPLX(  FN*SX  +  M1N*FMN*CX,  -FN*CX  +  M1N*FMN*SX  ) 

T1  « (2*N+1)*XINV 

T2  =>  T1  *FN  -  FNO 

FNO  *  FN 

FN*T2 

T2  *  -T1*FMN  -  FMNO 
FMNO  *  FMN 
FMN  =  T2 
MIN* -MIN 
4  CONTINUE 
RETURN 
END 
C 


SUBROUTINE  LEGNDR(THETA,NMAX,PMN) 

IMPLICIT  REAL*8  (A-H.O-Z) 

REAL* 8  PMN(-NMAX:NMAX,0:NMAX) 

C 

C  GIVEN  THE  VARIABLE  THETA,  AND  THE  MAXIMUM  ORDER  NMAX, 

C  THIS  ROUTINE  GENERATES  THE  ASSOC.  LEGENDRE  FUNCTIONS  PMN 
C  OF  THE  ARGUMENT  COS(THETA)  (THETA  MUST  BE  BETWEEN  0  &  PI) 

C  FOR  ALL  N  FROM  0  TO  NMAX  (INCLUSIVE) 

C  AND  FOR  ALL  M  FROM  -N  TO  N  (SOME  OTHERS  SET  TO  ZERO) 

C  INPUT: 

C  THETA  =  VARIABLE  (POLAR  ANGLE),  MUST  BE  BETWEEN  0  &  PI  (INCL.) 

C  NMAX  =  INTEGER  MAXIMUM  ORDER  OF  LEGENDRE  FUNCTIONS  DESIRED 
C  OUTPUT: 

C  PMN  *  DOUBLE  PREC.  ARRAY,  CONTAINS  ASSOC.  LEGENDRE  FNS 
C 

C  THIS  ROUTINE  IS  BASED  ON  THE  RECURSION  FORMULAE 
C  FROM  ABRAMOWITZ  &  STEGUN 
C 

X  =  DCOS(THETA) 

SINTHT  *  DSINfTHETA) 

IF  (  SINTHT  .GT.  0. )  THEN 
SININV  =  1.0D0/SINTHT 
ELSE 

SININV  =  0.0D0 
END  IF 

C  SET  VALUES  FOR  N  =  0, 1  (NMAX  MUST  BE  AT  LEAST  1) 

PMN(0,0)  =  1.0D0 
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PMN(1,0)  =  O.ODO 
PMN(-1 ,0)  ■  0.000 
PMN(0,1) »  X 
PMN(1,1)  =  -SINTHT 
PMN(-1,1)  -  SINTHT*0.5D0 

C  IN  LOOP,  TNP1  ■  2*N+1,  TNP2FC  =  (2*N+2)!,  MIN  =  (-1p(N+1) 

TNP1  =  1.0D0 
TNP2FC  >  2.000 
MIN  =>  -I.ODO 
DO  4  N  =  1,  NMAX-1 
TNP1  =  TNP1  +  2.000 
TNP2FC  »  TNP2FC  *  TNP1  *  (TNP1+1) 

MIN  =  -MIN 
DO  3  M  *  -N.  N 

PMN(M,N+i)  =  (TNP1  *X*PMN(M,N)  -  (N*M)*PMN(M,N-1))/(N-M+1) 
3  CONTINUE 

PMN(N+1  ,N)  *  0.000 
PMN(-N-I.N)  =  O.ODO 

PMN(N+1,N+1)  =  (X*PMN(N,N+1)  -  TNP1*PMN(N,N))  *  SININV 
PMN(-N-1,N+1)  =  M1N*PMN(N+1,N+1)/TNP2FC 
4  CONTINUE 
C  DO  120  N=0,NMAX 
C  DO  120  M=-N,N 
Cl 20  WRITE(6,130)  N,M,PMN(M,N) 

Cl  30  FORMAT(1X,’  N=  ',14, IX/  M='  tl4,1X,’  PMN=',F13.6) 

RETURN 

END 
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APPENDIX  B 


ft**************************************************** 

*  ATILA  MODEL  ESHELL4:  3-D  ELASTIC  SPHERICAL  SHELL,  1  CM  THICK 

*  AUTHOR:  BAKER  SHELL  RADIUS  0.5  M,  MADE  OF  NEW  CURVED  SHELL  ELEMENTS 

*  DATE:  MAR  1994  FLUID  RADIUS  2.5  M 

*  EACH  OCTANT  DIVIDED  INTO  9  REGIONS 

*  2  FLUID  LAYERS,  EACH  0.25  M  THICK, 

*  3  FLUID  LAYERS,  EACH  0.5  M  THICK 
****************************************************** 

PRINTING  =  2 
RADIATION  DIPOLAR 
ANALYSIS  HARMONIC 
MATERIAL 

25CD4SH  *  PROPERTIES  FOR  CURVED  COMPOSITE  SHELL  ELEMENTS 
0.215E+12  0.330E+00  0.750E+04  0.000E+00  O.OOOE+OO  O.OOOE+OO  & 

0.215E+12  0.330E+00  0.750E+04  O.OOOE+OO  1.0  0.000E+00 

LCPDDC  *  6 
NLOAD  =  1 
FREQUENCY  474. 

GEOMETRY 

1 

2.5  *  RADIUS  OF  RADIATION  BOUNDARY  =  2.5  M 

2 

1  0  0.01  *  THICKNESS  OF  SHELL  =  1  CM 

PRESSURE  SCATTERED 
SKYLINF  REAL 
PRECISION  DOUBLE 
NEWAXES  SPHERICAL 


0  0  0  0 

0  0 

NODES 

*  0001  *  0.490 

000 

000.0 

*  0002  *  0.490 

015 

000.0 

*  0003*  0.490 

015 

090.0 

*  0004  *  0.490 

015 

180.0 

*  0005  *  0.490 

015 

270.0 

*  0006  *  0.490 

030 

000.0 

*  0007  *  0.490 

030 

045.0 

*  0008  *  0.490 

030 

090.0 

*  0009  *  0.490 

030 

135.0 

*  0010  *  0.490 

030 

180.0 

*  0011  *  0.490 

030 

225.0 

*  0012  *  0.490 

030 

270.0 

*  0013*  0.490 

030 

315.0 

*  0014  *  0.490 

045 

000.0 

*  0015*  0.490 

045 

030.0 
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0016  *  0.490 

045 

060.0 

0017  *  0.490 

045 

090.0 

0018  *  0.490 

045 

120.0 

0019  *  0.490 

045 

150.0 

0020  *  0.490 

045 

180.0 

0021  *  0.490 

045 

210.0 

0022  *  0.490 

045 

240.0 

0023*  0.490 

045 

270.0 

0024*  0.490 

045 

300.0 

0025  *  0.490 

045 

330.0 

0026  *  0.490 

060 

000.0 

0027  *  0.490 

060 

022.5 

0028  *  0.490 

060 

045.0 

0029  *  0.490 

060 

067.5 

0030  *  0.490 

060 

090.0 

0031  *  0.490 

060 

112.5 

0032  *  0.490 

060 

135.0 

0033  *  0.490 

060 

157.5 

0034  *  0.490 

060 

180.0 

0035  *  0.490 

060 

202.5 

0036  *  0.490 

060 

225.0 

0037  *  0.490 

060 

247.5 

0038  *  0.490 

060 

270.0 

0039  *  0.490 

060 

292.5 

0040  *  0.490 

060 

315.0 

0041  *  0.490 

060 

337.5 

0042*  0.490 

075 

000.0 

0043  *  0.490 

075 

018.0 

0044  *  0.490 

075 

036.0 

0045  *  0.490 

075 

054.0 

0046  *  0.490 

075 

072.0 

0047  *  0.490 

075 

090.0 

0043  *  0.490 

075 

108.0 

0049  *  0.490 

075 

126.0 

0050  *  0.490 

075 

144.0 

0051  *  0.490 

075 

162.0 

0052  *  0.490 

075 

180.0 

0053  *  0.490 

075 

198.0 

0054  *  0.490 

075 

216.0 

0055  *  0.490 

075 

234.0 

0056  *  0.490 

075 

252.0 

0057  *  0.490 

075 

270.0 

0058  *  0.490 

075 

288.0 

0059*  0.490 

075 

306.0 

0060*  0.490 

075 

324.0 

0061  *  0.490 

075 

342.0 

0062*  0.490 

090 

000.0 

0063  *  0.490 

090 

015.0 

0064*  0.490 

090 

030.0 

0065  *  0.490 

090 

045.0 

« 

0066  *  0.490 

090 

060.0 
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*  0067  *  0.490  090  075.0 
’  0068  *  0.490  090  090.0 

*  0069  *  0.490  090  105.0 

*  0070  *  0.490  090  120.0 

*  0071  *  0.490  090  135.0 

*  0072  *  0.490  090  150.0 

*  0073  *  0.490  090  165.0 

*  0074  *  0.490  090  180.0 

*  0075  *  0.490  090  195.0 

*  0076  *  0.490  090  210.0 

*  0077  *  0.490  090  225.0 

*  0078  *  0.490  090  240.0 

*  0079  *  0.490  090  255.0 

*  0080  *  0.490  090  270.0 

*  0081  *  0.490  090  285.0 

*  0082  *  0.490  090  300.0 

*  0083  *  0.490  090  315.0 

*  0084  *  0.490  090  330.0 

*  0085  *  0.490  090  345.0 

*  0086  *  0.490  105  000.0 

*  0087  *  0.490  105  018.0 

*  0088  *  0.490  105  036.0 

*  0089  *  0.490  105  054.0 

*  0090  *  0.490  105  072.0 

*  0091  *  0.490  105  090.0 

*  0092  *  0.490  105  108.0 

*  0093  *  0.490  105  126.0 

*  0094  *  0.490  105  144.0 

*  0095  *  0.490  105  162.0 

*  0096  *  0.490  105  180.0 

*  0097  *  0.490  105  198.0 

*  0098  *  0.490  105  216.0 

*  0099  *  0.490  105  234.0 

*  0100  *  0.490  105  252.0 

*  0101  *  0.490  105  270.0 

*  0102  *  0.490  105  288.0 

*  0103  *  0.490  105  306.0 

*  0104  *  0.490  105  324.0 

*  0105  *  0.490  105  342.0 

*  0106  *  0.490  120  000.0 

*  0107  *  0.490  120  022.5 

*  0108  *  0.490  120  045.0 

*  0109  *  0.490  120  067.5 

*  0110  *  0.490  120  090.0 

*  0111  *  0.490  120  112.5 

*  0112  *  0.490  120  135.0 

*  0113  *  0.490  120  157.5 

*  0114*  0.490  120  180.0 

*  0115  *  0.490  120  202.5 

*  0116  *  0.490  120  225.0 

*  0117  *  0.490  120  247.5 


*  0118  *  0.490  120  270.0 

*  0119  *  0.490  120  292.5 

*  0120  *  0.490  120  315.0 

*  0121  *  0.490  120  337.5 

*  0122  *  0.490  135  000.0 

*  0123  *  0.490  135  030.0 

*  0124  *  0.490  135  060.0 

*  0125  *  0.490  135  090.0 

*  0126  *  0.490  135  120.0 

*  0127’  0.490  135  150.0 

*  0128  *  0.490  135  180.0 

*  0129  *  0.490  135  210.0 

*  0130  *  0.490  135  240.0 

*  0131  *  0.490  135  270.0 

*  0132  *  0.490  135  300.0 

*  0133  *  0.490  135  330.0 

*  0134  *  0.490  150  000.0 

*  0135’  0.490  150  045.0 

*  0136  *  0.490  150  090.0 

*  0137*  0.490  150  135.0 

*  0138  *  0.490  150  180.0 

*  0139  *  0.490  150  225.0 

*  0140  *  0.490  150  270.0 

*  0141  *  0.490  150  315.0 

*  0142  *  0.490  165  000.0 

*  0143  *  0.490  165  090.0 

*  0144*  0.490  165  180.0 

*  0145  *  0.490  165  270.0 

*  0146  *  0.490  180  000.0 

*  0147  *  0.000  000  000.0 

*  0148  *  0.000  000  000.0 

*  0149  *  0.000  000  000.0 

*  0150  *  0.000  000  000.0 

*  0151  *  0.000  000  000.0 

*  0152  *  0.000  000  000.0 

*  0153  *  0.000  000  000.0 

*  0154  *  0.000  000  000.0 

*  0155  *  0.000  000  000.0 

*  0156  *  0.000  000  000.0 

*  0157  *  0.000  000  000.0 

*  0158  *  0.000  000  000.0 

*  0159  *  0.000  000  000.0 

*  0160  *  0.000  000  000.0 
*  0161  *  0.000  000  000.0 
*  0162  *  0.000  000  000.0 

*  0163  *  0.000  000  000.0 

*  0164  *  0.000  000  000.0 

*  0165  *  0.000  000  000.0 

*  0166  *  0.000  000  000.0 

*  0167  *  0.000  000  000.0 

*  0168  *  0.000  000  000.0 
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0189* 

0.000 

000 

000.0 

0170* 

0.000 

000 

000.0 

0171  * 

0.000 

000 

000.0 

0172* 

0.000 

000 

000.0 

0173* 

0.000 

000 

000.0 

0174* 

0.000 

000 

000.0 

0175* 

0.000 

000 

000.0 

0176* 

0.000 

000 

000.0 

0177* 

0.000 

000 

000.0 

0178* 

0.000 

000 

000.0 

0179* 

0.000 

000 

000.0 

0180* 

0.000 

000 

000.0 

0181  * 

0.000 

000 

000.0 

0182* 

0.000 

000 

000.0 

0183* 

0.000 

000 

000.0 

0184* 

0.000 

000 

000.0 

0185* 

0.000 

000 

000.0 

0186* 

0.000 

000 

000.0 

0187* 

0.000 

000 

000.0 

0188* 

0.000 

000 

000.0 

0189* 

0.000 

000 

000.0 

0190* 

0.000 

000 

000.0 

0191  * 

0.000 

000 

000.0 

0192* 

0.000 

000 

000.0 

0193* 

0.000 

000 

000.0 

0194* 

0.000 

000 

000.0 

0195* 

0.000 

000 

000.0 

0196* 

0.000 

000 

000.0 

0197* 

0.000 

000 

000.0 

0198* 

0.000 

000 

000.0 

0199* 

0.000 

000 

000.0 

0200* 

0.000 

000 

000.0 

0201  * 

0.495 

000 

000.0 

0202* 

0.495 

015 

000.0 

0203* 

0.495 

015 

090.0 

0204* 

0.495 

015 

180.0 

0205* 

0.495 

015 

270.0 

0206* 

0.495 

030 

000.0 

0207* 

0.495 

030 

045.0 

0208* 

0.495 

030 

090.0 

0209* 

0.495 

030 

135.0 

0210* 

0.495 

030 

180.0 

0211  * 

0.495 

030 

225.0 

0212* 

0.495 

030 

270.0 

0213* 

0.495 

030 

315.0 

0214* 

0.495 

045 

000.0 

0215* 

0.495 

045 

030.0 

0216* 

0.495 

045 

060.0 

0217* 

0.495 

045 

090.0 

0218* 

0.495 

045 

120.0 

0219* 

0.495 

045 

150.0 
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0220  *  0.495 

045 

180.0 

0221  *  0.495 

045 

210.0 

0222  *  0.495 

045 

240.0 

0223  *  0.495 

045 

270.0 

0224  *  0.495 

045 

300.0 

0225  *  0.495 

045 

330.0 

0226  *  0.495 

060 

000.0 

0227  *  0.495 

060 

022.5 

0228  *  0.495 

060 

045.0 

0229  *  0.495 

060 

067.5 

0230  *  0.495 

060 

090.0 

0231  *  0.495 

060 

112.5 

0232  *  0.495 

060 

135.0 

0233  *  0.495 

060 

157.5 

0234  *  0.495 

060 

180.0 

0235*  0.495 

060 

202.5 

0236  *  0.495 

060 

225.0 

0237  *  0.495 

060 

247.5 

0238  *  0.495 

060 

270.0 

0239  *  0.495 

060 

292.5 

0240  *  0.495 

060 

315.0 

0241  *  0.495 

060 

337.5 

0242  *  0.495 

075 

000.0 

0243  *  0.495 

075 

018.0 

0244  *  0.495 

075 

036.0 

0245  *  0.495 

075 

054.0 

0246*  0.495 

075 

072.0 

0247  *  0.495 

075 

090.0 

0248  *  0.495 

075 

108.0 

0249  *  0.495 

075 

126.0 

0250*  0.495 

075 

144.0 

0251  *  0.495 

075 

162.0 

0252  *  0.495 

075 

180.0 

0253  *  0.495 

075 

198.0 

0254  *  0.495 

075 

216.0 

0255  *  0.495 

075 

234.0 

0256  *  0.495 

075 

252.0 

0257  *  0.495 

075 

270.0 

0258  *  0.495 

075 

288.0 

0259  *  0.495 

075 

306.0 

0260  *  0.495 

075 

324.0 

0261  *  0.495 

075 

342.0 

0262  *  0.495 

090 

000.0 

0263  *  0.495 

090 

015.0 

0264  *  0.495 

090 

030.0 

0265  *  0.495 

090 

045.7 

0266  *  0.495 

090 

060.0 

0267  *  0.495 

090 

075.0 

0268*  0.495 

090 

090.0 

0269  *  0.495 

090 

105.0 

0270  *  0.495 

090 

120.0 
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0271*  0.495 

090 

135.0 

0272  *  0.495 

090 

150.0 

0273  *  0.495 

090 

185.0 

0274*  0.495 

090 

180.0 

0275*  0.495 

090 

195.0 

0276*  0.495 

090 

210.0 

0277  *  0.495 

090 

225.0 

0278  *  0.495 

090 

240.0 

0279  *  0.495 

090 

255.0 

0280*  0.495 

090 

270.0 

0281  *  0.495 

090 

285.0 

0282  *  0.495 

090 

300.0 

0283  *  0.495 

090 

315.0 

0284  *  0.495 

090 

330.0 

0285  *  0.495 

090 

345.0 

0286  *  0.495 

105 

000.0 

0287  *  0.495 

105 

018.0 

0288  *  0.495 

105 

036.0 

0289  *  0.495 

105 

054.0 

0290  *  0.495 

105 

072.0 

0291  *  0.495 

105 

090.0 

0292  *  0.495 

105 

108.0 

0293  *  0.495 

105 

126.0 

0294  *  0.495 

105 

144.0 

0295  *  0.495 

105 

162.0 

0296  *  0.495 

105 

180.0 

0297*  0.495 

105 

198.0 

0298  *  0.495 

105 

216.0 

0299  *  0.495 

105 

234.0 

0300  *  0.495 

105 

252.0 

0301  *  0.495 

105 

270.0 

0302  *  0.495 

105 

288.0 

0303  *  0.495 

105 

306.0 

0304  *  0.495 

105 

324.0 

0305  *  0.495 

105 

342.0 

0306  *  0.495 

120 

000.0 

0307  *  0.495 

120 

022.5 

0308  *  0.495 

120 

045.0 

0309  *  0;495 

120 

067.5 

0310  *  0.495 

120 

090.0 

0311  *  0.495 

120 

112.5 

0312  *  0.495 

120 

135.0 

0313*  0.495 

120 

157.5 

0314  *  0.495 

120 

180.0 

0315*  0.495 

120 

202.5 

0316*  0.495 

120 

225.0 

0317  *  0.495 

120 

247.5 

0318  *  0.495 

120 

270.0 

0319*  0.495 

120 

292.5 

0320  *  0.495 

120 

315.0 

0321  *  0.495 

120 

337.5 
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*  0322*  0.486  135  000.0 

*  0323  *  0.495  135  030.0 

*  0324  *  0.495  135  060.0 

*  0325  *  0.495  135  090.0 

*  0326  *  0.495  135  120.0 

*  0327  *  0.495  135  150.0 

*  0328  *  0.495  135  180.0 

*  0329  *  0.495  135  210.0 

*  0330  *  0.495  135  240.0 

*  0331  *  0.495  135  270.0 

*  0332  *  0.495  135  300.0 

*  0333  *  0.495  135  330.0 

*  0334  *  0.495  150  000.0 

*  0335  *  0.495  150  045.0 

*  0336  *  0.495  150  090.0 

*  0337  *  0.495  150  135.0 

*  0338  *  0.495  150  180.0 

*  0339  *  0.495  150  225.0 

*  0340  *  0.495  150  270.0 

*  0341  *  0.495  150  315.0 

*  0342  *  0.495  165  000.0 

*  0343  *  0.495  165  090.0 

*  0344  *  0.495  165  180.0 

*  0345  *  0.495  165  270.0 

*  0346  *  0.495  180  000.0 

*  0347  *  0.000  000  000.0 

*  0348*  0.000  000  000.0 

*  0349  *  0.000  000  000.0 

*  0350  *  0.000  000  000.0 

*  0351  *  0.000  000  000.0 

*  0352  *  0.000  000  000.0 

*  0353  *  0.000  000  000.0 

*  0354  *  0.000  000  000.0 

*  0355  *  0.000  000  000.0 

*  0356  *  0.000  000  000.0 

*  0357  *  0.000  000  000.0 

*  0358  *  0.000  000  000.0 

*  0359  *  0.000  000  000.0 

*  0360  *  0.000  000  000.0 

*  0361  *  0.000  000  000.0 

*  0362  *  0.000  000  000.0 

*  0363  *  0.000  000  000.0 

*  0364  *  0.000  000  000.0 

*  0365  *  0.000  000  000.0 

*  0366  *  0.000  000  000.0 

*  0367  *  0.000  000  000.0 

*  0368  *  0.000  000  000.0 

*  0369  *  0.000  000  000.0 

*  0370  *  0.000  000  000.0 

*  0371  *  0.000  000  000.0 

*  0372  *  0.000  000  000.0 
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*  0373  *  0.000  000  000.0 

*  0374*  0.000  000  000.0 

*  0375*  0.000  000  000.0 

*  0376  *  0.000  000  000.0 

*  0377  *  0.000  000  000.0 

*  0378  *  0.000  000  000.0 

*  0379*  0.000  000  000.0 

*  0380  *  0.000  000  000.0 

*  0381  *  0.000  000  000.0 

*  0382  *  0.000  000  000.0 

*  0383  *  0.000  000  000.0 

*  0384  *  0.000  000  000.0 

*  0385  *  0.000  000  000.0 

*  0386  *  0.000  000  000.0 

*  0387  *  0.000  000  000.0 

*  0388  *  0.000  000  000.0 

*  0389*  0.000  000  000.0 

*  0390  *  0.000  000  000.0 

*  0391  *  0.000  000  000.0 

*  0392*  0.000  000  000.0 

*  0393  *  0.000  000  000.0 

*  0394  *  0.000  000  000.0 

*  0395  *  0.000  000  000.0 

*  0396  *  0.000  000  000.0 

*  0397  *  0.000  000  000.0 

*  0398  *  0.000  000  0000 

*  0399*  0.000  000  000.0 

*  0400*  0.000  000  000.0 

*  0401  *  0.500  000  000.0 

*  0402  *  0.500  015  000.0 

*  0403  *  0.500  015  090.0 

*  0404  *  0.500  015  180.0 

*  0405  *  0.500  015  270.0 

*  0406  *  0.500  030  000.0 

*  0407  *  0.500  030  045.0 

*  0408  *  0.500  030  090.0 

*  0409  *  0.500  030  135.0 

*  0410  *  0.500  030  180.0 

*  0411  *  0.500  030  225.0 

*  0412  *  0.500  030  270.0 

*  0413  *  0.500  030  315.0 

*  0414  *  0.500  045  000.0 

*  0415  *  0.500  045  030.0 

*  0416  *  0.500  045  060.0 

*  0417  *  0.500  045  090.0 

*  0418  *  0.500  045  120.0 

*  0419  *  0.500  045  150.0 

*  0420  *  0.500  045  180.0 

*  0421  *  0.500  045  210.0 

*  0422  *  0.500  045  240.0 

*  0423  *  0.500  045  270.0 
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0424  *  0.500 

045 

300.0 

0425  *  0.500 

045 

330.0 

0426  *  0.500 

060 

000.0 

0427  *  0.500 

060 

022.5 

0428  *  0.500 

060 

045.0 

0429  *  0.500 

060 

067.5 

0430  *  0.500 

060 

090.0 

0431  *  0.500 

060 

112.5 

0432  *  0.500 

060 

135.0 

0433  *  0.500 

060 

157.5 

0434  *  0.500 

060 

180.0 

0435  *  0.500 

060 

202.5 

0436  *  0.500 

060 

225.0 

0437  *  0.500 

060 

247.5 

0438  *  0.500 

060 

270.0 

0439  *  0.500 

060 

292.5 

0440  *  0.500 

060 

315.0 

0441  *  0.500 

060 

337.5 

0442  *  0.500 

075 

000.0 

0443  *  0.500 

075 

018.0 

0444  *  0.500 

075 

036.0 

0445  *  0.500 

075 

054.0 

0446  *  0.500 

075 

072.0 

0447  *  0.500 

075 

090.0 

0448  *  0.500 

075 

108.0 

0449  *  0.500 

075 

126.0 

0450  *  0.500 

075 

144.0 

0451  *  0.500 

075 

162.0 

0452  *  0.500 

075 

180.0 

0453  *  0.500 

075 

198.0 

0454  *  0.500 

075 

216.0 

0455  *  0.500 

075 

234.0 

0456  *  0.500 

075 

252.0 

0457  *  0.500 

075 

270.0 

0458  *  0.500 

075 

288.0 

0459  *  0.500 

075 

306.0 

0460  *  0.500 

075 

324.0 

0461  *  0.500 

075 

342.0 

0462  *  0.500 

090 

000.0 

0463  *  0.500 

090 

015.0 

0464  *  0.500 

090 

030.0 

0465  *  0.500 

090 

045.0 

0466  *  0.500 

090 

060.0 

0467  *  0.500 

090 

075.0 

0468  *  0.500 

090 

090.0 

0469  *  0.500 

090 

105.0 

0470  *  0.500 

090 

120.0 

0471  *  0.500 

090 

135.0 

0472  *  0.500 

090 

150.0 

0473  *  0.500 

090 

165.0 

0474  *  0.500 

090 

180.0 

0475* 

0.500 

090 

195.0 

0476* 

0.500 

090 

210.0 

0477* 

0.500 

090 

225.0 

0478* 

0.500 

090 

240.0 

0479* 

0.500 

090 

255.0 

0480* 

0.500 

090 

270.0 

0481  * 

0.500 

090 

285.0 

0482* 

0.500 

090 

300.0 

0483* 

0.500 

090 

315.0 

0484* 

0.500 

090 

330.0 

0485* 

0.500 

090 

345.0 

0486* 

0.500 

105 

000.0 

0487* 

0.500 

105 

018.0 

0488* 

0.500 

105 

036.0 

0489* 

0.500 

105 

054.0 

0490* 

0.500 

105 

072.0 

0491  * 

0.500 

105 

090.0 

0492* 

0.500 

105 

108.0 

0493* 

0.500 

105 

126.0 

0494* 

0.500 

105 

144.0 

0495’ 

0.500 

105 

162.0 

0496* 

0.500 

105 

180.0 

0497* 

0.500 

105 

198.0 

0498* 

0.500 

105 

216.0 

0499’ 

0.500 

105 

234.0 

0500* 

0.500 

105 

252.0 

0501  * 

0.500 

105 

270.0 

0502* 

0.500 

105 

288.0 

0503* 

0.500 

105 

306.0 

0504* 

0.500 

105 

324.0 

0505* 

0.500 

105 

342.0 

0506* 

0.500 

120 

000.0 

0507* 

0.500 

120 

022.5 

0508* 

0.500 

120 

045.0 

0509* 

0.500 

120 

067.5 

0510* 

0.500 

120 

090.0 

0511  * 

0.500 

120 

112.5 

0512* 

0.500 

120 

135.0 

0513* 

0.500 

120 

157.5 

0514* 

0.500 

120 

180.0 

0515* 

0.500 

120 

202.5 

0516* 

0.500 

120 

225.0 

0517* 

0.500 

120 

247.5 

0518* 

0.500 

120 

270.0 

0519* 

0.500 

120 

292.5 

0520* 

0.500 

120 

315.0 

0521  * 

0.500 

120 

337.5 

0522* 

0.500 

135 

000.0 

0523* 

0.500 

135 

030.0 

0524* 

0.500 

135 

060.0 

0525* 

0.500 

135 

090.0 
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0526  * 

0.500 

135 

120.0 

0527* 

0.500 

135 

150.0 

0528* 

0.500 

135 

180.0 

0529* 

0.500 

135 

210.0 

0530* 

0.500 

135 

240.0 

0531  * 

0.500 

135 

270.0 

0532* 

0.500 

135 

300.0 

0533* 

0.500 

135 

330.0 

0534* 

0.500 

150 

000.0 

0535* 

0.500 

150 

045.0 

0536* 

0.500 

150 

090.0 

0537* 

0.500 

150 

135.0 

0538' 

0.500 

150 

180.0 

0539* 

0.500 

150 

225.0 

0540* 

0.500 

150 

270.0 

0541  * 

0.500 

150 

315.0 

0542* 

0.500 

165 

000.0 

0543* 

0.500 

165 

090.0 

0544* 

0.500 

165 

180.0 

0545* 

0.500 

165 

270.0 

0546* 

0.500 

180 

000.0 

0547* 

0.000 

000 

000.0 

0548* 

0.000 

000 

000.0 

0549* 

0.000 

000 

000.0 

0550* 

0.000 

000 

000.0 

0551  * 

0.000 

000 

000.0 

0552* 

0.000 

000 

000.0 

0553* 

0.000 

000 

000.0 

0554* 

0.000 

000 

000.0 

0555* 

0.000 

000 

000.0 

0556* 

0.000 

000 

000.0 

0557* 

0.000 

000 

000.0 

0558* 

0.000 

000 

000.0 

0559* 

0.000 

000 

000.0 

0560* 

0.000 

000 

000.0 

0561  * 

0.000 

000 

000.0 

0562* 

0.000 

000 

000.0 

0563* 

0.000 

000 

000.0 

0564* 

0.000 

000 

000.0 

0565* 

0.000 

000 

000.0 

0566* 

0.000 

000 

000.0 

0567* 

0.000 

000 

000.0 

0568* 

0.000 

000 

000.0 

0569* 

0.000 

000 

000.0 

0570* 

0.000 

000 

000.0 

0571  * 

0.000 

000 

000.0 

0572* 

0.000 

000 

000.0 

0573* 

0.000 

000 

000.0 

0574* 

0.000 

000 

000.0 

0575* 

0.000 

000 

000.0 

0576* 

0.000 

000 

000.0 
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*  0577  *  0.000  000  000.0 

*  0578  *  0.000  000  000.0 

*  0579  *  0.000  000  000.0 

*  0580  *  0.000  000  000.0 

*  0581  *  0.000  000  000.0 

*  0582  *  0.000  000  000.0 

*  0583  *  0.000  000  000.0 

*  0584  *  0.000  000  000.0 

*  0585  *  0.000  000  000.0 

*  0586  *  0.000  000  000.0 

*  0587  *  0.000  000  000.0 

*  0588  *  0.000  000  000.0 

*  0589  *  0.000  000  000.0 

*  0590  *  0.000  000  000.0 

*  0591  *  0.000  000  000.0 

*  0592  *  0.000  000  000.0 

*  0593  *  0.000  000  000.0 

*  0594  *  0.000  000  000.0 

*  0595  *  0.000  000  000.0 

*  0596  *  0.000  000  000.0 

*  0597  *  0.000  000  000.0 

*  0598  *  0.000  000  000.0 

*  0599  *  0.000  000  000.0 

*  0600  *  0.000  000  000.0 

*  0601  *  0.500  000  000.0 

*  0602  *  0.500  015  000.0 

*  0603  *  0.500  015  090.0 

*  0604  *  0.500  015  180.0 

*  0605  *  0.500  015  270.0 

*  0606  *  0.500  030  000.0 

*  0607  *  0.500  030  045.0 

*  0608  *  0.500  030  090.0 

*  0609  *  0.500  030  135.0 

*  0610  *  0.500  030  180.0 

*  0611  *  0.500  030  225.0 

*  0612  *  0.500  030  270.0 

*  0613  *  0.500  030  315.0 

*  0614  *  0.500  045  000.0 

*  0615  *  0.500  045  030.0 

*  0616  *  0.500  045  060.0 

*  0617  *  0.500  045  090.0 

*  0618  *  0.500  045  120.0 

*  0619  *  0.500  045  150.0 

*  0620  *  0.500  045  180.0 

*  0621  *  0.500  045  210.0 

*  0622  *  0.500  045  240.0 

*  0623  *  0.500  045  270.0 

*  0624  *  0.500  045  300.0 

*  0625  *  0.500  045  330.0 

*  0626  *  0.500  060  000.0 

*  0627  *  0.500  060  022.5 
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*  0628  *  0.500  060  045.0 

*  0629  *  0.500  060  067.5 

*  0630  *  0.500  060  090.0 

*  0631  *  0.500  060  112.5 

*  0632  *  0.500  060  135.0 

*  0633  *  0.500  060  157.5 

*  0634  *  0.500  060  180.0 

*  0635  *  0.500  060  202.5 

*  0636  *  0.500  060  225.0 

*  0637  *  0.500  060  247.5 

*  0638  *  0.500  060  270.0 

*  0639  *  0.500  060  292.5 

*  0640  *  0.500  060  315.0 

*  0641  *  0.500  060  337.5 

*  0642  *  0.500  075  000.0 

*  0643  *  0.500  075  018.0 

*  0644  *  0.500  075  036.0 

*  0645  *  0.500  075  054.0 

*  0646  *  0.500  075  072.0 

*  0647  *  0.500  075  090.0 

*  0648  *  0.500  075  108.0 

*  0649  *  0.500  075  126.0 

*  0650  *  0.500  075  144.0 

*  0651  *  0.500  075  162.0 

*  0652  *  0.500  075  180.0 

*  0653  *  0.500  075  198.0 

*  0654  *  0.500  075  216.0 

*  0655  *  0.500  075  234.0 

*  0656  *  0.500  075  252.0 

*  0657  *  0.500  075  270.0 

*  0658  *  0.500  075  288.0 

*  0659  *  0.500  075  306.0 

*  0660  *  0.500  075  324.0 

*  0661  *  0.500  075  342.0 

*  0662  *  0.500  090  000.0 

*  0663  *  0.500  090  015.0 

*  0664  *  0.500  090  030.0 

*  0665  *  0.500  090  045.0 

*  0666  *  0.500  090  060.0 

*  0667  *  0.500  090  075.0 

*  0668  *  0.500  090  090.0 

*  0669  *  0.500  090  105.0 

*  0670  *  0.500  090  120.0 

*  0671  *  0.500  090  135.0 

*  0672  *  0.500  090  150.0 

*  0673  *  0.500  090  165.0 

*  0674  *  0.500  090  180.0 

*  0675  *  0.500  090  195.0 

*  0676  *  0.500  090  210.0 

*  0677  *  0.500  090  225.0 

*  0678  *  0.500  090  240.0 


*  0679  *  0.500  090  255.0 

*  0680  *  0.500  090  270.0 

*  0681  *  0.500  090  285.0 

*  0682  *  0.500  090  300.0 

*  0683  *  0.500  090  315.0 

*  0684  *  0.500  090  330.0 

*  0685  *  0.500  090  345.0 

*  0686  *  0.500  105  000.0 

*  0687  *  0.500  105  018.0 

*  0688  *  0.500  105  036.0 

*  0689  *  0.500  105  054.0 

*  0690  *  0.500  105  072.0 

*  0691  *  0.500  105  090.0 

*  0692  *  0.500  105  108.0 

*  0693  *  0.500  105  126.0 

*  0694  *  0.500  105  144.0 

*  0695  *  0.500  105  162.0 

*  0696  *  0.500  105  180.0 

*  0697  *  0.500  105  198.0 

*  0698  *  0.500  105  216.0 

*  0699  *  0.500  105  234.0 

*  0700  *  0.500  105  252.0 

*  0701  *  0.500  105  270.0 

*  0702  *  0.500  105  288.0 

*  0703  *  0.500  105  306.0 

*  0704  *  0.500  105  324.0 

*  0705  *  0.500  105  342.0 

*  0706  *  0.500  120  000.0 

*  0707  *  0.500  120  022.5 

*  0708  *  0.500  120  045.0 

*  0709  *  0.500  120  067.5 

*  0710  *  0.500  120  090.0 

*  0711  *  0.500  120  112.5 

*  0712*  0.500  120  135.0 

*  0713*  0.500  120  157.5 

*  0714*  0.500  120  180.0 

*  0715  *  0.500  120  202.5 

*  0716  *  0.500  120  225.0 

*  0717*  0.500  120  247.5 

*  0718*  0.500  120  270.0 

*  0719  *  0.500  120  292.5 

*  0720  *  0.500  120  315.0 

*  0721  *  0.500  120  337.5 

*  0722  *  0.500  135  000.0 

*  0723  *  0.500  135  030.0 

*  0724  *  0.500  135  060.0 

*  0725*  0.500  135  090.0 

*  0726*  0.500  135  120.0 

*  0727  *  0.500  135  150.0 

*  0728  *  0.500  135  180.0 

*  0729  *  0.500  135  210.0 
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0730  *  0.500 

135 

240.0 

0731  *  0.500 

135 

270.0 

0732  *  0.500 

135 

300.0 

0733  *  0.500 

135 

330.0 

0734  *  0.500 

150 

000.0 

0735  *  0.500 

150 

045.0 

0736  *  0.500 

150 

090.0 

0737  *  0.500 

150 

135.0 

0738  *  0.500 

150 

180.0 

0739  *  0.500 

150 

225.0 

0740  *  0.500 

150 

270.0 

0741  *  0.500 

150 

315.0 

0742  *  0.500 

165 

000.0 

0743  *  0.500 

165 

090.0 

0744  *  0.500 

165 

180.0 

0745  *  0.500 

165 

270.0 

0746  *  0.500 

180 

000.0 

0747  *  0.000 

000 

000.0 

0748  *  0.000 

000 

000.0 

0749  *  0.000 

000 

000.0 

0750  *  0.000 

000 

000.0 

0751  *  0.000 

000 

000.0 

0752  *  0.000 

000 

000.0 

0753  *  0.000 

000 

000.0 

0754  *  0.000 

000 

000.0 

0755  *  0.000 

000 

000.0 

0756  *  0.000 

000 

000.0 

0757  *  0.000 

000 

000.0 

0758  *  0.000 

000 

000.0 

0759  *  0.000 

000 

000.0 

0760  *  0.000 

000 

000.0 

0761  *  0.000 

000 

000.0 

0762  *  0.000 

000 

000.0 

0763  *  0.000 

000 

000.0 

0764  *  0.000 

000 

000.0 

0765  *  0.000 

000 

000.0 

0766  *  0.000 

000 

000.0 

0767  *  0.000 

000 

000.0 

0768  *  0.000 

000 

000.0 

0769  *  0.000 

000 

000.0 

0770  *  0.000 

000 

000.0 

0771  *  0.000 

000 

000.0 

0772  *  0.000 

000 

000.0 

0773  *  0.000 

000 

000.0 

0774  *  0.000 

000 

000.0 

0775  *  0.000 

000 

000.0 

0776  *  0.000 

000 

000.0 

0777*  0.000 

000 

000.0 

0778  *  0.000 

000 

000.0 

0779  *  0.000 

000 

000.0 

0780  *  0.000 

000 

000.0 

0781  * 

0.000 

000 

000.0 

0782* 

0.000 

000 

000.0 

0783* 

0.000 

000 

000.0 

0784* 

0.000 

000 

000.0 

0785* 

0.000 

000 

000.0 

0786* 

0.000 

000 

000.0 

0787* 

0.000 

000 

000.0 

0788* 

0.000 

000 

000.0 

0789* 

0.000 

000 

000.0 

0790* 

0.000 

000 

000.0 

0791  * 

0.000 

000 

000.0 

0792* 

0.000 

000 

000.0 

0793* 

0.000 

000 

000.0 

0794* 

0.000 

000 

000.0 

0795* 

0.000 

000 

000.0 

0796* 

0.000 

000 

000.0 

0797* 

0.000 

000 

000.0 

0798* 

0.000 

000 

000.0 

0799* 

0.000 

000 

000.0 

0800* 

0.000 

000 

000.0 

0801  * 

0.625 

000 

000.0 

0802* 

0.625 

015 

000.0 

0803* 

0.625 

015 

090.0 

0804* 

0.625 

015 

180.0 

0805* 

0.625 

015 

270.0 

0806* 

0.625 

030 

000.0 

0807* 

0.625 

030 

045.0 

0808* 

0.625 

030 

090.0 

0809* 

0.625 

030 

135.0 

0810* 

0.625 

030 

180.0 

0811  * 

0.625 

030 

225.0 

0812* 

0.625 

030 

270.0 

0813* 

0.625 

030 

315.0 

0814* 

0.625 

045 

000.0 

0815* 

0.625 

045 

030.0 

0816* 

0.625 

045 

060.0 

0817* 

0.625 

045 

090.0 

0818* 

0.625 

045 

120.0 

0819* 

0.625 

045 

150.0 

0820* 

0.625 

045 

180.0 

0821  * 

0.625 

045 

210.0 

0822* 

0.625 

045 

240.0 

0823* 

0.625 

045 

270.0 

0824* 

0.625 

045 

300.0 

0825* 

0.625 

045 

330.0 

0826* 

0.625 

060 

000.0 

0827* 

0.625 

060 

022.5 

0828* 

0.625 

060 

045.0 

0829* 

0.625 

060 

067.5 

0830* 

0.625 

060 

090.0 

0831  * 

0.625 

060 

112.5 
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0832  *  0.625 

060 

135.0 

0833  *  0.625 

060 

157.5 

0834  *  0.625 

060 

180.0 

0835  *  0.625 

060 

202.5 

0836  *  0.625 

060 

225.0 

0837  *  0.625 

060 

247.5 

0838  *  0.625 

060 

270.0 

0839  *  0  325 

060 

292.5 

0840  *  0.625 

060 

315.0 

0841  *  0.625 

060 

337.5 

0842  *  0.625 

075 

000.0 

0843  *  0.625 

075 

018.0 

0844  *  0.625 

075 

036.0 

0845  *  0.625 

075 

054.0 

0846  *  0.625 

075 

072.0 

0847  *  0.625 

075 

090.0 

0848  *  0.625 

075 

108.0 

0849  *  0.625 

075 

126.0 

0850  *  0.625 

075 

144.0 

0851  *  0.625 

075 

162.0 

0852  *  0.625 

075 

180.0 

0853  *  0.625 

075 

198.0 

0854  *  0.625 

075 

218.0 

0855  *  0.625 

075 

234.0 

0856  *  0.625 

075 

252.0 

0857  *  0.625 

075 

270.0 

0858  *  0.625 

075 

288.0 

0859  *  0.625 

075 

306.0 

0860  *  0.625 

075 

324.0 

0861  *  0.625 

075 

342.0 

0862  *  0.625 

090 

000.0 

0863  *  0.625 

090 

015.0 

0864  *  0.625 

090 

030.0 

0865*  0.625 

090 

045.0 

0866*  0.625 

090 

060.0 

0867  *  0.625 

090 

075.0 

0868  *  0.625 

090 

090.0 

0869*  0.625 

090 

105.0 

0870  *  0.625 

090 

120.0 

0871  *  0.625 

090 

135.0 

0872  *  0.625 

090 

150.0 

0873*  0.625 

090 

165.0 

0874  *  0.625 

090 

180.0 

0875*  0.625 

090 

195.0 

0876*  0.625 

090 

210.0 

0877  *  0.625 

090 

225.0 

0878*  0.625 

090 

240.0 

0879*  0.625 

090 

255.0 

0880*  0.625 

090 

270.0 

0881  *  0.625 

090 

285.0 

0882*  0.625 

090 

300.0 
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0683* 

0.625 

090 

315.0 

0884* 

0.625 

090 

330.0 

0885* 

0.625 

090 

345.0 

0886* 

0.625 

105 

000.0 

0887* 

0.625 

105 

018.0 

0888* 

0.625 

105 

036.0 

0889* 

0.625 

105 

054.0 

0890* 

0.625 

105 

072.0 

0891  * 

0.625 

105 

090.0 

0892* 

0.625 

105 

108.0 

0893* 

0.625 

105 

126.0 

0894* 

0.625 

105 

144.0 

0895* 

0.625 

105 

162.0 

0896* 

0.625 

105 

160.0 

0897* 

0.625 

105 

198.0 

0898* 

0.625 

105 

216.0 

0899* 

0.625 

105 

234.0 

0900* 

0.625 

105 

252.0 

0901  * 

0.625 

105 

270.0 

0902* 

0.625 

105 

288.0 

0903* 

0.625 

105 

306.0 

0904* 

0.625 

105 

324.0 

0905* 

0.625 

105 

342.0 

0906* 

0.625 

120 

000.0 

0907* 

0.625 

120 

022.5 

0908* 

0.625 

120 

045.0 

0909* 

0.625 

120 

067.5 

0910* 

0.625 

120 

090.0 

0911  * 

0.625 

120 

112.5 

0912* 

0.625 

120 

135.0 

0913* 

0.625 

120 

157.5 

0914* 

0.625 

120 

180.0 

0915* 

0.625 

120 

202.5 

0916* 

0.625 

120 

225.0 

0917* 

0.625 

120 

247.5 

0918* 

0.625 

120 

270.0 

0919* 

0.625 

120 

292.5 

0920* 

0.625 

120 

315.0 

0921  * 

0.625 

120 

337.5 

0922* 

0.625 

135 

000.0 

0923* 

0.625 

135 

030.0 

0924* 

0.625 

135 

060.0 

0925* 

0.625 

135 

090.0 

0926* 

0.625 

135 

120.0 

0927* 

0.625 

135 

150.0 

0928* 

0.625 

135 

180.0 

0929* 

0.625 

135 

210.0 

0930* 

0.625 

135 

240.0 

0931  * 

0.625 

135 

270.0 

0932* 

0.625 

135 

300.0 

0933* 

0.625 

135 

330.0 
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0934* 

0.625 

150 

000.0 

0935* 

0.625 

150 

045.0 

0936* 

0.625 

150 

090.0 

0937* 

0.625 

150 

135.0 

0938* 

0.625 

150 

180.0 

0939* 

0.625 

150 

225.0 

0940* 

0.625 

150 

270.0 

0941  * 

0.625 

150 

315.0 

0942* 

0.625 

165 

000.0 

0943* 

0.625 

165 

090.0 

0944* 

0.625 

165 

180.0 

0945* 

0.625 

165 

270.0 

0946* 

0.625 

180 

000.0 

0947* 

0.000 

000 

000.0 

0948* 

0.000 

000 

000.0 

0949* 

0.000 

000 

000.0 

0950* 

0.000 

000 

000.0 

0951  * 

0.000 

000 

000.0 

0952* 

0.000 

000 

000.0 

0953* 

0.000 

000 

000.0 

0954* 

0.000 

000 

000.0 

0955* 

0.000 

000 

000.0 

0956’ 

0.000 

000 

000.0 

0957* 

0.000 

000 

000.0 

0958* 

0.000 

000 

000.0 

0959* 

0.000 

000 

000.0 

0960* 

0.000 

000 

000.0 

0961  * 

0.000 

000 

000.0 

0962* 

0.000 

000 

000.0 

0963* 

0.000 

000 

000.0 

0964* 

0.000 

000 

000.0 

0965* 

0.000 

000 

000.0 

0966* 

0.000 

000 

000.0 

0967* 

0.000 

000 

000.0 

0968* 

0.000 

000 

000.0 

0969* 

0.000 

000 

000.0 

0970* 

0.000 

000 

000.0 

0971  * 

0.000 

000 

000.0 

0972’ 

0.000 

000 

000.0 

0973* 

0.000 

000 

000.0 

0974* 

0.000 

000 

000.0 

0975* 

0.000 

000 

000.0 

0976* 

0.000 

000 

000.0 

0977* 

0.000 

000 

000.0 

0978* 

0.000 

000 

000.0 

0979* 

0.000 

000 

000.0 

0980* 

0.000 

000 

000.0 

0981  * 

0.000 

oco 

000.0 

0982* 

0.000 

000 

000.0 

0983* 

0.000 

000 

000.0 

0984* 

0  000 

000 

000.0 
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0965  *  0.000 
0966  *  0.000 
0987  *  0.000 
0968  *  0.000 


0991  * 


0995  *  0.000 

000 

000.0 

0996  *  0.000 

000 

000.0 

0997  *  0.000 

000 

000.0 

0996  *  0.000 

000 

000.0 

0999  *  0.000 

000 

000.0 

1000  *  0.000 

000 

000.0 

1001  *  0.750 

000 

000.0 

1002*  0.750 

015 

000.0 

1003  *  0.750 

015 

090.0 

1004  *  0.750 

015 

180.0 

1005  *  0.750 

015 

270.0 

1006  *  0.750 

030 

000.0 

1007  *  0.750 

030 

045.0 

1008  *  0.750 

030 

090.0 

1009  *  0.750 

030 

135.0 

1010*  0.750 

030 

180.0 

1011*  0.750 

030 

225.0 

1012*  0.750 

030 

270.0 

1013*  0.750 

030 

315.0 

1014*  0.750 

045 

000.0 

1015*  0.750 

045 

030.0 

1016*  0.750 

045 

060.0 

1017*  0.750 

045 

090.0 

1018*  0.750 

045 

120.0 

1019*  0.750 

045 

150.0 

1020  *  0.750 

045 

180.0 

1021  *  0.750 

045 

210.0 

1022  *  0.750 

045 

240.0 

1023*  0.750 

045 

270.0 

1024  *  0.750 

045 

300.0 

1025*  0.750 

045 

330.0 

1026  *  0.750 

060 

000.0 

1027  *  0.750 

060 

022.5 

1028  *  0.750 

060 

045.0 

1029  *  0.750 

060 

067.5 

1030  *  0.750 

060 

090.0 

1031  *  0.750 

060 

112.5 

1032  *  0.750 

060 

135.0 

1033  *  0.750 

060 

157.5 

1034  *  0.750 

060 

180.0 

1035*  0.750 

060 

202.5 

*  1036  *  0.750  060  225.0 

*  1037  *  0.750  060  247.5 

*  1038  *  0.750  060  270.0 

*  1039  *  0.750  060  292.5 

*  1040  *  0.750  060  315.0 

*  1041  *  0.750  060  337.5 

*  1042  *  0.750  075  000.0 

*  1043  *  0.750  075  018.0 

*  1044  *  0.750  075  036.0 

*  1045  *  0.750  075  054.0 

*  1046  *  0.750  075  072.0 

*  1047  *  0.750  075  090.0 

*  1048  *  0.750  075  108.0 

*  1049  *  0.750  075  126.0 

*  1050  *  0.750  075  144.0 

*  1051  *  0.750  075  162.0 

*  1052  *  0.750  075  180.0 

*  1053  *  0.750  075  198.0 

*  1054  *  0.750  075  216.0 

*  1055  *  0.750  075  234.0 

*  1056  *  0.750  075  252.0 

*  1057  *  0.750  075  270.0 

*  1058  *  0.750  075  288.0 

*  1059  *  0.750  075  306.0 

*  1060  *  0.750  075  324.0 

*  1061  *  0.750  075  342.0 

*  1062  *  0.750  090  000.0 

*  1063  *  0.750  090  015.0 

*  1064  *  0.750  090  030.0 

*  1065  *  0.750  090  045.0 

*  1066  *  0.750  090  060.0 

*  1067  *  0.750  090  075.0 

*  1068  *  0.750  090  090.0 

*  1069  *  0.750  090  105.0 

*  1070  *  0.750  090  120.0 

*  1071  *  0.750  090  135.0 

*  1072  *  0.750  090  150.0 

*  1073  *  0.750  090  165.0 

*  1074  *  0.750  090  180.0 

*  1075  *  0.750  090  195.0 

*  1076  *  0.750  090  210.0 

*  1077  *  0.750  090  225.0 

*  1078  *  0.750  090  240.0 

*  1079  *  0.750  090  255.0 

*  1080  *  0.750  090  270.0 

*  1081  *  0.750  090  285.0 

*  1082  *  0.750  090  300.0 

*  1083  *  0.750  090  315.0 

*  1084  *  0.750  090  330.0 

*  1085  *  0.750  090  345.0 

*  1086  *  0.750  105  000.0 
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r 


1087* 

0.750 

105 

018.0 

1088* 

0.750 

105 

036.0 

1089* 

0.750 

105 

054.0 

1090* 

0.750 

105 

072.0 

1091  * 

0.750 

105 

090.0 

1092* 

0.750 

105 

108.0 

1093* 

0.750 

105 

126.0 

1094* 

0.750 

105 

144.0 

1095* 

0.750 

105 

162.0 

1096* 

0.750 

105 

180.0 

1097* 

0.750 

105 

198.0 

1098* 

0.750 

105 

216.0 

1099* 

0.750 

105 

234.0 

1100* 

0.750 

105 

252.0 

1101  * 

0.750 

105 

270.0 

1102* 

0.750 

105 

288.0 

1103* 

0.750 

105 

306.0 

1104* 

0.750 

105 

324.0 

1105* 

0.750 

105 

342.0 

1106* 

0.750 

120 

000.0 

1107* 

0.750 

120 

022.5 

1108* 

0.750 

120 

045.0 

1109* 

0.750 

120 

067.5 

1110* 

0.750 

120 

090.0 

1111  * 

0.750 

120 

112.5 

1112* 

0.750 

120 

135.0 

1113* 

0.750 

120 

157.5 

1114* 

0.750 

120 

180.0 

1115* 

0.750 

120 

202.5 

1116* 

0.750 

120 

225.0 

1117* 

0.750 

120 

247.5 

1118* 

0.750 

120 

270.0 

1119* 

0.750 

120 

292.5 

1120* 

0.750 

120 

315.0 

1121  * 

0.750 

120 

337.5 

1122* 

0.750 

135 

000.0 

1123* 

0.750 

135 

030.0 

1124* 

0.750 

135 

060.0 

1125* 

0.750 

135 

090.0 

1126* 

0.750 

135 

120.0 

1127* 

0.750 

135 

150.0 

1128* 

0.750 

135 

180.0 

1129* 

0.750 

135 

210.0 

1130* 

0.750 

135 

240.0 

1131  * 

0.750 

135 

270.0 

1132* 

0.750 

135 

300.0 

1133* 

0.750 

135 

330.0 

1134* 

0.750 

150 

000.0 

1135* 

0.750 

150 

045.0 

1136* 

0.750 

150 

090.0 

1137* 

0.750 

150 

135.0 
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*  1138  *  0.750  150  180.0 

*  1139  *  0.750  150  225.0 

*  1140  *  0.750  150  270.0 

*  1141  *  0.750  150  315.0 

*  1142  *  0.750  165  000.0 

*  1143  *  0.750  165  090.0 

*  1144  *  0.750  165  180.0 

*  1145  *  0.750  165  270.0 

*  1146  *  0.750  180  000.0 

*  1147  *  0.000  000  000.0 

*  1148  *  0.000  000  000.0 

*  1149  *  0.000  000  000.0 

*  1150  *  0.000  000  000.0 

*  1151  *  0.000  000  000.0 

*  1152  *  0.000  000  000.0 

*  1153  *  0.000  000  000.0 

*  1154  *  0.000  000  000.0 

*  1155  *  0.000  000  000.0 

*  1156  *  0.000  000  000.0 

*  1157  *  0.000  000  000.0 

*  1158  *  0.000  000  000.0 

*  1159  *  0.000  000  000.0 

*  1160  *  0.000  000  000.0 
*  1161  *  0.000  000  000.0 
*  1162  *  0.000  000  000.0 

*  1163  *  0.000  000  000.0 

*  1164*  0.000  000  000.0 

*  1165  *  0.000  000  000.0 

*  1166  *  0.000  000  000.0 

*  1167  *  0.000  000  000.0 

*  1168  *  0.000  000  000.0 

*  1169  *  0.000  000  000.0 

*  1170  *  0.000  000  000.0 

*  1171  *  0.000  000  000.0 

*  1172  *  0.000  000  000.0 

*  1173  *  0.000  000  000.0 

*  1174  *  0.000  000  000.0 

*  1175  *  0.000  000  000.0 

*  1176  *  0.000  000  000.0 

*  1177  *  0.000  000  000.0 

*  1178  *  0.000  000  000.0 

*  1179  *  0.000  000  000.0 

*  1180  *  0.000  000  000.0 
*  1181  *  0.000  000  000.0 
*  1182  *  0.000  000  000.0 

*  1183  *  0.000  000  000.0 

*  1184  *  0.000  000  000.0 

*  1185  *  0.000  000  000.0 

*  1186  *  0.000  000  000.0 

*  1187  *  0.000  000  000.0 

*  1188  *  0.000  000  000.0 
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1189* 

0.000 

000 

000.0 

1190* 

0.000 

000 

000.0 

1191  * 

0.000 

000 

000.0 

1192* 

0.000 

000 

000.0 

1193* 

0.000 

000 

000.0 

1194* 

0.000 

000 

000.0 

1195* 

0.000 

000 

000.0 

1196* 

0.000 

000 

000.0 

1197* 

0.000 

000 

000.0 

1198* 

0.000 

000 

000.0 

1199* 

0.000 

000 

000.0 

1200* 

0.000 

000 

000.0 

1201  * 

0.875 

000 

000.0 

1202* 

0.875 

015 

000.0 

1203* 

0.875 

015 

090.0 

1204* 

0.875 

015 

180.0 

1205* 

0.875 

015 

270.0 

1206* 

0.875 

030 

000.0 

1207* 

0.875 

030 

045.0 

1208* 

0.875 

030 

090.0 

1209* 

0.875 

030 

135.0 

1210* 

0.875 

030 

180.0 

1211  * 

0.875 

030 

225.0 

1212* 

0.875 

030 

270.0 

1213* 

0.875 

030 

315.0 

1214* 

0.875 

045 

000.0 

1215* 

0.875 

045 

030.0 

1216* 

0.875 

045 

060.0 

1217* 

0.875 

045 

090.0 

1218* 

0.875 

045 

120.0 

1219* 

0.875 

045 

150.0 

1220* 

0.875 

045 

180.0 

1221  * 

0.875 

045 

210.0 

1222* 

0.875 

045 

240.0 

1223* 

0.875 

045 

270.0 

1224* 

0.875 

045 

300.0 

1225* 

0.875 

045 

330.0 

1226* 

0.875 

060 

000.0 

1227* 

0.875 

060 

022.5 

1228* 

0.875 

060 

045.0 

1229* 

0.875 

060 

067.5 

1230* 

0.875 

060 

090.0 

1231  * 

0.875 

060 

112.5 

1232* 

0.875 

060 

135.0 

1233* 

0.875 

060 

157.5 

1234* 

0.875 

060 

180.0 

1235* 

0.875 

060 

202.5 

1236* 

0.875 

060 

225.0 

1237* 

0.875 

060 

247.5 

1238* 

0.875 

060 

270.0 

1239* 

0.875 

060 

292.5 
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1240*  0.875 

060 

315.0 

1241*  0.875 

060 

337.5 

1242  *  0.875 

075 

000.0 

1243  *  0.875 

075 

018.0 

1244  *  0.875 

075 

036.0 

1245  *  0.875 

075 

054.0 

1246  *  0.875 

075 

072.0 

1247  *  0.875 

075 

090.0 

1248  *  0.875 

075 

108.0 

1249  *  0.875 

075 

126.0 

1250*  0.875 

075 

144.0 

1251  *  0.875 

075 

162.0 

1252  *  0.875 

075 

180.0 

1253  *  0.875 

075 

198.0 

1254  *  0.875 

075 

216.0 

1255  *  0.875 

075 

234.0 

1256  *  0.875 

075 

252.0 

1257  *  0.875 

075 

270.0 

1258  *  0.875 

075 

288.0 

1259  *  0.875 

075 

306.0 

1260  *  0.875 

075 

324.0 

1261  *  0.875 

075 

342.0 

1262  *  0.875 

090 

000.0 

1263  *  0.875 

090 

015.0 

1264*  0.875 

090 

030.0 

1265  *  0.875 

090 

045.0 

1266*  0.875 

090 

060.0 

1267*  0.875 

090 

075.0 

1268  *  0.875 

090 

090.0 

1269  *  0.875 

090 

105.0 

1270  *  0.875 

090 

120.0 

1271  *  0.875 

090 

135.0 

1272*  0.875 

090 

150.0 

1273  *  0.875 

090 

165.0 

1274  *  0.875 

090 

180.0 

1275  *  0.875 

090 

195.0 

1276  *  0.875 

090 

210.0 

1277  *  0.875 

090 

225.0 

1278  *  0.875 

090 

240.0 

1279  *  0.875 

090 

255.0 

1280  *  0.875 

090 

270.0 

1281  *  0.875 

090 

285.0 

1282*  0.875 

090 

300.0 

1283*  0.875 

090 

315.0 

1284  *  0.875 

090 

330.0 

1285  *  0.875 

090 

345.0 

1286*  0.875 

105 

000.0 

1287  *  0.875 

105 

018.0 

1288  *  0.875 

105 

036.0 

1289  *  0.875 

105 

054.0 

1290*  0.875 

105 

072.0 
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1291  *  0.875 

105 

090.0 

1292  *  0.875 

105 

108.0 

1293  *  0.875 

105 

126.0 

1294  *  0.875 

105 

144.0 

1295  *  0.875 

105 

162.0 

1296  *  0.875 

105 

180.0 

1297  *  0.875 

105 

198.0 

1298*  0.875 

105 

216.0 

1299  *  0.875 

105 

234.0 

1300  *  0.875 

105 

252.0 

1301  *  0.875 

105 

270.0 

1302*  0.875 

105 

288.0 

1303  *  0.875 

105 

306.0 

1304  *  0.875 

105 

324.0 

1305  *  0.875 

105 

342.0 

1306  *  0.875 

120 

000.0 

1307  *  0.875 

120 

022.5 

1308  *  0.875 

120 

045.0 

1309  *  0.875 

120 

067.5 

1310*  0.875 

120 

090.0 

1311  *  0.875 

120 

112.5 

1312*  0.875 

120 

135.0 

1313*  0.875 

120 

157.5 

1314*  0.875 

120 

180.0 

1315*  0.875 

120 

202.5 

1316*  0.875 

120 

225.0 

1317*  0.875 

120 

247.5 

1318*  0.875 

120 

270.0 

1319*  0.875 

120 

292.5 

1320  *  0.875 

120 

315.0 

1321  *  0.875 

120 

337.5 

1322  *  0.875 

135 

000.0 

1323  *  0.875 

135 

030.0 

1324  *  0.875 

135 

060.0 

1325  *  0.875 

135 

090.0 

1326  *  0.875 

135 

120.0 

1327  *  0.875 

135 

150.0 

1328*  0.875 

135 

180.0 

1329  *  0.875 

135 

210.0 

1330  *  0.875 

135 

240.0 

1331  *  0.875 

135 

270.0 

1332*  0.875 

135 

300.0 

1333  *  0.875 

135 

330.0 

1334  *  0.875 

150 

000.0 

1335  *  0.875 

150 

045.0 

1336  *  0.875 

150 

090.0 

1337*  0.875 

150 

135.0 

1338  *  0.875 

150 

180.0 

1339*  0.875 

150 

225.0 

1340  *  0.875 

150 

270.0 

1341  *  0.875 

150 

315.0 
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1342* 

0.875 

165 

000.0 

1343* 

0.875 

165 

090.0 

1344* 

0.875 

165 

180.0 

1345* 

0.875 

165 

270.0 

1346* 

0.875 

180 

000.0 

1347* 

0.000 

000 

000.0 

1348* 

0.000 

000 

000.0 

1349* 

0.000 

000 

000.0 

1350* 

0.000 

000 

000.0 

1351  * 

0.000 

000 

000.0 

1352* 

0.000 

000 

000.0 

1353* 

0.000 

000 

000.0 

1354* 

0.000 

000 

000.0 

1355* 

0.000 

000 

000.0 

1356* 

0.000 

000 

000.0 

1357* 

0.000 

000 

000.0 

1358* 

0.000 

000 

000.0 

1359* 

0.000 

000 

000.0 

1360* 

0.000 

000 

000.0 

1361  * 

0.000 

000 

000.0 

1362* 

0.000 

000 

000.0 

1363* 

0.000 

000 

000.0 

1364* 

0.000 

000 

000.0 

1365* 

0.000 

000 

000.0 

1366* 

0.000 

000 

000.0 

1367* 

0.000 

000 

000.0 

1368* 

0.000 

000 

000.0 

1369* 

0.000 

000 

000.0 

1370* 

0.000 

000 

000.0 

1371  * 

0.000 

000 

000.0 

1372* 

0.000 

000 

000.0 

1373* 

0.000 

000 

000.0 

1374* 

0.000 

000 

000.0 

1375* 

0.000 

000 

000.0 

1376* 

0.000 

000 

000.0 

1377* 

0.000 

000 

000.0 

1378* 

0.000 

000 

000.0 

1379* 

0.000 

000 

000.0 

1380* 

0.000 

000 

000.0 

1381  * 

0.000 

000 

000.0 

1382* 

0.000 

000 

000.0 

1383* 

0.000 

000 

000.0 

1384* 

0.000 

000 

000.0 

1385* 

0.000 

000 

000.0 

1386* 

0.000 

000 

000.0 

1387* 

0.000 

000 

000.0 

1388* 

0.000 

000 

000.0 

1389* 

0.000 

000 

000.0 

1390* 

0.000 

000 

000.0 

1391  * 

0.000 

000 

000.0 

1392* 

0.000 

000 

000.0 

1393* 

0.000 

000 

000.0 

1394* 

0.000 

000 

000.0 

1395* 

0.000 

000 

000.0 

1396* 

0.000 

000 

000.0 

1397* 

0.000 

000 

000.0 

1398* 

0.000 

000 

000.0 

1399* 

0.000 

000 

000.0 

1400* 

0.000 

000 

000.0 

1401  * 

1.000 

000 

000.0 

1402* 

1.000 

015 

000.0 

1403* 

1.000 

015 

090.0 

1404* 

1.000 

015 

180.0 

1405* 

1.000 

015 

270.0 

1406* 

1.000 

030 

000.0 

1407* 

1.000 

030 

045.0 

1408* 

1.000 

030 

090.0 

1409* 

1.000 

030 

135.0 

1410* 

1.000 

030 

180.0 

1411  * 

1.000 

030 

225.0 

1412* 

1.000 

030 

270.0 

1413* 

1.000 

030 

315.0 

1414* 

1.000 

045 

000.0 

1415* 

1.000 

045 

030.0 

1416* 

1.000 

045 

060.0 

1417* 

1.000 

045 

090.0 

1418* 

1.000 

045 

120.0 

1419* 

1.000 

045 

150.0 

1420* 

1.000 

045 

180.0 

1421  * 

1.000 

045 

210.0 

1422* 

1.000 

045 

240.0 

1423* 

1.000 

045 

270.0 

1424* 

1.000 

045 

300.0 

1425* 

1.000 

045 

330.0 

1426* 

1.000 

060 

000.0 

1427* 

1.000 

060 

022.5 

1428* 

1.000 

060 

045.0 

1429* 

1.000 

060 

067.5 

1430* 

1.000 

060 

090.0 

1431  * 

1.000 

060 

112.5 

1432* 

1.000 

060 

135.0 

1433* 

1.000 

060 

157.5 

1434* 

1.000 

060 

190.0 

1435* 

1.000 

060 

202.5 

1436  * 

1.000 

060 

225.0 

1437* 

1.000 

060 

247.5 

1438* 

1.000 

060 

270.0 

1439* 

1.000 

060 

292.5 

1440* 

1.000 

060 

315.0 

1441  * 

1.000 

060 

337.5 

1442* 

1.000 

075 

000.0 

1443* 

1.000 

075 

018.0 
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1444* 

1.000 

075 

036.0 

1445* 

1.000 

075 

054.0 

1446* 

1.000 

075 

072.0 

1447* 

1.000 

075 

090.0 

1448* 

1.000 

075 

108.0 

1449* 

1.000 

075 

126.0 

1450* 

1.000 

075 

144.0 

1451  * 

1.000 

075 

162.0 

1452* 

1.000 

075 

180.0 

1453* 

1.000 

075 

198.0 

1454* 

1.000 

075 

216.0 

1455* 

1.000 

075 

234.0 

1456* 

1.000 

075 

252.0 

1457* 

1.000 

075 

270.0 

1458* 

1.000 

075 

288.0 

1459* 

1.000 

075 

306.0 

1460* 

1.000 

075 

324.0 

1461  * 

1.000 

075 

342.0 

1462* 

1.000 

090 

000.0 

1463* 

1.000 

090 

015.0 

1464* 

1.000 

090 

030.0 

1465* 

1.000 

090 

045.0 

1466* 

1.000 

090 

060.0 

1467* 

1.000 

090 

075.0 

1468* 

1.000 

090 

090.0 

1469* 

1.000 

090 

105.0 

1470* 

1.000 

090 

120.0 

1471  * 

1.000 

090 

135.0 

1472* 

1.000 

090 

150.0 

1473* 

1.000 

090 

165.0 

1474* 

1.000 

090 

180.0 

1475* 

1.000 

090 

195.0 

1476* 

1.000 

090 

210.0 

1477* 

1.000 

090 

225.0 

1478* 

1.000 

090 

240.0 

1479* 

1.000 

090 

255.0 

1480* 

1.000 

090 

270.0 

1481  * 

1.000 

090 

285.0 

1482* 

1.000 

090 

300.0 

1483* 

1.000 

090 

315.0 

1484* 

1.000 

090 

330.0 

1485* 

1.000 

090 

345.0 

1486* 

1.000 

105 

000.0 

1487* 

1.000 

105 

018.0 

1488* 

1.000 

105 

036.0 

1489* 

1.000 

105 

054.0 

1490* 

1.000 

105 

072.0 

1491  * 

1.000 

105 

090.0 

1492* 

1.000 

105 

108.0 

1493* 

1.000 

105 

126.0 

1494* 

1.000 

105 

144.0 
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1495* 

1.000 

105 

162.0 

1496* 

1.000 

105 

180.0 

1497* 

1.000 

105 

198.0 

1496* 

1.000 

105 

216.0 

1499* 

1.000 

105 

234.0 

1500* 

1.000 

105 

252.0 

1501  * 

1.000 

105 

270.0 

1502* 

1.000 

105 

288.0 

1503* 

1.000 

105 

306.0 

1504* 

1.000 

105 

324.0 

1505* 

1.000 

105 

342.0 

1506* 

1.000 

120 

000.0 

1507* 

1.000 

120 

022.5 

1508* 

1.000 

120 

045.0 

1509* 

1.000 

120 

067.5 

1510* 

1.000 

120 

090.0 

1511  * 

1.000 

120 

112.5 

1512* 

1.000 

120 

135.0 

1513* 

1.000 

120 

157.5 

1514* 

1.000 

120 

180.0 

1515* 

1.000 

120 

202.5 

1516* 

1.000 

120 

225.0 

1517* 

1.000 

120 

247.5 

1518* 

1.000 

120 

270.0 

1519* 

1.000 

120 

292.5 

1520* 

1.000 

120 

315.0 

1521  * 

1.000 

120 

337.5 

1522* 

1.000 

135 

000.0 

1523* 

1.000 

135 

030.0 

1524* 

1.000 

135 

060.0 

1525* 

1.000 

135 

090.0 

1526* 

1.000 

135 

120.0 

1527* 

1.000 

135 

150.0 

1528* 

1.000 

135 

180.0 

1529* 

1.000 

135 

210.0 

1530* 

1.000 

135 

240.0 

1531  * 

1.000 

135 

270.0 

1532* 

1.000 

135 

300.0 

1533* 

1.000 

135 

330.0 

1534* 

1.000 

150 

000.0 

1535* 

1.000 

150 

045.0 

1536* 

1.000 

150 

090.0 

1537* 

1.000 

150 

135.0 

1538* 

1.000 

150 

180.0 

1539* 

1.000 

150 

225.0 

1540* 

1.000 

150 

270.0 

1541  * 

1.000 

150 

315.0 

1542* 

1.000 

165 

000.0 

1543* 

1.000 

165 

090.0 

1544* 

1.000 

165 

180.0 

1545* 

1.000 

165 

270.0 
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1546* 

1.000 

180 

000.0 

1547* 

0.000 

000 

000.0 

1548* 

0.000 

000 

000.0 

1549* 

0.000 

000 

000.0 

1550* 

0.000 

000 

000.0 

1551  * 

0.000 

000 

000.0 

1552* 

0.000 

000 

000.0 

1553* 

0.000 

000 

000.0 

1554* 

0.000 

000 

000.0 

1555* 

0.000 

000 

000.0 

1556* 

0.000 

000 

000.0 

1557* 

0.000 

000 

000.0 

1558* 

0.000 

000 

000.0 

1559* 

0.000 

000 

000.0 

1560* 

0.000 

000 

000.0 

1561  * 

0.000 

000 

000.0 

1562* 

0.000 

000 

000.0 

1563* 

0.000 

000 

000.0 

1564* 

0.000 

000 

000.0 

1565* 

0.000 

000 

000.0 

1566* 

0.000 

000 

000.0 

1567* 

0.000 

000 

000.0 

1568* 

0.000 

000 

000.0 

1569* 

0.000 

000 

000.0 

1570* 

0.000 

000 

000.0 

1571  * 

0.000 

000 

000.0 

1572* 

0.000 

000 

000.0 

1573* 

0.000 

000 

000.0 

1574* 

0.000 

000 

000.0 

1575* 

0.000 

000 

000.0 

1576* 

0.000 

000 

000.0 

1577* 

0.000 

000 

000.0 

1578* 

0.000 

000 

000.0 

1579* 

0.000 

000 

000.0 

1580* 

0.000 

000 

000.0 

1581  * 

0.000 

000 

000.0 

1582* 

0.000 

000 

000.0 

1583* 

0.000 

000 

000.0 

1584* 

0.000 

000 

000.0 

1585* 

0.000 

000 

000.0 

1586* 

0.000 

000 

000.0 

1587* 

0.000 

000 

000.0 

1588* 

0.000 

000 

000.0 

1589* 

0.000 

000 

000.0 

1590* 

0.000 

000 

000.0 

1591  * 

0.000 

000 

000.0 

1592* 

0.000 

000 

000.0 

1593* 

0.000 

000 

000.0 

1594* 

0.000 

000 

000.0 

1595* 

0.000 

000 

000.0 

1596* 

0.000 

000 

000.0 
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*  1597* 

*  1598* 

*  1599* 

*  1600* 
*  1601  * 
*  1602* 

*  1603* 

*  1604* 

*  1605  * 

*  1606* 

*  1607* 

*  1608* 

*  1609* 

*  1610* 
*  1611  * 
*  1612* 

*  1613* 

*  1614* 

*  1615* 

*  1616* 

*  1617* 

*  1618* 

*  1619* 

*  1620* 
*  1621  * 
*  1622* 

*  1623* 

*  1624* 

*  1625* 

*  1626* 

*  1627* 

*  1628* 

*  1629* 

*  1630* 

*  1631  * 

*  1632* 

*  1633* 

*  1634* 

*  1635* 

*  1636* 

*  1637* 

*  1638* 

*  1639* 

*  1640* 

*  1641  * 

*  1642* 

*  1643* 

*  1644* 

*  1645* 

*  1646* 

*  1647* 


0.000 

000 

000.0 

0.000 

000 

000.0 

0.000 

000 

000.0 

0.000 

000 

000.0 

1.250 

000 

000.0 

1.250 

015 

000.0 

1.250 

015 

090.0 

1.250 

015 

180.0 

1.250 

015 

270.0 

1.250 

030 

000.0 

1.250 

030 

045.0 

1.250 

030 

090.0 

1.250 

030 

135.0 

1.250 

030 

180.0 

1.250 

030 

225.0 

1.250 

030 

270.0 

1.250 

030 

315.0 

1.250 

045 

000.0 

1.250 

045 

030.0 

1.250 

045 

060.0 

1.250 

045 

090.0 

1.250 

045 

120.0 

1.250 

045 

150.0 

1.250 

045 

180.0 

1.250 

045 

210.0 

1.250 

045 

240.0 

1.250 

045 

270.0 

1.250 

045 

300.0 

1.250 

045 

330.0 

1.250 

060 

000.0 

1.250 

060 

022.5 

1.250 

060 

045.0 

1.250 

060 

067.5 

1.250 

060 

090.0 

1.250 

060 

112.5 

1.250 

060 

135.0 

1.250 

060 

157.5 

1.250 

060 

180.0 

1.250 

060 

202.5 

1.250 

060 

225.0 

1.250 

060 

247.5 

1.250 

060 

270.0 

1.250 

060 

292.5 

1.250 

060 

315.0 

1.250 

060 

337.* 

1.250 

075 

000.0 

1.250 

075 

018.0 

1.250 

075 

036.0 

1.250 

075 

054.0 

1.250 

075 

072.0 

1.250 

075 

090.0 
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1648* 

1.250 

075 

108.0 

1649* 

1.250 

075 

126.0 

1650* 

1.250 

075 

144.0 

1651  * 

1.250 

075 

162.0 

1652* 

1.250 

075 

180.0 

1653* 

1.250 

075 

198.0 

1654* 

1.250 

075 

216.0 

1655* 

1.250 

075 

234.0 

1656* 

1.250 

075 

252.0 

1657* 

1.250 

075 

270.0 

1658* 

1.250 

075 

288.0 

1659* 

1.250 

075 

306.0 

1660* 

1.250 

075 

324.0 

1661  * 

1.250 

075 

342.0 

1662* 

1.250 

090 

000.0 

1663* 

1.250 

090 

015.0 

1664* 

1.250 

090 

030.0 

1665* 

1.250 

090 

045.0 

1666* 

1.250 

090 

060.0 

1667* 

1.250 

090 

075.0 

1668* 

1.250 

090 

090.0 

1669* 

1.250 

090 

105.0 

1670* 

1.250 

090 

120.0 

1671  * 

1.250 

090 

135.0 

1672* 

1.250 

090 

150.0 

1673* 

1.250 

090 

165.0 

1674* 

1.250 

090 

180.0 

1675* 

1.250 

090 

195.0 

1676* 

1.250 

090 

210.0 

1677* 

1.250 

090 

225.0 

1678* 

1.250 

090 

240.0 

1679* 

1.250 

090 

255.0 

1680* 

1.250 

090 

270.0 

1681  * 

1.250 

090 

285.0 

1682* 

1.250 

090 

300.0 

1683* 

1.250 

090 

315.0 

1684* 

1.250 

090 

330.0 

1685* 

1.250 

090 

345.0 

1686* 

1.250 

105 

000.0 

1687* 

1.250 

105 

018.0 

1688* 

1.250 

105 

036.0 

1689* 

1.250 

105 

054.0 

1690* 

1.250 

105 

072.0 

1691  * 

1.250 

105 

090.0 

1692* 

1.250 

105 

108.0 

1693* 

1.250 

105 

126.0 

1694* 

1.250 

105 

144.0 

1695* 

1.250 

105 

162.0 

1696* 

1.250 

105 

180.0 

1697* 

1.250 

105 

198.0 

1698* 

1.250 

105 

216.0 
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*  1699*  1.250  105  234.0 

*  1700*  1.250  105  252.0 

*  1701  *  1.250  105  270.0 

*  1702*  1.250  105  288.0 

*  1703*  1.250  105  306.0 

*  1704*  1.250  105  324.0 

*  1705*  1.250  105  342.0 

*  1706*  1.250  120  000.0 

*  1707*  1.250  120  022.5 

*  1708*  1.250  120  045.0 

*  1709*  1.250  120  067.5 

*  1710*  1.250  120  090.0 

*  1711  *  1.250  120  112.5 

*  1712*  1.250  120  135.0 

*  1713*  1.250  120  157.5 

*  1714*  1.250  120  180.0 

*  1715*  1.250  120  202.5 

*  1716*  1.250  120  225.0 

*  1717*  1.250  120  247.5 

*  1718*  1.250  120  270.0 

*  1719*  1.250  120  292.5 

*  1720*  1.250  120  315.0 

*  1721  *  1.250  120  337.5 

*  1722*  1.250  135  000.0 

*  1723*  1.250  135  030.0 

*  1724*  1.250  135  060.0 

*  1725*  1.250  135  090.0 

*  1726*  1.250  135  120.0 

*  1727*  1.250  135  150.0 

*  1728*  1.250  135  180.0 

*  1729*  1.250  135  210.0 

*  1730*  1.250  135  240.0 

*  1731  *  1.250  135  270.0 

*  1732*  1.250  135  300.0 

*  1733*  1.250  135  330.0 

*  1734*  1.250  150  000.0 

*  1735*  1.250  150  045.0 

*  1736*  1.250  150  090.0 

*  1737*  1.250  150  135.0 

*  1738*  1.250  150  180.0 

*  1739*  1.250  150  225.0 

*  1740*  1.250  150  270.0 

*  1741  *  1.250  150  315.0 

*  1742*  1.250  165  000.0 

*  1743*  1.250  165  090.0 

*  1744*  1.250  165  180.0 

*  1745*  1.250  165  270.0 

*  1746*  1.250  180  000.0 

*  1747  *  0.000  000  000.0 

*  1748  *  0.000  000  .  30.0 

*  1749  *  0.000  000  000.0 
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1750*  0.000 

000 

000.0 

1751  *  0.000 

000 

000.0 

1752  *  0.000 

000 

000.0 

1753  *  0.000 

000 

000.0 

1754  *  0.000 

000 

000.0 

1755  *  0.000 

000 

000.0 

1756  *  0.000 

000 

000.0 

1757  *  0.000 

000 

000.0 

1758  *  0.000 

000 

000.0 

1759  *  0.000 

000 

000.0 

1760  *  0.000 

000 

000.0 

1761  *  0.000 

000 

000.0 

1762  *  0.000 

000 

000.0 

1763  *  0.000 

000 

000.0 

1764  *  0.000 

000 

000.0 

1765  *  0.000 

000 

000.0 

1766  *  0.000 

000 

000.0 

1767  *  0.000 

000 

000.0 

1768  *  0.000 

000 

000.0 

1769  *  0.000 

000 

000.0 

1770  *  0.000 

000 

000.0 

1771  *  0.000 

000 

000.0 

1772  *  0.000 

000 

000.0 

1773  *  0.000 

000 

000.0 

1774  *  0.000 

000 

000.0 

1775  *  0.000 

000 

000.0 

1776  *  0.000 

000 

000.0 

1777  *  0.000 

000 

000.0 

1778*  0.000 

000 

000.0 

1779*  C  000 

000 

000.0 

1780  *  0.000 

000 

000.0 

1781  *  0.000 

000 

000.0 

1782  *  0.000 

000 

000.0 

1783  *  0.000 

000 

000.0 

1784  *  0.000 

000 

000.0 

1785*  0.000 

000 

000.0 

1786  *  0.000 

000 

000.0 

1787  *  0.000 

000 

000.0 

1788  *  0.000 

000 

000.0 

1789  *  0.000 

000 

000.0 

1790  *  0.000 

000 

000.0 

1791  *  0.000 

000 

000.0 

1792  *  0.000 

000 

000.0 

1793  *  0.000 

000 

000.0 

1794  *  0.000 

000 

000.0 

1795  *  0.000 

000 

000.0 

1796*  0.000 

000 

000.0 

1797  *  0.000 

000 

000.0 

1798  *  0.000 

000 

000.0 

1799  *  0.000 

000 

000.0 

1800  *  0.000 

000 

000.0 
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1801  * 

1.500 

000 

000.0 

1802* 

1.500 

015 

000.0 

1803* 

1.500 

015 

090.0 

1804* 

1.500 

015 

180.0 

1805* 

1.500 

015 

270.0 

1806* 

1.500 

030 

000.0 

1807* 

1.500 

030 

045.0 

1808* 

1.500 

030 

090.0 

1809* 

1.500 

030 

135.0 

1810* 

1.500 

030 

180.0 

1811  * 

1.500 

030 

225.0 

1812* 

1.500 

030 

270.0 

1813* 

1.500 

030 

315.0 

1814* 

1.500 

045 

000.0 

1815* 

1.500 

045 

030.0 

1816* 

1.500 

045 

060.0 

1817* 

1.500 

045 

090.0 

1818* 

1.500 

045 

120.0 

1819* 

1.500 

045 

150.0 

1820* 

1.500 

045 

180.0 

1821  * 

1.500 

045 

210.0 

1822* 

1.500 

045 

240.0 

1823* 

1.500 

045 

270.0 

1824* 

1.500 

045 

300.0 

1825* 

1.500 

045 

330.0 

1826* 

1.500 

060 

000.0 

1827* 

1.500 

060 

022.5 

1823* 

1.500 

060 

045.0 

1829* 

1.500 

060 

067.5 

1830* 

1.500 

060 

090.0 

1831  * 

1.500 

060 

112.5 

1832* 

1.500 

060 

135.0 

1833* 

1.500 

060 

157.5 

1834* 

1.500 

060 

180.0 

1835* 

1.500 

060 

202.5 

1836* 

1.500 

060 

225.0 

1837* 

1.500 

060 

247.5 

1838* 

1.500 

060 

270.0 

1839* 

1.500 

060 

292.5 

1840* 

1.500 

060 

315.0 

1841  * 

1.500 

060 

337.5 

1842* 

1.500 

075 

000.0 

1843* 

1.500 

075 

018.0 

1844* 

1.500 

075 

036.0 

1845* 

1.500 

075 

054.0 

1846* 

1.500 

075 

072.0 

1847* 

1.500 

075 

090.0 

1848* 

1.500 

075 

108.0 

1849* 

1.500 

075 

126.0 

1850* 

1.500 

075 

144.0 

1851  * 

1.500 

075 

162.0 
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*  1903  *  1.500  105  306.0 

*  1904*  1.500  105  324.0 

*  1905*  1.500  105  342.0 

*  1906*  1.500  120  000.0 

*  1907*  1.500  120  022.5 

*  1908*  1.500  120  045.0 

*  1909*  1.500  120  067.5 

*  1910*  1.500  120  090.0 

*  1911  *  1.500  120  112.5 

*  1912*  1.500  120  135.0 

*  1913*  1.500  120  157.5 

*  1914*  1.500  120  180.0 

*  1915*  1.500  120  202.5 

*  1916*  1.500  120  225.0 

*  1917*  1.500  120  247.5 

*  1918*  1.500  120  270.0 

*  1919*  1.500  120  292.5 

*  1920*  1.500  120  315.0 

*  1921  *  1.500  120  337.5 

*  1922*  1.500  135  000.0 

*  1923*  1.500  135  030.0 

*  1924*  1.500  135  060.0 

*  1925*  1.500  135  090.0 

*  1926*  1.500  135  120.0 

*  1927*  1.500  135  150.0 

*  1928*  1.500  135  180.0 

*  1929*  1.500  135  210.0 

*  1930*  1.500  135  240.0 

*  1931  *  1.500  135  270.0 

*  1932*  1.500  135  300.0 

*  1933*  1.500  135  330.0 

*  1934*  1.500  150  000.0 

*  1935*  1.500  150  045.0 

*  1936*  1.500  150  090.0 

*  1937*  1.500  150  135.0 

*  1938*  1.500  150  180.0 

*  1939*  1.500  150  225.0 

*  1940*  1.500  150  270.0 

*  1941  *  1.500  150  315.0 

*  1942*  1.500  165  000.0 

*  1943*  1.500  165  090.0 

*  1944*  1.500  165  180.0 

*  1945*  1.500  165  270.0 

*  1946*  1.500  180  000.0 

*  1947  *  0.000  000  000.0 

*  1948  *  0.000  000  000.0 

*  1949  *  0.000  000  000.0 

*  1950  *  0.000  000  000.0 

*  1951  *  0.000  000  000.0 

*  1952  *  0.000  000  000.0 

*  1953  *  0.000  000  000.0 
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1954  *  0.000 

000 

000.0 

1955  *  0.000 

000 

000.0 

1956  *  0.000 

000 

000.0 

1957  *  0.000 

000 

000.0 

1958  *  0.000 

000 

000.0 

1959  *  0.000 

000 

000.0 

1960  *  0.000 

000 

000.0 

1961  *  0.000 

000 

000.0 

1962  *  0.000 

000 

000.0 

1963  *  0.000 

000 

000.0 

1964  *  0.000 

000 

000.0 

1965  *  0.000 

000 

000.0 

1966  *  0.000 

000 

000.0 

1967  *  0.000 

000 

000.0 

1968  *  0.000 

000 

000.0 

1969  *  0.000 

000 

000.0 

1970  *  0.000 

000 

000.0 

1971  *  0.000 

000 

000.0 

1972  *  0.000 

000 

000.0 

1973  *  0.000 

000 

000.0 

1974  *  0.000 

000 

000.0 

1975  *  0.000 

000 

000.0 

1976  *  0.000 

000 

000.0 

1977  *  0.000 

000 

000.0 

1978  *  0.000 

000 

000.0 

1979  *  0.000 

000 

000.0 

1980  *  0.000 

000 

000.0 

1981  *  0.000 

000 

000.0 

1982  *  0.000 

000 

000.0 

1983  *  0.000 

000 

000.0 

1984  *  0.000 

000 

000.0 

1985  *  0.000 

000 

000.0 

1986  *  0.000 

000 

000.0 

1987  *  0.000 

000 

000.0 

1988  *  0.000 

000 

000.0 

1989  *  0.000 

000 

000.0 

1990  *  0.000 

000 

000.0 

1991  *  0.000 

000 

000.0 

1992  *  0.000 

000 

000.0 

1993  *  0.000 

000 

000.0 

1994  *  0.000 

000 

000.0 

1995  *  0.000 

000 

000.0 

1996  *  0.000 

000 

000.0 

1997  *  0.000 

000 

000.0 

1998  *  0.000 

000 

000.0 

1999  *  0.000 

000 

000.0 

2000  *  0.000 

000 

000.0 

2001  *  1.750 

000 

000.0 

2002*  1.750 

015 

000.0 

2003*  1.750 

015 

090.0 

2004*  1.750 

015 

180.0 

*  2005*  1.750  015  270.0 

*  2006’  1.750  030  000.0 

*  2007*  1.750  030  045.0 

*  2008*  1.750  030  090.0 

*  2009*  1.750  030  135.0 

*  2010*  1.750  030  180.0 

*  2011  *  1.750  030  225.0 

*  2012*  1.750  030  270.0 

*  2013*  1.750  030  315.0 

*  2014*  1.750  045  000.0 

*  2015*  1.750  045  030.0 

*  2016*  1.750  045  060.0 

*  2017*  1.750  045  090.0 

*  2018*  1.750  045  120.0 

*  2019*  1.750  045  150.0 

*  2020*  1.750  045  180.0 

*  2021  *  1.750  045  210.0 

*  2022*  1.750  045  240.0 

*  2023*  1.750  045  270.0 

*  2024*  1.750  045  300.0 

*  2025*  1.750  045  330.0 

*  2026*  1.750  060  000.0 

*  2027*  1.750  060  022.5 

*  2028*  1.750  060  045.0 

*  2029*  1.750  060  067.5 

*  2030*  1.750  060  090.0 

*  2031  *  1.750  060  112.5 

*  2032*  1.750  060  135.0 

*  2033*  1.750  060  157.5 

*  2034*  1.750  060  180.0 

*  2035*  1.750  060  202.5 

*  2036*  1.750  060  225.0 

*  2037*  1.750  060  247.5 

*  2038*  1.750  060  270.0 

*  2039*  1.750  060  292.5 

*  2040*  1.750  060  315.0 

*  2041  *  1.750  060  337.5 

*  2042*  1.750  075  000.0 

*  2043*  1.750  075  018.0 

*  2044*  1.750  075  036.0 

*  2045*  1.750  075  054.0 

*  2046*  1.750  075  072.0 

*  2047*  1.750  075  090.0 

*  2048*  1.750  075  108.0 

*  2049*  1.750  075  126.0 

*  2050*  1.750  075  144.0 

*  2051  *  1.750  075  162.0 

*  2052*  1.750  075  180.0 

*  2053*  1.750  075  198.0 

*  2054*  1.750  075  216.0 

*  2055*  1.750  075  234.0 
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*  2056*  1.750  075  252.0 

*  2057*  1.750  075  270.0 

*  2058*  1.750  075  288.0 

*  2059*  1.750  075  306.0 

*  2060*  1.750  075  324.0 

*  2061  *  1.750  075  342.0 

*  2062*  1.750  090  000.0 

*  2063*  1.750  090  015.0 

*  2064*  1.750  090  030.0 

*  2065*  1.750  090  045.0 

*  2066*  1.750  090  060.0 

*  2067*  1.750  090  075.0 

*  2068*  1.750  090  090.0 

*  2069*  1.750  090  105.0 

*  2070*  1.750  090  120.0 

*  2071  *  1.750  090  135.0 

*  2072*  1.750  090  150.0 

*  2073*  1.750  090  165.0 

*  2074*  1.750  090  180.0 

*  2075*  1.750  090  195.0 

*  2076*  1.750  090  210.0 

*  2077*  1.750  090  225.0 

*  2078*  1.750  090  240.0 

*  2079*  1.750  090  255.0 

*  2080*  1.750  090  270.0 

*  2081  *  1.750  090  285.0 

*  2082*  1.750  090  300.0 

*  2083*  1.750  090  315.0 

*  2084*  1.750  090  330.0 

*  2085*  1.750  090  345.0 

*  2086*  1.750  105  000.0 

*  2087*  1.750  105  018.0 

*  2088*  1.750  105  036.0 

*  2089*  1.750  105  054.0 

*  2090*  1.750  105  072.0 

*  2091  *  1.750  105  090.0 

*  2092*  1.750  105  108.0 

*  2093*  1.750  105  126.0 

*  2094*  1.750  105  144.0 

*  2095*  1.750  105  162.0 

*  2096*  1.750  105  180.0 

*  2097*  1.750  105  198.0 

*  2098*  1.750  105  216.0 

*  2099*  1.750  105  234.0 

*  2100*  1.750  105  252.0 

*  2101  *  1.750  105  270.0 

*  2102*  1.750  105  288.0 

*  2103*  1.750  105  306.0 

*  2104*  1.750  105  324.0 

*  2105*  1.750  105  342.0 

*  2106*  1.750  120  000.0 
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2107* 

1.750 

120 

022.5 

2108* 

1.750 

120 

045.0 

2109* 

1.750 

120 

067.5 

2110* 

1.750 

120 

090.0 

2111  * 

1.750 

120 

112.5 

2112* 

1.750 

120 

135.0 

2113* 

1.750 

120 

157.5 

2114* 

1.750 

120 

180.0 

2115* 

1.750 

120 

202.5 

2116* 

1.750 

120 

225.0 

2117* 

1.750 

120 

247.5 

2118* 

1.750 

120 

270.0 

2119* 

1.750 

120 

292.5 

2120* 

1.750 

120 

315.0 

2121  * 

1.750 

120 

337.5 

2122* 

1.750 

135 

000.0 

2123* 

1.750 

135 

030.0 

2124* 

1.750 

135 

060.0 

2125* 

1.750 

135 

090.0 

2126* 

1.750 

135 

120.0 

2127* 

1.750 

135 

150.0 

2128* 

1.750 

135 

180.0 

2129* 

1.750 

135 

210.0 

2130* 

1.750 

135 

240.0 

2131  * 

1.750 

135 

270.0 

2132* 

1.750 

135 

300.0 

2133* 

1.750 

135 

330.0 

2134* 

1.750 

150 

000.0 

2135* 

1.750 

150 

045.0 

2136* 

1.750 

150 

090.0 

2137* 

1.750 

150 

135.0 

2138* 

1.750 

150 

180.0 

2139* 

1.750 

150 

225.0 

2140* 

1.750 

150 

270.0 

2141  * 

1.750 

150 

315.0 

2142* 

1.750 

165 

000.0 

2143* 

1.750 

165 

090.0 

2144* 

1.750 

165 

180.0 

2145* 

1.750 

165 

270.0 

2146* 

1.750 

180 

000.0 

2147* 

0.000 

000 

000.0 

2148* 

0.000 

000 

000.0 

2149* 

0.000 

000 

000.0 

2150* 

0.000 

000 

000.0 

2151  * 

0.000 

000 

000.0 

2152* 

0.000 

000 

000.0 

2153* 

0.000 

000 

000.0 

2154* 

0.000 

000 

000.0 

2155* 

0.000 

000 

000.0 

2156* 

0.000 

000 

000.0 

2157* 

0.000 

000 

000.0 
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2158  *  0.000 

000 

000.0 

2159  *  0.000 

000 

000.0 

2160  *  0.000 

000 

000.0 

2161  *  0.000 

000 

000.0 

2162  *  0.000 

000 

000.0 

2163  *  0.000 

000 

000.0 

2164  *  0.000 

000 

000.0 

2165*  0.000 

000 

000.0 

2166  *  0.000 

000 

000.0 

2167  *  0.000 

000 

000.0 

2168  *  0.000 

000 

000.0 

2169  *  0.000 

000 

000.0 

2170  *  0.000 

000 

000.0 

2171*  0.000 

000 

000.0 

2172  *  0.000 

000 

000.0 

2173  *  0.000 

000 

000.0 

2174  *  0.000 

000 

000.0 

2175  *  0.000 

000 

000.0 

2176*  0.000 

000 

000.0 

2177  *  0.000 

000 

000.0 

2178*  0.000 

000 

000.0 

2179  *  0.000 

000 

000.0 

2180  *  0.000 

000 

000.0 

2181  *  0.000 

000 

000.0 

2182  *  0.000 

000 

000.0 

2183*  0.000 

000 

000.0 

2184  *  0.000 

000 

000.0 

2185*  0.000 

000 

000.0 

2186  *  0.000 

000 

000.0 

2187  *  0.000 

000 

000.0 

2188  *  0.000 

000 

000.0 

2189  *  0.000 

000 

000.0 

2190  *  0.000 

000 

000.0 

2191  *  0.000 

000 

000.0 

2192  *  0.000 

000 

000.0 

2193  *  0.000 

000 

000.0 

2194  *  0.000 

000 

000.0 

2195  *  0.000 

000 

000.0 

2196  *  0.000 

000 

000.0 

2197*  0.000 

000 

000.0 

2198  *  0.000 

000 

000.0 

2199  *  0.000 

000 

000.0 

2200  *  0.000 

000 

000.0 

2201  *  2.000 

000 

000.0 

2202  *  2.000 

015 

000.0 

2203  *  2.000 

015 

090.0 

2204  *  2.000 

015 

180.0 

2205  *  2.000 

015 

270.0 

2206  *  2.000 

030 

000.0 

2207  *  2.000 

030 

045.0 

2208  *  2.000 

030 

090.0 
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2209* 

2.000 

030 

135.0 

2210* 

2.000 

030 

180.0 

2211  * 

2.000 

030 

225.0 

2212* 

2.000 

030 

270.0 

2213* 

2.000 

030 

315.0 

2214* 

2.000 

045 

000.0 

2215* 

2.000 

045 

030.0 

2216* 

2.000 

045 

060.0 

2217* 

2.000 

045 

090.0 

2218* 

2.000 

045 

120.0 

2219* 

2.000 

045 

150.0 

2220* 

2.000 

045 

180.0 

2221  * 

2.000 

045 

210.0 

0000  * 

2.000 

045 

240.0 

2223* 

2.000 

045 

270.0 

2224* 

2.000 

045 

300.0 

2225* 

2.000 

045 

330.0 

2226* 

2.000 

060 

000.0 

2227* 

2.000 

060 

022.5 

2228* 

2.000 

060 

045.0 

2229* 

2.000 

060 

067.5 

2230* 

2.000 

060 

090.0 

2231  * 

2.000 

060 

112.5 

2232* 

2.000 

060 

135.0 

2233* 

2.000 

060 

157.5 

2234* 

2.000 

060 

180.0 

2235* 

2.000 

060 

202.5 

2236* 

2.000 

060 

225.0 

2237* 

2.000 

060 

247.5 

2238* 

2.000 

060 

270.0 

2239* 

2.000 

060 

292.5 

2240* 

2.000 

060 

315.0 

2241  * 

2.000 

060 

337.5 

2242* 

2.000 

075 

000.0 

2243* 

2.000 

075 

018.0 

2244* 

2.000 

075 

036.0 

2245* 

2.000 

075 

054.0 

2246* 

2.000 

075 

072.0 

2247* 

2.000 

075 

090.0 

2248* 

2.000 

075 

108.0 

2249* 

2.000 

075 

126.0 

2250* 

2.000 

075 

144.0 

2251  * 

2.000 

075 

162.0 

2252* 

2.000 

075 

180.0 

2253* 

2.000 

075 

198.0 

2254* 

2.000 

075 

216.0 

2255* 

2.000 

075 

234.0 

2256* 

2.000 

075 

252.0 

2257* 

2.000 

075 

270.0 

2258* 

2.000 

075 

288.0 

2259* 

2.000 

075 

306.0 
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*  2260  *  2.000 

075 

324.0 

*  2261  *  2.000 

075 

342.0 

*  2262  *  2.000 

090 

000.0 

*  2263  *  2.000 

090 

015.0 

*  2264  *  2.000 

090 

030.0 

*  2265  *  2.000 

090 

045.0 

*  2266  *  2.000 

090 

060.0 

*  2267  *  2.000 

090 

075.0 

*  2268  *  2.000 

090 

090.0 

*  2269  *  2.000 

090 

105.0 

*  2270  *  2.000 

090 

120.0 

*  2271  *  2.000 

090 

135.0 

*  2272*  2.000 

090 

150.0 

*  2273  *  2.000 

090 

165.0 

*  2274  *  2.000 

090 

180.0 

*  2275  *  2.000 

090 

195.0 

*  2276  *  2.000 

090 

210.0 

*  2277  *  2.000 

090 

225.0 

*  2278  *  2.000 

090 

240.0 

*  2279  *  2.000 

090 

255.0 

*  2280  *  2.000 

090 

270.0 

*  2281  *  2.000 

090 

285.0 

*  2282  *  2.000 

090 

300.0 

*  2283  *  2.000 

090 

315.0 

*  2284  *  2.000 

090 

330.0 

*  2285  *  2.000 

090 

345.0 

*  2286  *  2.000 

105 

000.0 

*  2287  *  2.000 

105 

018.0 

*  2288  *  2.000 

105 

036.0 

*  2289  *  2.000 

105 

054.0 

*  2290  *  2.000 

105 

072.0 

*  2291  *  2.000 

105 

090.0 

*  2292  *  2.000 

105 

108.0 

*  2293  *  2.000 

105 

126.0 

*  2294  *  2.000 

105 

144.0 

*  2295  *  2.000 

105 

162.0 

*  2296  *  2.000 

105 

180.0 

*  2297  *  2.000 

105 

198.0 

*  2298  *  2.000 

105 

216.0 

*  2299  *  2.000 

105 

234.0 

*  2300  *  2.000 

105 

252.0 

*  2301  *  2.000  105  270.0 

*  2302  *  2.000  105  288.0 

*  2303  *  2.000  105  306.0 

*  2304  *  2.000  105  324.0 

*  2305  *  2.000  105  342.0 

*  2306  *  2.000  120  000.0 

*  2307  *  2.000  120  022.5 

*  2308  *  2.000  120  045.0 

*  2309  *  2.000  120  067.5 

*  2310  *  2.000  120  090.0 
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2311  * 

2.000 

120 

112.5 

2312* 

2.000 

120 

135.0 

2313* 

2.000 

120 

157.5 

2314* 

2.000 

120 

180.0 

2315* 

2.000 

120 

202.5 

2316* 

2.000 

120 

225.0 

2317* 

2.000 

120 

247.5 

2318* 

2.000 

120 

270.0 

2319* 

2.000 

120 

292.5 

2320* 

2.000 

120 

315.0 

2321  * 

2.000 

120 

337.5 

2322* 

2.000 

135 

000.0 

2323* 

2.000 

135 

030.0 

2324* 

2.000 

135 

060.0 

2325* 

2.000 

135 

090.0 

2326* 

2.000 

135 

120.0 

2327* 

2.000 

135 

150.0 

2328* 

2.000 

135 

180.0 

2329* 

2.000 

135 

210.0 

2330* 

2.000 

135 

240.0 

2331  * 

2.000 

135 

270.0 

2332* 

2.000 

135 

300.0 

2333* 

2.000 

135 

330.0 

2334* 

2.000 

150 

000.0 

2335* 

2.000 

150 

045.0 

2336* 

2.000 

150 

090.0 

2337* 

2.000 

150 

135.0 

2338* 

2.000 

150 

180.0 

2339* 

2.000 

150 

225.0 

2340* 

2.000 

150 

270.0 

2341  * 

2.000 

150 

315.0 

2342* 

2.000 

165 

000.0 

2343* 

2.000 

165 

090.0 

2344* 

2.000 

165 

180.0 

2345* 

2.000 

165 

270.0 

2346* 

2.000 

180 

000.0 

2347* 

0.000 

000 

000.0 

2348* 

0.000 

000 

000.0 

2349* 

0.000 

000 

000.0 

2350* 

0.000 

000 

000.0 

2351  * 

0.000 

000 

000.0 

2352* 

0.000 

000 

000.0 

2353* 

0.000 

000 

000.0 

2354* 

0.000 

000 

000.0 

2355* 

0.000 

000 

000.0 

2356* 

0.000 

000 

000.0 

2357* 

0.000 

000 

000.0 

2358* 

0.000 

000 

000.0 

# 

2359* 

0.000 

000 

000.0 

# 

2360* 

0.000 

000 

000.0 

* 

2361  * 

0.000 

000 

000.0 
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*  2362  *  0.000  000  000.0 

*  2363  *  0.000  000  000.0 

*  2364  *  0.000  000  000.0 

*  2365  *  0.000  000  000.0 

*  2366  *  0.000  000  000.0 

*  2367  *  0.000  000  000.0 

*  2368  *  0.000  000  000.0 

*  2369  *  0.000  000  000.0 

*  2370  *  0.000  000  000.0 

*  2371  *  0.000  000  000.0 

*  2372  *  0.000  000  000.0 

*  2373  *  0.000  000  000.0 

*  2374  *  0.000  000  000.0 

*  2375  *  0.000  000  000.0 

*  2376  *  0.000  000  000.0 

*  2377  *  0.000  000  000.0 

*  2378  *  0.000  000  000.0 

*  2379  *  0.000  000  000.0 

*  2380  *  0.000  000  000.0 

*  2381  *  0.000  000  000.0 

*  2382  *  0.000  000  000.0 

*  2383  *  0.000  000  000.0 

*  2384  *  0.000  000  000.0 

*  2385  *  0.000  000  000.0 

*  2386  *  0.000  000  000.0 

*  2387  *  0.000  000  000.0 

*  2388  *  0.000  000  000.0 

*  2389  *  0.000  000  000.0 

*  2390  *  0.000  000  000.0 

*  2391  *  0.000  000  000.0 

*  2392  *  0.000  000  000.0 

*  2393  *  0.000  000  000.0 

*  2394  *  0.000  000  000.0 

*  2395  *  0.000  000  000.0 

*  2396  *  0.000  000  000.0 

*  2397  *  0.000  000  000.0 

*  2398  *  0.000  000  000.0 

*  2399  *  0.000  000  000.0 

*  2400  *  0.000  000  000.0 

*  2401  *  2.250  000  000.0 

*  2402  *  2.250  015  000.0 

*  2403  *  2.250  015  090.0 

*  2404  *  2.250  015  180.0 

*  2405*  2.250  015  270.0 

*  2406  *  2.250  030  000.0 

*  2407  *  2.250  030  045.0 

*  2408  *  2.250  030  090.0 

*  2409  *  2.250  030  135.0 

*  2410  *  2.250  030  180.0 

*  2411  *  2.250  030  225.0 

*  2412  *  2.250  030  270.0 
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2413* 

2.250 

030 

315.0 

2414* 

2.250 

045 

000.0 

2415* 

2.250 

045 

030.0 

2416* 

2.250 

045 

060.0 

2417* 

2.250 

045 

090.0 

2418* 

2.250 

045 

120.0 

2419* 

2.250 

045 

150.0 

2420* 

2.250 

045 

180.0 

2421  * 

2.250 

045 

210.0 

2422* 

2.250 

045 

240.0 

2423* 

2.250 

045 

270.0 

2424* 

2.250 

045 

300.0 

2425* 

2.250 

045 

330.0 

2426* 

2.250 

060 

000.0 

2427* 

2.250 

060 

022.5 

2428* 

2.250 

060 

045.0 

2429* 

2.250 

060 

067.5 

2430* 

2.250 

060 

090.0 

2431  * 

2.250 

060 

112.5 

2432* 

2.250 

060 

135.0 

2433* 

2.250 

060 

157.5 

2434* 

2.250 

060 

180.0 

2435* 

2.250 

060 

202.5 

2436* 

2.250 

060 

225.0 

2437* 

2.250 

060 

247.5 

2438* 

2.250 

060 

270.0 

2439* 

2.250 

060 

292.5 

2440* 

2.250 

060 

315.0 

2441  * 

2.250 

060 

337.5 

2442* 

2.250 

075 

000.0 

2443* 

2.250 

075 

018.0 

2444* 

2.250 

075 

036.0 

2445* 

2.250 

075 

054.0 

2446* 

2.250 

075 

072.0 

2447* 

2.250 

075 

090.0 

2448* 

2.250 

075 

108.0 

2449* 

2.250 

075 

126.0 

2450* 

2.250 

075 

144.0 

2451  * 

2.250 

075 

162.0 

2452* 

2.250 

075 

180.0 

2453* 

2.250 

075 

198.0 

2454* 

2.250 

075 

216.0 

2455* 

2.250 

075 

234.0 

2456* 

2.250 

075 

252.0 

2457* 

2.250 

075 

270.0 

2458* 

2.250 

075 

288.0 

2459* 

2.250 

075 

306.0 

2460* 

2.250 

075 

324.0 

2461  * 

2.250 

075 

342.0 

2462* 

2.250 

090 

000.0 

2463* 

2.250 

090 

015.0 
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*  2464  *  2250  090  030.0 

*  2465  *  2.250  090  045.0 

*  2466  *  2.250  090  060.0 

*  2467  *  2.250  090  075.0 

*  2466  *  2.250  090  090.0 

*  2469  *  2.250  090  105.0 

*  2470  *  2.250  090  120.0 

*  2471  *  2.250  090  135.0 

*  2472  *  2.250  090  150.0 

*  2473  *  2.250  090  165.0 

*  2474  *  2.250  090  180.0 

*  2475  *  2.250  090  195.0 

*  2476  *  2.250  090  210.0 

*  2477  *  2.250  090  225.0 

*  2478  *  2.250  090  240.0 

*  2479  *  2.250  090  255.0 

*  2480  *  2.250  090  270.0 

*  2481  *  2.250  090  285.0 

*  2462  *  2.250  090  300.0 

*  2483  *  2.250  090  315.0 

*  2484  *  2.250  090  330.0 

*  2485*  2.250  090  345.0 

*  2486  *  2.250  105  000.0 

*  2487  *  2.250  105  018.0 

*  2488  *  2.250  105  036.0 

*  2489  *  2.250  105  054.0 

*  2490  *  2.250  105  072.0 

*  2491  *  2.250  105  090.0 

*  2492  *  2.250  105  108.0 

*  2493  *  2.250  105  126.0 

*  2494  *  2.250  105  144.0 

*  2495  *  2.250  105  162.0 

*  2496  *  2.250  105  180.0 

*  2497  *  2.250  105  198.0 

*  2498  *  2.250  105  216.0 

*  2499  *  2.250  105  234.0 

*  2500  *  2.250  105  252.0 

*  2501  *  2.250  105  270.0 

*  2502  *  2.250  105  288.0 

*  2503  *  2.250  105  306.0 

*  2504  *  2.250  105  324.0 

*  2505  *  2.250  105  342.0 

*  2506  *  2.250  120  000.0 

*  2507  *  2.250  120  022.5 

*  2508  *  2.250  120  045.0 

*  2509  *  2.250  120  067.5 

*  2510  *  2.250  120  090.0 

*  2511  *  2.250  120  112.5 

*  2512  *  2.250  120  135.0 

*  2513  *  2.250  120  157.5 

*  2514  *  2.250  120  180.0 
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*  2515  *  2.250  120  202.5 

*  2510  *  2.250  120  225.0 

*  2517  *  2.250  120  247.5 

*  2518  *  2.250  120  270.0 

*  2519  *  2.250  120  292.5 

*  2520  *  2.250  120  315.0 

*  2521  *  2.250  120  337.5 

*  2522  *  2.250  135  000.0 

*  2523  *  2.250  135  030.0 

*  2524  *  2.250  135  060.0 

*  2525  *  2.250  135  090.0 

*  2526  *  2.250  135  120.0 

*  2527  *  2.250  135  150.0 

*  2528  *  2.250  135  180.0 

*  2529  *  2.250  135  210.0 

*  2530  *  2.250  135  240.0 

*  2531  *  2.250  135  270.0 

*  2532  *  2.250  135  300.0 

*  2533  *  2.250  135  330.0 

*  2534  *  2.250  150  000.0 

*  2535  *  2.250  150  045.0 

*  2536  *  2.250  150  090.0 

*  2537  *  2.250  150  135.0 

*  2538  *  2.250  150  180.0 

*  2539  *  2.250  150  225.0 

*  2540  *  2.250  150  270.0 

*  2541  *  2.250  150  315.0 

*  2542  *  2.250  165  000.0 

*  2543  *  2.250  165  090.0 

*  2544  *  2.250  165  180.0 

*  2545  *  2.250  165  270.0 

*  2546  *  2.250  180  000.0 

*  2547  *  0.000  000  000.0 

*  2548  *  0.000  000  000.0 

*  2549  *  0.000  000  000.0 

*  2550  *  0.000  000  000.0 

*  2551  *  0.000  000  000.0 

*  2552  *  0.000  000  000.0 

*  2553  *  0.000  000  000.0 

*  2554  *  0.000  000  000.0 

*  2555  *  0.000  000  000.0 

*  2556  *  0.000  000  000.0 

*  2557  *  0.000  000  000.0 

*  2558  *  0.000  000  000.0 

*  2559  *  0.000  000  000.0 

*  2560  *  0.000  000  000.0 

*  2561  *  0.000  000  000.0 

*  2562  *  0.000  000  000.0 

*  2563  *  0.000  000  000.0 

*  2564  *  0.000  000  000.0 

*  2565  *  0.000  000  000.0 
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*  2566  *  0.000  000  000.0 

*  2567  *  0.000  000  000.0 

*  2568  *  0.000  000  000.0 

*  2569  *  0.000  000  000.0 

*  2570  *  0.000  000  000.0 

*  2571  *  0.000  000  000.0 

*  2572  *  0.000  000  000.0 

*  2573  *  0.000  000  000.0 

*  2574  *  0.000  000  000.0 

*  2575  *  0.000  000  000.0 

*  2576  *  0.000  000  000.0 

*  2577  *  0.000  000  000.0 

*  2578  *  0.000  000  000.0 

*  2579  *  0.000  000  000.0 

*  2580  *  0.000  000  000.0 

*  2581  *  0.000  000  000.0 

*  2582  *  0.000  000  000.0 

*  2583  *  0.000  000  000.0 

*  2584  *  0.000  000  000.0 

*  2585  *  0.000  000  000.0 

*  2586  *  0.000  000  000.0 

*  2587  *  0.000  000  000.0 

*  2588  *  0.000  000  000.0 

*  2589  *  0.000  000  000.0 

*  2590  *  0.000  000  000.0 

*  2591  *  0.000  000  000.0 

*  2592  *  0.000  000  000.0 

*  2593  *  0.000  000  000.0 

*  2594  *  0.000  000  000.0 

*  2595  *  0.000  000  000.0 

*  2596  *  0.000  000  000.0 

*  2597  *  0.000  000  000.0 

*  2598  *  0.000  000  000.0 

*  2599  *  0.000  000  000.0 

*  2600  *  0.000  000  000.0 

*  2601  *  2.500  000  000.0 

*  2602  *  2.500  015  000.0 

*  2603  *  2.500  015  090.0 

*  2604  *  2.500  015  180.0 

*  2605  *  2.500  015  270.0 

*  2606  *  2.500  030  000.0 

*  2607  *  2.500  030  045.0 

*  2608  *  2.500  030  090.0 

*  2609  *  2.500  030  135.0 

*  2610  *  2.500  030  180.0 

*  2611  *  2.500  030  225.0 

*  2612  *  2.500  030  270.0 

*  2613  *  2.500  030  315.0 

*  2614  *  2.500  045  000.0 

*  2615  *  2.500  045  030.0 

*  2616  *  2.500  045  060.0 
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2617* 

2.500 

045 

090.0 

2618* 

2.500 

045 

120.0 

2619* 

2.500 

045 

150.0 

2620* 

2.500 

045 

180.0 

2621  * 

2.500 

045 

210.0 

2622* 

2.500 

045 

240.0 

2623* 

2.500 

045 

270.0 

2624* 

2.500 

045 

300.0 

2625* 

2.500 

045 

330.0 

2626* 

2.500 

060 

000.0 

2627* 

2.500 

060 

022.5 

2628* 

2.500 

060 

045.0 

2629* 

2.500 

060 

067.5 

2630* 

2.500 

060 

090.0 

2631  * 

2.500 

060 

112.5 

2632* 

2.500 

060 

135.0 

2633* 

2.500 

060 

157.5 

2634* 

2.500 

060 

180.0 

2635* 

2.500 

060 

202.5 

2636* 

2.500 

060 

225.0 

2637* 

2.500 

060 

247.5 

2638* 

2.500 

060 

270.0 

2639* 

2.500 

060 

292.5 

2640* 

2.500 

060 

315.0 

2641  * 

2.500 

060 

337.5 

2642* 

2.500 

075 

000.0 

2643* 

2.500 

075 

018.0 

2644* 

2.500 

075 

036.0 

2645* 

2.500 

075 

054.0 

2646* 

2.500 

075 

072.0 

2647* 

2.500 

075 

090.0 

2648* 

2.500 

075 

108.0 

2649* 

2.500 

075 

126.0 

2650* 

2.500 

075 

144.0 

2651  * 

2.500 

075 

162.0 

2652* 

2.500 

075 

180.0 

2653* 

2.500 

075 

198.0 

2654* 

2.500 

075 

216.0 

2655* 

2.500 

075 

234.0 

2656* 

2.500 

075 

252.0 

2657* 

2.500 

075 

270.0 

2658* 

2.500 

075 

288.0 

2659* 

2.500 

075 

306.0 

2660* 

2.500 

075 

324.0 

2661  * 

2.500 

075 

342.0 

2662* 

2.500 

090 

000.0 

2663* 

2.500 

090 

015.0 

2664* 

2.500 

090 

030.0 

2665* 

2.500 

090 

045.0 

2666* 

2.500 

090 

060.0 

2667* 

2.500 

090 

075.0 
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*  2668  *  2.500  090  090.0 

*  2669  *  2.500  090  105.0 

*  2670  *  2.500  090  120.0 

*  2671  *  2.500  090  135.0 

*  2672  *  2.500  090  150.0 

*  2673  *  2.500  090  165.0 

*  2674  *  2.500  090  180.0 

*  2675  *  2.500  090  195.0 

*  2676  *  2.500  090  210.0 

*  2677  *  2.500  090  225.0 

*  2678  *  2.500  090  240.0 

*  2679  *  2.500  090  255.0 

*  2680  *  2.500  090  270.0 

*  2681  *  2.500  090  285.0 

*  2682  *  2.500  090  300.0 

*  2683  *  2.500  090  315.0 

*  2684  *  2.500  090  330.0 

*  2685  *  2.500  090  345.0 

*  2686  *  2.500  105  000.0 

*  2687  *  2.500  105  018.0 

*  2688  *  2.500  105  036.0 

*  2689  *  2.500  105  054.0 

*  2690  *  2.500  105  072.0 

*  2691  *  2.500  105  090.0 

*  2692  *  2.500  105  108.0 

*  2693  *  2.500  105  126.0 

*  2694  *  2.500  105  144.0 

*  2695  *  2.500  105  162.0 

*  2696  *  2.500  105  180.0 

*  2697  *  2.500  105  198.0 

*  2698  *  2.500  105  216.0 

*  2699  *  2.500  105  234.0 

*  2700  *  2.500  105  252.0 

*  2701  *  2.500  105  270.0 

*  2702*  2.500  105  288.0 

*  2703  *  2.500  105  306.0 

*  2704  *  2.500  105  324.0 

*  2705  *  2.500  105  342.0 

*  2706  *  2.500  120  000.0 

*  2707  *  2.500  120  022.5 

*  2708  *  2.500  120  045.0 

*  2709  *  2.500  120  067.5 

*  2710  *  2.500  120  090.0 

*  2711  *  2.500  120  112.5 

*  2712  *  2.500  120  135.0 

*  2713  *  2.500  120  157.5 

*  2714  *  2.500  120  180.0 

*  2715  *  2.500  120  202.5 

*  2716  *  2.500  120  225.0 

*  2717  *  2.500  120  247.5 

*  2718*  2.500  120  270.0 
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*  2719  *  2.500  120  292.5 

*  2720  *  2.500  120  315.0 

*  2721  *  2.500  120  337.5 

*  2722  *  2.500  135  000.0 

*  2723  *  2.500  135  030.0 

*  2724  *  2.500  135  060.0 

*  2725  *  2.500  135  090.0 

*  2726  *  2.500  135  120.0 

*  2727  *  2.500  135  150.0 

*  2728  *  2.500  135  180.0 

*  2729  *  2.500  135  210.0 

*  2730  *  2.500  135  240.0 

*  2731  *  2.500  135  270.0 

*  2732  *  2.500  135  300.0 

*  2733  *  2.500  135  330.0 

*  2734  *  2.500  150  000.0 

*  2735  *  2.500  150  045.0 

*  2736  *  2.500  150  090.0 

*  2737  *  2.500  150  135.0 

*  2738  *  2.500  150  180.0 

*  2739  *  2.500  150  225.0 

*  2740  *  2.500  150  270.0 

*  2741  *  2.500  150  315.0 

*  2742  *  2.500  165  000.0 

*  2743  *  2.500  165  090.0 

*  2744  *  2.500  165  180.0 

*  2745  *  2.500  165  270.0 

*  2746  *  2.500  180  000.0 
ELEMENTS 

SHEL06C  25CD4SH  2 
401  406  408  402  407  403 
401  408  410  403  409  404 
401  410412  404  411  405 
401  412  406  405  413  402 
406  426  428  414  427  415 
428  408  406  416  407  415 
408  428  430  416  429  417 
408  430  432  417  431  418 
432  410  408  419  409  418 
410  432  434  419  433  420 
410  434  436  420  435  421 
436  412  410  422  411  421 
412  436  438  422  437  423 
412  438  440  423  439  424 
440  406  412  425  413  424 
406  440  426  425  441  414 
426  462  464  442  463  443 
464  428  426  444  427  443 
428  464  466  444  465  445 
466  430  428  446  429  445 
430  466  468  446  467  447 
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430  468  470  447  469  448 
470  432  430  449  431  448 
432  470  472  449  471  460 
472  434  432  451  433  450 
434  472  474  451  473  452 
434  474  476  452  475  453 
476  436  434  454  435  453 
436  476  478  454  477  455 
478  438  436  456  437  455 
438  478  480  456  479  457 
438  480  482  457  481  458 
482  440  438  459  439  458 
440  482  484  459  483  460 
484  426  440  461  441  460 
426  484  462  461  485  442 
506  464  462  487  463  486 
464  506  508  487  507  488 
508  466  464  489  465  488 
466  508  510  489  509  490 
510  468  466  491  467  490 
510  470  468  492  469  491 
470  510  512492  511  493 
512  472  470  494  471  493 
472  512  514  494  513  495 
514  474  472  496  473  495 
514  476  474  497  475  496 
476  514  516  497  515  498 
516  478  476  499  477  498 
478  516  518  499  517  500 
518  480  478  501  479  500 
518  482  480  502  481  501 
482  518  520  502  519  503 
520  484  482  504  483  503 
484  520  506  504  521  505 
506  462  484  486  485  505 
534  508  506  523  507  52 2 
508  534  536  523  535  524 
536  510  508  525  509  524 
536  512  510  526  511  525 
512  536  538  526  537  527 
538  514  512  528  513  527 
538  516  514  529  515  528 
516  538  540  529  539  530 
540  518  516  531  517  530 
540  520  518  532  519  531 
520  540  534  532  541  533 
534  506  520  522  521  533 
546  536  534  543  535  542 
546  538  536  544  537  543 
546  540  538  545  539  544 
546  534  540  542  541  545 
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TRIA12I 

601  608  606  401  408  406  603  607  602  403  407  402 
601  610  608  401  410  408  604  609603  404  409  403 
601  612  610  401  412  410  605  611  604  405  411  404 
601  606  612  401  406  412  602  613  605  402  413  405 
606  628  626  406  428  426  615  627  614  415  427  414 
628  606  608  428  406  408  615  607  616  415  407  416 
608  630  628  408  430  428  617  629  616  417  429  416 
608  632  630  408  432  430  618  631  617  418  431  417 
632  608  610  432  408  410  618  609  619  418  409  419 
610  634  632  410  434  432  620  633  619  420  433  419 
610  636  634  410  436  434  621  635  620  421  435  420 
636  610  612  436  410412  621  611  622  421  411422 
612  638  636  412  438  436  623  637  622  423  437  422 
612  640  638  412  440  438  624  639  623  424  439  423 
640  612  606  440  412  406  624  613  625  424  413  425 
606  626  640  406  426  440  614  641  625  414  441  425 
626  664  662  426  464  462  643  663  642  443  463  442 
664  626  628  464  426  428  643  627  644  443  427  444 
628  666  664  428  466  464  645  665  644  445  465  444 
666  628  630  466  428  430  645  629  646  445  429  446 
630  668  666  430  468  466  647  667  646  447  467  446 
630  670  668  430  470  468  648  669  647  448  469  447 
670  630  632  470  430  432  648  631  649448431  449 
632  672  670  432  472  470  650  671  649  450  471  449 
672  632  634  472  432  434  650  633  651  450433  451 
634  674  672  434  474  472  652  673  651  452  473  451 
634  676  674  434  476  474  653  675  652  453  475  452 
676  634  636  476  434  436  653  635  654  453435  454 
636  678  676  436  478  476  655  677  654  455  477  454 
678  636  638  478  436  438  655  637  656  455  437  456 
638  680  678  438  480  478  657  679  656  457  479  456 
638  682  680  438  482  480  658  681  657  458  481  457 
682  638  640  482  438  440  658  639  659  458  439  459 
640  684  682  440  484  482  660  683  659  460  483  459 
684  640  626  484  440  426  660  641  661  460  441  461 
626  662  684  426  462  484  642  685  661  442  485  461 
706  662  664  506  462  464  686  663  687  486  463  487 
664  708  706  464  508  506  688  707  687  488  507  487 
708  664  666  508  464  466  688  665  689  488  465  489 
666  710  708  466  510  508  690  709  689  490  509  489 
710  666  668  510  466  468  690  667  691  490  467  491 
710  668  670  510  468  470  691  669  692  491  469  492 
670  712  710  470  512  510  693  711  692  493  511  492 
712  670  672  512  470  472  693  671  694  493  471  494 
672  714  712  472  514  512  695  713  694  495  513  494 
714  672  674  514  472  474  695  673  696  495  473  496 
71 4  674  676  51 4  474  476  696  675  697  496  475  497 
676  716  714  476  516  514  698  715  697  498  515  497 
716  676  678  516  476  478  698  677  699  498  477  499 
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678  718  716  478  518  516  700  717  699  500  517  499 
718  678  680  518  478  480  700  679  701  500  479  501 
718  680  682  518  480482  701  681  702501  481  502 
682  720  718  482  520  518  703  719  702  503  519  502 
720  682  684  520  482  484  703  683  704  503  483  504 
684  706  720  484  506  520  705  721  704  505  521  504 
706  684  662  506  484  462  705  685  686  505  485  486 
734  706  708  534  506  508  722  707  723  522  507  523 
708  736  734  508  536  534  724  735  723  524  535  523 
736  708  710  536  508  510  724  709  725  524  509  525 
736  710  712  536  510  512  725  71 1  726  525  51 1  526 
712  738  736  512  538  536  727  737  726  527  537  526 
738  712  714  538  512  514  727  713  728  527  513  528 
738  714  716  538  514  516  728  715  729  528  515  529 
716  740  738  516  540  538  730  739  729  530  539  529 
740  716  718  540  516  518  730  717  731  530  517  531 
740  718  720  540  518  520  731  719  732  531  519  532 
720  734  740  520  534  540  733  741  732  533  541  532 
734  720  706  534  520  506  733  721  722  533  521  522 
746  734  736  546  534  536  742  735  743  542  535  543 
746  736  738  546  536  538  743  737  744  543  537  544 
746  738  740  546  538  540  744  739  745  544  539  545 
746  740  734  546  540  534  745  741  742  545  541  542 

PRIS15F  EAU 

601  606  608  1001  1006  1008  602  607  603  801  806  808  1002  1007  1003 

601  608  610  1001  1008  1010  603  609  604  801  808  810  1003  1009  1004 

601  610  6121001  10101012  604  611  605  801  810  8121004  1011  1005 
601  612  606  1001  1012  1006  605  613  602  801  812  806  1005  1013  1002 

606  626  628  1006  1026  1028  614  627  615  806  826  828  1014  1027  1015 

628  608  606  1028  1008  1006  616  607  615  828  808  806  1016  1007  1015 

608  628  630  1008  1028  1030  616  629  617  808  828  830  1016  1029  1017 

608  630  632  1008  1030  1032  617  631  618  808  830  832  1017  1031  1018 

632  610  608  1032  1010  1008  619  609  618  832  810  808  1019  1009  1018 

610  632  634  1010  1032  1034  619  633  620  810  832  834  1019  1033  1020 

610  634  636  1010  1034  1036  620  635  621  810  834  836  1020  1035  1021 

636  612  610  1036  1012  1010  622  611  621  836  812  810  1022  1011  1021 

612  636  638  1012  1036  1038  622  637  623  812  836  838  1022  1037  1023 

612  638  640  1012  1038  1040  623  639  624  812  838  840  1023  1039  1024 

640  606  612  1040  1006  1012  625  613  624  840  806  812  1025  1013  1024 

606  640  626  1006  1040  1026  625  641  614  806  840  826  1025  1041  1014 

626  662  664  1026  1062  1064  642  663  643  826  862  864  1042  1063  1043 

664  628  626  1064  1028  1026  644  627  643  864  828  826  1044  1027  1043 

628  664  666  1028  1064  1066  644  665  645  828  864  866  1044  1065  1045 

666  630  628  1066  1030  1028  646  629  645  866  830  828  1046  1029  1045 

630  666  668  1030  1066  1068  646  667  647  830  866  868  1046  1067  1047 

630  668  670  1030  1068  1070  647  669  648  830  868  870  1047  1069  1048 

670  632  630107010321030  649  631  648  870  832  83010491031  1048 

632  670  672103210701072  649  671  650  832  870  872  10491071  1050 

672  634  632  1072  1034  1032  651  633  650  872  834  832  1051  1033  1050 

634  672  674103410721074  651  673  652  834  872  8741051  1073  1052 
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634  674  676  1034  1074  1076  652  675  653  834  874  876  1052  1075  1053 

676  636  634  1076  1036  1034  654  635  653  876  836  834  1054  1035  1053 

636  676  678  1036  1076  1078  654  677  655  836  876  878  1054  1077  1055 

678  638  636  1078  1038  1036  656  637  655  878  838  836  1056  1037  1055 

638  678  680  1038  1078  1080  656  679  657  838  878  880  1056  1079  1057 

638  680  682  1038  1080  1082  657  681  658  838  880  882  1057  1081  1058 

682  640  638  1082  1040  1038  659  639  658  882  840  838  1059  1039  1058 

640  682  684  1040  1082  1084  659  683  660  840  882  884  1059  1083  1060 

684  626  640  1084  1026  1040  661  641  660  884  826  840  1061  1041  1060 

626  684  662  1026  1084  1062  661  685  642  826  884  862  1061  1085  1042 

706  664  662  1106  1064  1062  687  663  686  906  864  862  1087  1063  1086 

664  706  708106411061108  687  707  688  864  906  908  108711071088 

708  666  664  1108  1066  1064  689  665  688  908  866  864  1089  1065  1088 

666  708  710106611081110  689  709  690  866  908  910108911091090 
710  668  666  1110  1068  1066  691  667  690  910  868  866  1091  1067  1090 

710  670  668111010701068  692  669  691  910  870  868  1092  1069  1091 

670  710  712  107011101112  692  711  693  870  910  912  10921111  1093 

712  672  6701112  10721070  694  671  693  912  872  8701094  1071  1093 

672  712  714  107211121114  694  713  695  872  912  914  109411131095 
714  674  672  1114  1074  1072  696  673  695  914  874  872  1096  1073  1095 

714  676  67411141076  1074  697  675  696  914  876  874  10971075  1096 

676  714  716107611141116  697  715  698  876  914  916  10971115  1098 

716  678  676111610781076  699  677  698  916  878  876  1099  10771098 

678  716  718  107811161118  699  717  700  878  916  918109911171100 

718  680  67811181080  1078  701  679  700  918  880  8781101  10791100 

718  682  680111810821080  702  681  701  918  882  88011021081  1101 

682  718  720108211181120  702  719  703  882  918  920110211191103 

720  684  6821120  10841082  704  683  703  920  884  882110410831103 

684  720  706  108411201106  704  721  705  884  920  90611041121  1105 

706  662  684110610621084  686  685  705  906  862  884  1086  10851105 

734  708  706  1134  1108  1106  723  707  722  934  908  906  1123  1107  1122 

708  734  736110811341136  723  735  724  908  934  936112311351124 
736  710  708113611101108  725  709  724  936  910  908112511091124 
736  712  710113611121110  726  711  725  936  912  91011261111  1125 
712  736  738111211361138  726  737  727  912  936  938112611371127 
738  714  712113811141112  728  713  727  938  914  912112811131127 
738  716  714113811161114  729  715  728  938  916  914112911151128 
716  738  740111611381140  729  739  730  916  938  940112911391130 
740  718  716114011181116  731  717  730  940  918  9161131  11171130 
740  720  718114011201118  732  719  731  940  920  918113211191131 
720  740  734112011401134  732  741  733  920  940  93411321141  1133 

734  706  7201134  11061120  722  721  733  934  906  92011221121  1133 

746  736  734114611361134  743  735  742  946  936  934114311351142 
746  738  736114611381136  744  737  743  946  938  936114411371143 
746  740  738114611401138  745  739  744  946  940  938114511391144 
746  734  740114611341140  742  741  745  946  934  94011421141  1145 
1001  1006  1008  1401  1406  1408  1002  1007  1003  1201  1206  1208  1402  1407  1403 

1001  1008  1010  1401  1408  1410  1003  1009  1004  1201  1208  1210  1403  1409  1404 

1001  1010  1012  1401  1410  1412  1004  1011  1005  1201  1210  1212  1404  1411  1405 

1001  1012  1006  1401  1412  1406  1005  1013  1002  1201  1212  1206  1405  1413  1402 

1006  1026  1028  1406  1426  1428  1014  1027  1015  1206 1226 1228  1414  1427  1415 
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1028  1008  1006  1428  1408  1406  1016  1007  1015  1228  1208  1206  1416  1407  1415 
1008  1028  1030  1408  1428  1430  1016  1029  1017  1208  1228  1230  1416  1429  1417 
1008  1030  1032  1408  1430  1432  1017  1031  1018  1208  1230  1232  1417  1431  1418 
1032  1010  1008  1432  1410  1408  1019  1009 1018  1232  1210  1208  1419  1409  1418 
1010  1032  1034  1410  1432  1434  1019  1033  1020  1210  1232  1234  1419  1433  1420 
1010  1034  1036  1410  1434  1436  1020  1035  1021  1210  1234  1236  1420  1435  1421 
1036  1012  1010  1436  1412  1410  1022  1011  1021  1236  1212  1210  1422  1411  1421 
1012  1036  1038  1412  1436  1438  1022  1037  1023  1212  1236  1238  1422  1437  1423 
1012  1038 1040 1412  1438  1440  1023  1039  1024  1212  1238  1240  1423  1439  1424 
1040  1006  1012  1440  1406  1412  1025  1013  1024  1240  1206  1212  1425  1413  1424 
1006  1040  1026  1406  1440  1426  1025  1041  1014  1206  1240  1226  1425  1441  1414 
1026  1062  1064  1426  1462  1464  1042  1063  1043  1226  1262  1264  1442  1463  1443 
1064  1028  1026  1464  1428  1426  1044  1027  1043  1264  1228  1226  1444  1427  1443 
1028  1064  1066  1428  1464  1466  1044  1065  1045  1228  1264  1266  1444  1465  1445 
1066  1030  1028  1466  1430  1428  1046  1029  1045  1266  1230  1228  1446  1429  1445 
1030  1066  1068  1430  1466  1468  1046  1067  1047  1230  1266  1268  1446  1467  1447 
1030  1068  1070 1430  1468  1470  1047  1069  1048  1230  1268  1270  1447  1469  1448 
1070  1032  1030  1470  1432  1430  1049  1031  1048  1270  1232  1230  1449  1431  1448 
1032  1070  1072  1432  1470  1472  1049  1071  1050  1232  1270  1272  1449  1471  1450 
1072  1034  1032 1472  1434  1432  1051  1033  1050  1272  1234  1232  1451  1433  1450 
1034  1072  1074 1434  1472  1474  1051  1073  1052  1234  1272  1274  1451  1473  1452 
1034  1074  1076  1434  1474  1476  1052  1075  1053  1234  1274  1276  1452  1475  1453 
1076  1036  1034  1476  1436  1434  1054  1035  1053  1276  1236  1234  1454  1435  1453 
1036  1076  1078 1436  1476  1478  1054  1077  1055  1236  1276  1278  1454  1477  1455 
1078  1038  1036  1478  1438  1436  1056  1037  1055  1278  1238  1236  1456  1437  1455 
1038  1078  1080  1438  1478  1480  1056  1079  1057  1238  1278  1280  1456  1479  1457 
1038  1080  1082  1438 1480  1482  1057  1081  1058 1238  1280 1282  1457  1481  1458 
1082  1040  1038  1482  1440  1438  1059  1039  1058  1282  1240  1238  1459  1439  1458 
1040  1082  1084  1440  1482  1484  1059  1083  1060 1240  1282  1284  1459  1483  1460 
1084  1026  1040  1484  1426  1440  1061  1041  1060  1284  1226  1240  1461  1441  1460 
1026  1084  1062  1426  1484  1462  1061  1085  1042  1226  1284  1262  1461  1485  1442 
1106  1064  1062  1506 1464  1462  1087  1063  1086  1306  1264  1262  1487  1463  1486 
1064  1106  1108  1464  1506  1508  1087  1107  1088  1264  1306  1308  1487  1507  1488 
1108  1066  1064  1508  1466  1464  1089  1065  1088  1308  1266  1264  1489  1465  1488 
1066  1108  1110  1466  1508  1510  1089  1109  1090  1266  1308  1310  1489  1509  1490 
1110  1068  1066  1510  1468  1466  1091  1067  1090  1310  1268  1266  1491  1467  1490 
1110  1070  1068  1510  1470  1468  1092  1069  1091  1310  1270  1268  1492  1469  1491 
1070  1110  1112  1470  1510  1512  1092  1111  1093  1270  1310  1312  1492  1511  1493 
1112  1072  1070  1512  1472  1470  1094  1071  1093  1312  1272  1270  1494  1471  1493 
1072  1112  1114  1472  1512  1514  1094  1113  1095  1272  1312  1314  1494  1513  1495 
1114  1074  1072  1514  1474  1472  1096  1073  1095  1314  1274  1272  1496  1473  1495 
1114  1076  1074  1514  1476  1474  1097  1075  1096  1314  1276  1274  1497  1475  1496 
1076  1114  1116  1476  1514  1516  1097  1115  1098  1276  1314  1316  1497  1515  1498 
1116  1078  1076  1516  1478  1476  1099  1077  1098  1316  1278  1276  1499  1477  1498 
1078  1116  1118  1478  1516  1518  1099  1117  1100  1278  1316  1318  1499  1517  1500 
1118  1080  1078  1518  1480  1478  1101  1079  1100  1318  1280  1278  1501  1479  1500 
1118  1082  1080  1518  1482  1480  1102  1081  1101  1318  1282  1280  1502  1481  1501 
1082  1118  1120  1482  1518  1520  1102  1119  1103  1282  1318  1320  1502  1519  1503 
1120  1084  1082  1520  1484  1482  1104  1083  1103  1320  1284  1282  1504  1483  1503 
1084  1120  1106  1484  1520  1506  1104  1121  1105  1284  1320  1306  1504  1521  1505 
1106  1062  1084  1506  1462  1484  1086  1085  1105  1306  1262  1284  1486  1485  1505 
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1426  1484  1462  1826  1884  1862  1461  1485 
1506  1464  1462  1906  1864  1862  1487  1463 
1464  1506  1508  1864  1906  1908  1487  1507 
1508  1466  1464  1908  1866  1864  1489  1465 
1466  1508  1510  1866  1908  1910  1489  1509 
1510  1468  1466  1910  1868  1866  1491  1467 
1510  1470  1468  1910  1870  1868  1492  1469 
1470  1510  1512  1870  1910  1912  1492  1511 
1512  1472  1470  1912  1872  1870  1494  1471 
1472  1512  1514  1872  1912  1914  1494  1513 
1514  1474  1472  1914  1874  1872  1496  1473 
1514  1476  1474  1914  1876  1874  1497  1475 
1476  1514  1516  1876  1914  1916  1497  1515 
1516  1478  1476  1916  1878  1876  1499  1477 
1478  1516  1518  1878  1916  1918  1499  1517 
1518  1480  1478  1918  1880  1878  1501  1479 
1518  1482  1480  1918  1882  1880  1502  1481 
1482  1518  1520  1882  1918  1920  1502  1519 
1520  1484 1482  1920 1884  1882  1504  1483 
1484  1520  1506  1884  1920  1906  1504  1521 
1506  1462  1484  1906  1862  1884  1486  1485 
1534  1508  1506  1934  1908  1906  1523  1507 
1508  1534  1536  1908  1934  1936  1523  1535 
1536  1510  1508  1936  1910  1908  1525  1509 
1536  1512  1510  1936  1912  1910  1526  1511 
1512  1536  1538  1912  1936  1938  1526  1537 
1538  1514  1512  1938 1914 1912  1528  1513 
1538  1516  1514  1938  1916  1914  1529 1515 
1516  1538  1540 1916  1938  1940  1529  1539 
1540  1518  1516  1940  1918  1916  1531  1517 
1540  1520  1518  1940  1920  1918  1532  1519 
1520  1540  1534  1920  1940  1934  1532  1541 
1534  1506  1520  1934  1906  1920  1522  1521 
1546  1536  1534  1946  1936  1934  1543  1535 
1546  1538  1536  1948  1938  1936  1544  1537 
1546  1540  1538  1946  1940  1938  1545  1539 
1546  1534  1540  1946  1934  1940  1542  1541 
1801  1806  1808  2201  2206  2208  1802  1807 
1801  1808  181022012208  221018031809 
1801  1810  1812  2201  2210  2212  1804  1811 
1801  1812  1806  2201  2212  2206  1805  1813 
1806  1826  1828  2206  2226  2228  1814  1827 
1828  1808  1806  2228  2208  2206  1816  1807 
1808  1828  1830  2208  2228  2230  1816  1829 
1808  1830  1832  2208  2230  2232  1817  1831 
1832  1810  1808  2232  2210  2208  1819  1809 
1810  1832  1834  2210  2232  2234  1819  1833 
1810  1834  1836  2210  2234  2236  1820  1835 
1836  1812  1810  2236  2212  2210  1822  1811 
1812  1836  1838  2212  2236  2238  1822  1837 
1812  1838  1840  2212  2238  2240  1823  1839 


1442  1626  1684  1662  1861  1885  1842 
1486  1706  1664  1662  1887  1863  1886 
1488  1864  1706  1708  1887  1907  1888 
1488 1708  1666  1664  1889  1865  1888 
1490  1666  1708  1710  1889  1909  1890 

1490  1710  1668  1666  1891  1867  1890 

1491  1710  1670  1668  1892  1869  1891 
1493  1670  1710  1712  1892  1911  1893 
1493  1712  1672  1670  1894  1871  1893 
1495  1672  1712  1714  1894  1913  1895 

1495  1714  1674  1672  1896  1873  1895 

1496  1714  1676  1674  1897  1875  1896 
1498  1676  1714 1716  1897  1915  1898 
1498  1716  1678  1676  1899  1877  1898 
1500  1678  1716  1718  1899  1917  1900 

1500  1718  1680  1678  1901  1879  1900 

1501  1718  1682  1680  1902  1881  1901 
1503  1682  1718  1720  1902  1919  1903 
1503  1720  1684  1682  1904  1883  1903 
1505  1684  1720  1706  1904  1921  1905 
1505  1706 1662  1684  1886  1885  1905 
1522 1734  1708  1706  1923  1907  1922 
1524  1708  1734  1736  1923  1935  1924 

1524  1736  1710  1708  1925  1909  1924 

1525  1736  1712  1710  1926  1911  1925 
1527  1712  1736  1738  1926  1937  1927 
1527  1738  1714  1712  1928  1913  1927 
1528 1738  1716  1714  1929  1915  1928 
1530 1716  1738  1740  1929 1939  1930 

1530  1740  1718  1716  1931  1917  1930 

1531  1740  1720  1718  1932  1919  1931 
1533  1720  1740  1734  1932  1941  1933 
1533  1734  1706  1720  1922  1921  1933 

1542  1746  1736  1734  1943  1935  1942 

1543  1746  1738  1736  1944  1937  1943 

1544  1746  1740  1738  1945  1939  1944 

1545  1746  1734  1740  1942  1941  1945 

1803  2001  2006  2008  2202  2207  2203 

1 804  2001  2008  2010  2203  2209  2204 

1805  2001  2010  2012  2204  2211  2205 
1802  2001  2012  2006  2205  2213  2202 
1815  2006  2026  2028  2214  2227  2215 
1815  2028  2008  2006  2216  2207  2215 

1817  2008  2028  2030  2216  2229  2217 

1818  2008  2030  2032  2217  2231  2218 
1818  2032  2010  2008  2219  2209  2218 

1820  2010  2032  2034  2219  2233  2220 

1821  2010  2034  2036  2220  2235  2221 
1821  2036  2012  2010  2222  2211  2221 

1823  2012  2036  2038  2222  2237  2223 

1824  2012  2038  2040  2223  2239  2224 
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1840  1806  1812  2240  2206  2212  1825  1813  1824  2040  2006  2012  2225  2213  2224 
1806  1840  1826  2206  2240  2226  1825  1841  1814  2006  2040  2026  2225  2241  2214 
1826  1862  1864  2228  2262  22641842186318432026  2062  2064  2242  2263  2243 
1864  1828  1826  2264  2228  2226  1844  1827  1843  2064  2028  2028  2244  2227  2243 
182818841866  22282284  22661844186518452028  2064  2066  2244  2265  2245 
1866  1830  1828  2266  2230  2228  1846  1829  1845  2066  2030  2028  2246  2229  2245 
1830  1866  1868  2230  2266  2268  1846  1867  1847  2030  2066  2068  2246  2267  2247 
1830  1868  1870  2230  2268  2270  1847  1869  1848  2030  2068  2070  2247  2269  2248 
1870  1832  1830  2270  2232  2230  1849  1831  1848  2070  2032  2030  2249  2231  2248 
1832  1870  1872  2232  2270  2272  1849  1871  1850  2032  2070  2072  2249  2271  2250 
1872  1834  1832  2272  2234  2232  1851  1833  1850  2072  2034  2032  2251  2233  2250 
1834  1872  1874  2234  2272  2274  1851  1873  1852  2034  2072  2074  2251  2273  2252 
1834  1874  1876  2234  2274  2276  1852  1875  1853  2034  2074  2076  2252  2275  2253 
1876  1836  1834  2276  2236  2234  1854  1835  1853  2076  2036  2034  2254  2235  2253 
1836  1876  1878  2236  2276  2278  1854  1877  1855  2038  2076  2078  2254  2277  2255 
1878  1838  1836  2278  2238  2238  1856  1837  1855  2078  2038  2036  2256  2237  2255 
1838  1878  1880  2238  2278  2280  1856  1879  1857  2038  2078  2080  2256  2279  2257 
1838  1880  1882  2238  2280  2282  1857  1881  1858  2038  2080  2082  2257  2281  2258 
1882  1840  1838  2282  2240  2238  1859  1839  1858  2082  2040  2038  2259  2239  2258 
1840  1882  1884  2240  2282  2284  1859  1883  1860  2040  2082  2084  2259  2283  2260 
1884  1826  1840  2284  2226  2240  1861  1841  1860  2084  2026  2040  2261  2241  2260 
1826188418622226  2284  22621861  18851842  2026  2084  2062  22612285  2242 
1906  1864  1862  2306  2264  2262  1887  1863 1886  2106  2064  2062  2287  2263  2286 
1864  1906  1908  2264  2306  2308  1887  1907  1888  2064  2106  2108  2287  2307  2288 
1908  1866  1864  2308  2266  2264  1889  1865  1888  2108  2066  2064  2289  2265  2288 
1866  1908  1910  2266  2308  2310  1889  1909 1890  2066  2108  2110  2289  2309  2290 
1910  1868  1866  2310  2268  2266  1891  1887  1890  21102068  2066  2291  2267  2290 
1910  1870  1868  2310  2270  2268  1892  1869  1891  2110  2070  2068  2292  2269  2291 
1870  1910  1912  2270  231023121892  1911  189320702110  21122292  2311  2293 
1912  1872  1870  2312  2272  2270  1894  1871  1893  2112  2072  2070  2294  2271  2293 
1872  1912  1914  2272  2312  2314  1894  1913  1895  2072  2112  2114  2294  2313  2295 
1914  1874  1872  2314  2274  2272  1896  1873  1895  2114  2074  2072  2296  2273  2295 
1914  1876  1874  2314  2276  2274  1897  1875  1U06  2114  2076  2074  2297  2275  2296 
1876  1914  1916  2276  2314  2316  1897  1915  1898  2076  2114  2116  2297  2315  2298 
1916  1878  1876  2316  2278  2276  1899  1877  1898  2116  2078  2076  2299  2277  2298 
1878  1916  1918  2278  2316  2318  1899  1917  1900  2078  21 16  2118  2299  2317  2300 
1918  1880  1878  2318  2280  2278  1901  1879  1900  2118  2080  2078  2301  2279  2300 
1918  1882  1880  2318  2282  2280  1902  1881  1901  2118  2082  2080  2302  2281  2301 
1882  1918  1920  2282  2318  2320  1902  1919  1903  2082  2118  2120  2302  2319  2303 
1920  1884  1882  2320  2284  2282  1904  1883  1903  2120  2084  2082  2304  2283  2303 
1884  1920  1906  2284  2320  2306  1904  1921  1905  2084  2120  2106  2304  2321  2305 
1906  1862  1884  2306  2262  2284  1886  1885  1905  2106  2062  2084  2286  2285  2305 
1934  1908  1906  2334  2308  2306  1923  1907  1922  2134  2108  2106  2323  2307  2322 
1908  1934  1936  2308  2334  2336  1923  1935  1924  2108  2134  2136  2323  2335  2324 
1936  1910  1908  2336  2310  2308  1925  1909  1924  2136  2110  2108  2325  2309  2324 
1936  1912  1910  2336  2312  2310  1926  1911  1925  2136  2112  21102326  2311  2325 
1912  1936  1938  2312  2336  2338  1926  1937  1927  2112  2136  2138  2326  2337  2327 
1938  1914  1912  2338  2314  2312  1928  19131927  2138  2114  2112  2328  2313  2327 
1938  1916  1914  2338  2316  2314  1929  1915  1928  2138  2116  2114  2329  2315  2328 
1916  1938  1940  2316  2338  2340  1929  1939  1930  2116  2138  2140  2329  2339  2330 
1940  1918  1916  2340  2318  2316  1931  1917  19302140  2118  2116  2331  2317  2330 
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1940  1920  1918  2340  2320  2318  1932  1919  1931  2140  2120  21 18  2332  2319  2331 
1920  1940  1934  2320  2340  2334  19321941  19332120  2140  21342332  2341  2333 
1934  1906  1920  2334  2306  2320  1922  1921  1933  2134  2106  2120  2322  2321  2333 
1946  1936  1934  2346  2336  2334  1943  1935  1942  2146  2136  2134  2343  2335  2342 
1946  1938  1936  2346  2338  2336  1944  1937  1943  2146  2138  2136  2344  2337  2343 
1946  1940  1938  2346  2340  2338  1945  1939  1944  2146  2140  2138  2345  2339  2344 
1946  1934  1940  2346  2334  2340  1942  1941  1945  2146  2134  2140  2342  2341  2345 
2201  2206  2208  2601  2606  2608  2202  2207  2203  2401  2406  2408  2602  2607  2603 
2201  2208  2210  2601  2608  2610  2203  2209  2204  2401  2408  2410  2603  2609  2604 
2201  2210  2212  2601  2610  2612  2204  2211  2205  2401  2410  2412  2604  261 1  2605 
2201  2212  2206  2601  2612  2606  2205  2213  2202  2401  2412  2406  2605  2613  2602 
2206  2226  2228  2606  2626  2628  2214  2227  2215  2406  2426  2428  2614  2627  2615 
2228  2208  2206  2628  2608  2606  2216  2207  2215  2428  2408  2406  2616  2607  2615 
2208  2228  2230  2608  2628  2630  2216  2229  2217  2408  2428  2430  2616  2629  2617 
2208  2230  2232  2608  2630  2632  2217  2231  2218  2408  2430  2432  2617  2631  2618 
2232  2210  2208  2632  2610  2608  2219  2209  2218  2432  2410  2408  2619  2609  2616 
2210  2232  2234  2610  2632  2634  2219  2233  2220  2410  2432  2434  2619  2633  2620 
2210  2234  2236  2610  2634  2636  2220  2235  2221  2410  2434  2436  2620  2635  2621 
2236  2212  2210  2636  2612  2610  2222  2211  2221  2436  2412  24102622  2611  2621 
2212  2236  2238  2612  2636  2638  2222  2237  2223  2412  2436  2438  2622  2637  2623 
2212  2238  2240  2612  2638  2640  2223  2239  2224  2412  2438  2440  2623  2639  2624 
2240  2206  2212  2640  2606  2612  2225  2213  2224  2440  2406  2412  2625  2613  2624 
2206  2240  2226  2606  2640  2626  2225  2241  2214  2406  2440  2426  2625  2641  2614 
2226  2262  2264  2626  2662  2664  2242  2263  2243  2426  2462  2464  2642  2663  2643 
2264  2228  2226  2664  2628  2626  2244  2227  2243  2464  2428  2426  2644  2627  2643 
2228  2264  2266  2628  2664  2666  2244  2265  2245  2428  2464  2466  2644  2665  2645 
2266  2230  2228  2666  2630  2628  2246  2229  2245  2466  2430  2428  2646  2629  2645 
2230  2266  2268  2630  2666  2668  2246  2267  2247  2430  2466  2468  2646  2667  2647 
2230  2268  2270  2630  2668  2670  2247  2269  2248  2430  2468  2470  2647  2669  2648 
2270  2232  2230  2670  2632  2630  2249  2231  2248  2470  2432  2430  2649  2631  2648 
2232  2270  2272  2632  2670  2672  2249  2271  2250  2432  2470  2472  2649  2671  2650 
2272  2234  2232  2672  2634  2632  2251  2233  2250  2472  2434  2432  2651  2633  2650 
2234  2272  2274  2634  2672  2674  2251  2273  2252  2434  2472  2474  2651  2673  2652 
2234  2 7* '  2276  2634  2674  2676  2252  2275  2253  2434  2474  2476  2652  2675  2653 
2276  2236  2234  2676  2636  2634  2254  2235  2253  2476  2436  2434  2654  2635  2653 
2236  2276  2278  2636  2676  2678  2254  2277  2255  2436  2476  2478  2654  2677  2655 
2278  2238  2236  2678  2638  2636  2256  2237  2255  2478  2438  2436  2656  2637  2655 
2238  2278  2280  2638  2678  2680  2256  2279  2257  2438  2478  2480  2656  2679  2657 
2238  2280  2282  2638  2680  2682  2257  2281  2258  2438  2480  2482  2657  2681  2658 
2282  2240  2238  2682  2640  2638  2259  2239  2258  2482  2440  2438  2659  2639  2658 
2240  2282  2284  2640  2682  2684  2259  2283  2260  2440  2482  2484  2659  2683  2660 
2284  2226  2240  2684  2626  2640  2261  2241  2260  2484  2426  2440  2661  2641  2660 
2226  2284  2262  2626  2684  2662  2261  2285  2242  2426  2484  2462  2661  2685  2642 
2306  2264  2262  2706  2664  2662  2287  2263  2286  2506  2464  2462  2687  2663  2686 
2264  2306  2308  2664  2706  2708  2287  2307  2288  2464  2506  2508  2687  2707  2688 
2308  2266  2264  2708  2666  2664  2289  2265  2288  2508  2466  2464  2689  2665  2688 
2266  2308  2310  2666  2708  2710  2289  2309  2290  2466  2508  2510  2689  2709  2690 
2310  2268  2266  2710  2668  2666  2291  2267  2290  2510  2468  2466  2691  2667  2690 
2310  2270  2268  2710  2670  2668  2292  2269  2291  2510  2470  2468  2692  2669  2691 
2270  2310  2312  2670  2710  2712  2292  231 1  2293  2470  2510  2512  2692  271 1  2693 
2312  2272  2270  2712  2672  2670  2294  2271  2293  2512  2472  2470  2694  2671  2693 
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2272  2312  2314  2672  2712  2714  2294  2313  2295  2472  2512  2514  2694  2713  2695 
2314  2274  2272  2714  2674  2672  2296  2273  2295  2514  2474  2472  2696  2673  2695 
2314  2276  2274  2714  2676  2674  2297  2275  2296  2514  2476  2474  2697  2675  2696 
2276  2314  2316  2676  2714  2716  2297  2315  2298  2476  2514  2516  2697  2715  2698 
2316  2278  2276  2716  2678  2676  2299  2277  2298  2516  2478  2476  2699  2677  2698 
2278  2316  2318  2678  2716  2718  2299  2317  2300  2478  2516  2518  2699  2717  2700 
2318  2280  2278  2718  2680  2678  2301  2279  2300  2518  2480  2478  2701  2679  2700 
2318  2282  2280  2718  2682  2680  2302  2281  2301  2518  2482  2480  2702  2681  2701 
2282  2318  2320  2682  2718  2720  2302  2319  2303  2482  2518  2520  2702  2719  2703 
2320  2284  2282  2720  2684  2682  2304  2283  2303  2520  2484  2482  2704  2683  2703 
2284  2320  2306  2684  2720  2706  2304  2321  2305  2484  2520  2506  2704  2721  2705 
2306  2262  2284  2706  2662  2684  2286  2285  2305  2506  2462  2484  2686  2685  2705 
2334  2308  2306  2734  2708  2706  2323  2307  2322  2534  2508  2506  2723  2707  2722 
2308  2334  2336  2708  2734  2736  2323  2335  2324  2508  2534  2536  2723  2735  2724 
2336  2310  2308  2736  2710  2708  2325  2309  2324  2536  2510  2508  2725  2709  2724 
2336  2312  2310  2736  2712  2710  2326  231 1  2325  2536  2512  2510  2726  271 1  2725 
2312  2336  2338  2712  2736  2738  2326  2337  2327  2512  2536  2538  2726  2737  2727 
2338  2314  2312  2738  2714  2712  2328  2313  2327  2538  2514  2512  2728  2713  2727 
2338  2316  2314  2738  2716  2714  2329  2315  2328  2538  2516  2514  2729  2715  2728 
2316  2338  2340  2716  2738  2740  2329  2339  2330  2516  2538  2540  2729  2739  2730 
2340  2318  2316  2740  2718  2716  2331  2317  2330  2540  2518  2516  2731  2717  2730 
2340  2320  2318  2740  2720  2718  2332  2319  2331  2540  2520  2518  2732  2719  2731 
2320  2340  2334  2720  2740  2734  2332  2341  2333  2520  2540  2534  2732  2741  2733 
2334  2306  2320  2734  2706  2720  2322  2321  2333  2534  2506  2520  2722  2721  2733 
2346  2336  2334  2746  2736  2734  2343  2335  2342  2546  2536  2534  2743  2735  2742 
2346  2338  2336  2746  2738  2736  2344  2337  2343  2546  2538  2536  2744  2737  2743 
2346  2340  2338  2746  2740  2738  2345  2339  2344  2546  2540  2538  2745  2739  2744 
2346  2334  2340  2746  2734  2740  2342  2341  2345  2546  2534  2540  2742  2741  2745 

TRIA06R  EAU  1 

2601  2606  2608  2602  2607  2603 
2601  2608  2610  2603  2609  2604 
2601  2610  2612  2604  261 1  2605 
2601  2612  2606  2605  2613  2602 
2606  2626  2628  2614  2627  2615 
2628  2608  2606  261 6  2607  261 5 
2608  2628  2630  2616  2629  2617 
2608  2630  2632  2617  2631  2618 
2632  2610  2608  2619  2609  2618 
2610  2632  2634  2619  2633  2620 
261 C  2634  2636  2620  2635  2621 
2636  2612  2610  2622  261 1  2621 
2612  2636  2638  2622  2637  2623 
2612  2638  2640  2623  2639  2624 
2640  2606  2612  2625  2613  2624 
2606  2640  2626  2625  2641  2614 
2626  2662  2664  2642  2663  2643 
2664  2628  2626  2644  2627  2643 
2628  2664  2666  2644  2665  2645 
2666  2630  2628  2646  2629  2645 
2630  2666  2668  2646  2667  2647 
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2630  2668  2670  2647  2669  2648 
2670  2632  2630  2649  2631  2648 
2632  2670  2672  2649  2671  2650 
2672  2634  2632  2651  2633  2650 
2634  2672  2674  2651  2673  2652 
2634  2674  2676  2652  2675  2653 
2676  2636  2634  2654  2635  2653 
2636  2676  2678  2654  2677  2655 
2678  2638  2636  2656  2637  2655 
2638  2678  2680  2656  2679  2657 
2638  2680  2682  2657  2681  2658 
2682  2640  2638  2659  2639  2658 
2640  2682  2684  2659  2683  2660 
2684  2626  2640  2661  2641  2660 
2626  2684  2662  2661  2685  2642 
2706  2664  2662  2687  2663  2686 
2664  2706  2708  2687  2707  2688 
2708  2666  2664  2689  2665  2688 
2666  2708  2710  2689  2709  2690 
2710  2668  2666  2691  2667  2690 
2710  2670  2668  2692  2669  2691 
2670  2710  2712  2692  271 1  2693 
2712  2672  2670  2694  2671  2693 
2672  2712  2714  2694  2713  2695 
2714  2674  2672  2696  2673  2695 
2714  2676  2674  2697  2675  2696 
2676  2714  2716  2697  2715  2698 
2716  2678  2676  2699  2677  2698 
2678  2716  2718  2699  2717  2700 
2718  2680  2678  2701  2679  2700 
2718  2682  2680  2702  2681  2701 
2682  2718  2720  2702  2719  2703 
2720  2684  2682  2704  2683  2703 
2684  2720  2706  2704  2721  2705 
2706  2662  2684  2686  2685  2705 
2734  2708  2706  2723  2707  2722 
2708  2734  2736  2723  2735  2724 
2736  2710  2708  2725  2709  2724 
2736  2712  2710  2726  271 1  2725 
2712  2736  2738  2726  2737  2727 
2738  2714  2712  2728  2713  2727 
2738  2716  2714  2729  2715  2728 
2716  2738  2740  2729  2739  2730 
2740  2718  2716  2731  2717  2730 
2740  2720  2718  2732  2719  2731 
2720  2740  2734  2732  2741  2733 
2734  2706  2720  2722  2721  2733 
2746  2736  2734  2743  2735  2742 
2746  2738  2736  2744  2737  2743 
2746  2740  2738  2745  2739  2744 
2746  2734  2740  2742  2741  2745 


ENO 
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APPENDIX  C 


PROGRAM  BNDMAT 
CALL  SVDSUB 
STOP 
END 


*•••*•*** 


SUBROUTINE  SVDSUB 

Q  *#*•*••**••**••***•*•««•**••••#•••*•••***••••*••*****•**•• 

PARAMETER  (NHARM=2,MSPHER=3, 

1  IBLOCK=(NHARM+1  )**2,NPOINT =146) 

IMPLICIT  REAL*8  (A-H.O-Z) 

COMPLEX'1 6  CZERO,CONE,CIIHRADB(0:NHARM),MATINF(NPOINT1IBLOCK), 

1  WK(NPOINT)lSIG(IBLOCK),EXTRA(IBLOCK),VM(IBLOCK,IBLOCK), 

2  UM(NPOINTtNPOINT),UIINV(IBLOCK,NPOINT),PHIFAC,PDATA(NPOINT)1 

3  ACHK(IBLOCK),RES(NPOINT),MATCHK(NPOINT1IBLOCK),HRADR(0:NHARM), 
4CTEM1.ZZCHK 

REAL'S  KK,KA,PMX(-NHARM:NHARM,0:NHARM),KR 
COMMON/COMO/CZERO,CONE.CI,ORIGIN(0:MSPHER,3) 

PM.ODO'DATAN(I.ODO) 

C 

C  THIS  ROUTINE  INITIALIZES  ARRAYS  FOR  THE  MAIN  ROUTINE 
C  SPECIFICALLY,  IT: 

C  SETS  COMPLEX  CONSTANTS  FOR  0,  1 ,  AND  I 

C  COMPUTES  HANKEL  &  BESSEL  FUNCTIONS  (&  DERIVATIVES)  FOR  RADIUS 
C 

C  INPUT: 

C  KK  =  WAVENUMBER 
C  RADSPH  =  RADIUS  OF  SPHERE 
C  CFLD=FLUID  WAVE  SPEED  1524  METERS/SEC 
C  FLDDEN=FLUID  DENSITY  1 000  KG/METER**3 
C  RHOC  =  RHO'C,  RADIATION  IMPEDANCE  OF  MEDIUM 
C  X,Y,Z  =  COORDINATES  FOR  EACH  SPHERE 
C  NPOINT  =  NUMBER  OF  DATA  POINTS  USED  IN  THE  CALCULATION 
C 

C  SET  COMPLEX  CONSTANTS 
CZERO  =  DCMPLX(O.ODO.O.ODO) 

CONE  =  DCMPLX(1.0D0,0.0D0) 

Cl  =  DCMPLX(0.0D0, 1 .0D0) 

C  INPUT  THE  EXCITATION  MODE  MMI  AND  NNI 
MMI=0 
NNI=2 

lll=NNI*NNI+NNI+MMI+1 
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C  INPUT  WAVENUMBER,  SPHERE  RADIUS.  AND  MEDIUM'  S  DENSITY*SOUND  SPEED 
RADSPH»0.5D0 
CFLDsl 492.94D0 
HERTZ=474.0D0 
KK  =  2*PI*HERTZ/CFLD 
FLDDENsI  OOO.ODO 
RHOC=CFLD*FLDDEN 
KA  =  KK*RADSPH 
WRITE(6,21 )  KK.RADSPH.RHOC 
21  FORMAT(5X,’  K=\1E10.4,'  A='  ,E10.4,*  RHO*C  ='  .E10.4) 

CALL  HANKEL(KA,NHARM,HRADB) 


C  THE  DATA  FILE  MUST  RESIDE  IN  UNIT  12  TO  DO  THE  FOLLOWING 
DO  10  Isrl.NPOINT 
READ(12,15)  X,Y,Z,PDATA(I) 

1 5  FORMAT(E1 1 .3.E1 2.3, El  2.3, El  3.3, El  2.3) 

IF(I.EQ.I)  THEN 
RBND=DSQRT  (X**2+Y**2+Z**2) 

KR=KK*RBND 
PRINT  *,KR 

CALL  HANKEL(KR.NHARM.HRADR) 

ENDIF 

THETX=DATAN2(DSQRT(X**2+Y**2),Z) 

CALL  LEGNDR(THETX,NHARM,PMX) 

PHIX=DATAN2(Y,X) 

THE  FOLLOWING  IS  TO  CHECK  OUTPUTS  FROM  ATILA  N=1,  M=-1 

ZZCHK=DCMPLX(  0.97307D0,  0.23051  D0)*DCONJG(HRADR(0)) 
ZZCHK=DCMPLX(0.98795D0,  0.1 5478DO)*DCONJG(HRADR(NNI))* 
1  PMX(MMI,NNI)*DCMPLX(DCOS(MMPPHIX),DSIN(MMI*PHIX)) 
ZZCHK=DCMPLX(0.99876D0,-0.049693D0)*DCONJG(HRADR(NNI))* 

1  PMX(MMI,NNI)*DCMPLX(DCOS(MMI*PHIX),DSIN(MMI*PHIX)) 


PDATA(I)=ZZCHK 

PHASED=DATAN2(DIMAG(PDATA(I)),DREAL(PDATA(I))) 
PHASEA=DATAN2(DIMAG(ZZCHK),DREAL(ZZCHK)) 

WRITE(6,16)  X,Y,Z,PDATA(I),CDABS(PDATA(I)), 

1  ZZCHK.CDAPS(ZZCHK) 

ZZCHK=ZZCHK-PDATA(I) 

WRITE(6,617)  ZZCHK.CDABS(ZZCHK) 

617  FORMAT(1 3X.2E1 7.5, El  7.5) 

WRITE(6,616)  PHASED, PHASEA.PHASEA-PHASED 
616  FORMAT(7X,3E17.6) 

16  FORMAT(2X,6E12.3y.38X,3E12.3) 

WRITE(25,125)  PHIX, PHASED 
125  FORMAT(2X,2F13.5) 

DO  10  N=0,NHARM 
DO  10  M=-N,N 

PHIFAC=DCMPLX(DCOS(M*PHIX),DSIN(M*PHIX)) 
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JJaN*N+N+M+1 

MATINF(l,JJ)«DCONJG(HRADR(N))*PMX(M,N)*PHIFAC 

MATCHK(I,JJ)=MAT1NF(I,JJ) 

10  CONTINUE 

C4LLZSVDC(MATINF, NPOINT, NPOINT, IBLOCK, SIG, EXTRA, UM, NPOINT, VM, 

1  IBLOCK.WK.21  .INFOX) 

WRITE(6,30)  (SIG(I),  1=1  .IBLOCK) 

30  FORMAT ( 1  X.2E20. 5) 

DO  40  M,  IBLOCK 
DO  40  J=1  .NPOINT 
UIINV(l,J)=CZERO 
DO  40  K*1  .IBLOCK 

IF(CDABS(SIG(K)).LT.10E-8)  GO  TO  40 
UIINV(l,J)=UIINV(l,J)+VM(l,K)*DCONJG(UM{J,K))/SIG(K) 

40  CONTINUE 
DO  50  1=1,  IBLOCK 
ACHK(l)=CZERO 
DO  50  K=1  .NPOINT 

ACHK(I)=ACHK(I)+UIINV(I,K)*PDATA(K) 

50  CONTINUE 
WRITE(6.160) 

160  FORMAT(1X,’  THE  FOLLOWING  ARE  THE  UNMODIFIED  T  MATRIX  ELEMENTS’ ,/, 

1  IX,’  OR  THE  VALUES  OF  THE  RADIATED  PRESSURES’ ) 

DO  55  N=0,NHARM 
DO  55  M=-N,N 
l=N*N+N+M+1 
CTEM1  =ACHK(I) 

55  WRITE(6,60)  l,CTEM1,CDABS(CTEM1) 

60  FORMAT(1  X.13.1  X.3E20.5) 

WRITE(6,161) 

161  FORMAT(1X,’  THE  FOLLOWING  ARE  THE  TRUE  T  MATRIX  ELEMENTS’ ) 

C  III  IS  THE  EXCITATION  MODE  EG  FOR  M=0,N=0 111=1 

DO  155  N=0,NHARM 
DO  155  M=-N,N 
l=N*N+N+M+1 

IF((N.EQ.NNI).AND.(M.EQ.MMI))  THEN 
CTEM1  =0.5D0*(ACHK(I)-1 .000) 

ELSE 

CTEM1=ACHK(I) 

ENDIF 

155  WRITE(6,60)  l,CTEM1,CDABS(CTEM1) 

C  WANT  TO  GET  A  HANDLE  ON  THE  RELATIVE  ERROR  INVOLVED  IN  THE  SVD  CALC 
DO  70  J=1  .NPOINT 
RES(J)=DCMPLX(0.0D0,0.0D0) 

DO  80  1=1  .IBLOCK  80  RES(J)=RES(J)+MATCHK(J,I)*ACHK(I) 
RES(J)=RES(J)-PDATA(J) 

70  CONTINUE 
DO  100  1=1  .NPOINT 
RESM=RESM+RES(l)*DCONJG(RES(l)) 

DATAM=DATAM+PDATA(l)*CONJG(PDATA(l» 

100  CONTINUE 
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RELERR-RESM/DATAM 
WRITE(6,110)  RESM.DAT AM, RELERR 
110  FORMAT(1X,'  MAGNITUDE  OF  RESIDUAL  =  ’ ,  E15A/.1X, 

1  1  MAGNITUDE  OF  DATA  =  ' ,  El  5.6,/, IX, 

2  ’  RELATIVE  ERROR  IN  THE  CALCULATION  =  * ,  El 5.6) 
RETURN 

END 


SUBROUTINE  HANKEL(X,NMAX,H) 

IMPLICIT  REAL*8  (A-H.O-Z) 

COMPLEX*16  H(0:NMAX) 

C 

C  GIVEN  THE  VARIABLE  X,  AND  THE  MAXIMUM  ORDER  NMAX, 

C  THIS  ROUTINE  GENERATES  THE  SPHERICAL  HANKEL  FUNCTION  HN 
C  FOR  ALL  N  FROM  0  TO  NMAX  (INCLUSIVE) 

C  INPUT: 

C  X  =  DOUBLE  PREC.  VARIABLE  (RADIUS) 

C  NMAX  *  INTEGER  MAXIMUM  ORDER  OF  BESSEL  FUNCTIONS  DESIRED 
C  OUTPUT: 

C  H(N)  =  ARRAY  OF  SPHERICAL  HANKEL  FUNCTIONS  HN(X),  WHERE 
C  HN  =  JN  +  I  YN 
C 

C  THIS  ROUTINE  IS  BASED  ON  THE  RECURSION  FORMULAE 
C  FROM  ABRAMOWITZ  &  STEGUN:  10.1.10  &  10.1.15,  PP.438-9 
C  THE  F  S  ARE  THE  COEFFICIENTS  OF  ORDER  N  &  -(N+1 ), 

C  THE  FO*  S  ARE  OLD  F  S,  FOR  RECURSION 
C 

IF  (  X  .LE.  0.0D0 )  THEN 
H(0)  =  DCMPLX(1.0D0,-1.0D35) 

DO  2  N  =  1 ,  NMAX 
H(N)  =  CMPLX(0.0D0,-1 .0D35) 

2  CONTINUE 
RETURN 
END  IF 
SX  =  DSIN(X) 

CX  =  DCOS(X) 

XINV  =  1 .0D0/X 
MIN  =  -1.0D0 
FN  =  XINV 
FMN  =  0.0D0 
FNO  =  FMN 
FMNO  =  FN 
DO  4  N  =  0,  NMAX 

H(N)  =  CMPLX(  FN*SX  +  M1N*FMN*CX,  -FN'CX  +  M1N*FMN*SX ) 

T1  =(2*N+1)*XINV 
T2  =  T1*FN  -  FNO 
FNO  =  FN 
FN  =  T2 

T2  =  -T1*FMN  -  FMNO 
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FMNO  -  FMN 
FMN-T2 
MIN* -MIN 
4  CONTINUE 
RETURN 
END 


SUBROUTINE  LEGNDR(THETA,NMAX,PMN) 

IMPLICIT  REAL*8  (A-H.O-Z) 

REAL‘8  PMN(-NMAX:NMAX,0:NMAX) 

C 

C  GIVEN  THE  VARIABLE  THETA,  AND  THE  MAXIMUM  ORDER  NMAX, 

C  THIS  ROUTINE  GENERATES  THE  ASSOC.  LEGENDRE  FUNCTIONS  PMN 
C  OF  THE  ARGUMENT  COSfTHETA)  (THETA  MUST  BE  BETWEEN  0  &  PI) 

C  FOR  ALL  N  FROM  0  TO  NMAX  (INCLUSIVE) 

C  AND  FOR  ALL  M  FROM  -N  TO  N  (SOME  OTHERS  SET  TO  ZERO) 

C  INPUT: 

C  THETA  =  VARIABLE  (POLAR  ANGLE),  MUST  BE  BETWEEN  0  &  PI  (INCL) 

C  NMAX  *  INTEGER  MAXIMUM  ORDER  OF  LEGENDRE  FUNCTIONS  DESIRED 
C  OUTPUT: 

C  PMN  *  DOUBLE  PREC.  ARRAY,  CONTAINS  ASSOC.  LEGENDRE  FNS 
C 

C  THIS  ROUTINE  IS  BASED  ON  THE  RECURSION  FORMULAE 
C  FROM  ABRAMOWITZ  &  STEGUN 
C 

X  »  DCOS(THETA) 

SINTHT  *  DSINfTHETA) 

IF  ( SINTHT  .GT.  0. )  THEN 
SININV  *  1.0D0/SINTHT 
ELSE 

SININV  =  0.0D0 
END  IF 

C  SET  VALUES  FOR  N  =  0, 1  (NMAX  MUST  BE  AT  LEAST  1) 

PMN(O.O)  *  1.0D0 
PMN(1 ,0)  *  0.0D0 
PMN(-1 ,0)  *  0.0D0 
PMN(0,1)  *  X 
PMN(1,1)  = -SINTHT 
PMN(-1,1)  *  SINTHT0.5D0 

C  IN  LOOP,  TNP1  *  2*N+1,  TNP2FC  =  (2*N+2)I,  MIN  =  (-1)**(N+1) 

TNP1  =  1.0D0 
TNP2FC  =  2.0D0 
MIN  =  -1.000 
DO  4  N  =  1,  NMAX-1 
TNP1  =TNP1  +2.0D0 
TNP2FC  *  TNP2FC  *  TNP1  *  (TNP1+1) 

MIN  =  -MIN 
DO  3  M  =  -N,  N 

PMN(M,N+1)  =  (TNP1*X*PMN(M,N)  -  (N+M)*PMN(M,N-1))/(N-M+1) 
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3  CONTINUE 
PMN(N+1,N)- 0.000 
PMN(-N-1  ,N) «  0.000 

PMN(N+1,N+1)  * (X*PMN(N,N+1)  -  TNP1*PMN(N,N))  *  SININV 
PMN(-N-1,N+1)  -  M1N*PMN(N+1,N+1)/TNP2FC 

4  CONTINUE 

C  DO  120  NaO.NMAX 

C  DO  120  M=-N,N 

Cl  20  WR!TE(6,130)  N.M,PMN(M.N) 

Cl 30  FORMATOX,1  N=\I4,1X,’  M^’.W.IX,’  PMN»*,F13.6) 
RETURN 
END 


l 
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